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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. They 
should be accompanied by a summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions. The summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed zoo words 
in . Authors are requested to submit MSS. in duplicate. These should be 

double spacing and leaving a margin of not less than one inch on the 
left-hand side. Corrections to the MSS. should be made in the text and not in the 
Unless a paper reaches the Secretaries more than seven days before a Council 
meeting it will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zs. f. Astrophys., 15, 239, 1938. 
H. Jeffreys, Theory of Probability, 2nd edn., section 5.45, p. 258, Oxford, 1948. 


3. Notation.—For technical astronomical terms, authors should conform closely to the 
recommendations of Commission 3 of the International Astronomical Union (Trans. 
I.A.U., Vol. VI, p. 345, 1938). Council has decided to adopt the I.A.U. 4-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, p. 221, 1932). In general 
matters, authors should follow the recommendations in Symbols, Signs and Abbreviations 
(London: Royal Society, 1951) except where these conflict with I.A.U. practice. 

4. Diagrams.—These should be designed to appear upright on the page, drawn 
about twice the size required in print and prepared for direct photographic 
reproduction except for the lettering, which should be inserted in pencil. 

should be given in the manuscript indicating where in the text the 
figure should appear. Blocks are retained by the Society for 10 years; unless the author 
requires them before the end of this period they are then destroyed. 


5. Tables.—These should be arranged so that they can be printed upright on 
the page. 

6. Proofs.—Costs of alterations exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximum conciseness and simplifi- 
cation of symbols. and equations consistent with clarity. 

7. Revised Manuscripts.—When papers are submitted in revised form it is especially 
requested that they be accompanied by the original MS. 


Reading of Papers at Meetings 


8, When submitting papers authors are requested to indicate whether they will be 
willing and able to read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


g. Postcards giving the programme of each meeting are issued some days before the 
ing concerned. Fellows wishing to receive such cards whether Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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MEETING OF 1952 OCTOBER 10 
Professor H. Dingle, President, in the Chair 


The President announced that Her Majesty Queen Elizabeth II had 
graciously consented to become the Patron of the Society. 

The election by the Council of the following Fellows was duly confirmed :— 

Ronald George Cox, 23 Ashleigh Gardens, Upminster, Essex (proposed by 
A. C. Stevenson); 

James Edmund Gibbs, B.Sc., A.R.C.S., D.I.C., 20 Norton Avenue, Surbiton, 
Surrey (proposed by R. W. B. Pearse); and 

Panos D. Solomonides, M.A., 84 St. Mary Abbot’s Court, London, W.14 
(proposed by Robert J. ‘Trumpler). 


One hundred and seventy-seven presents were announced as having been 

received since the last Meeting, including :- 

F. Becker, Einfiihrung in die Astronomie (presented by the author); 

Barbara Bell, A study of Doppler and damping effects in the solar atmosphere 
(presented by Harvard Observatory); 

L. Bragard, La géodésie et la méthode gravimétrique (presented by the author); 

Paul Coudere, The expansion of the universe, translated by J. B. Sidgwick 
(presented by the translator); 

M. N. Stoyko and Mme P. Dubois, La deuxiéme opération des longitudes 
Octobre-Novembre 1933 (presented by the International Union of Geodesy 
and Geophysics) ; 

Louise F. Jenkins, General catalogue of trigonometric stellar parallaxes 
(presented by Yale University Observatory); 

J. G. Porter, Comets and meteor streams (presented by the author); 

J. F. Scott, The scientific work of René Descartes (presented by the author); 

W. M. Smart, Astronomy, second edition (presented by the author); and 

Geschichte und Literature des Lichtwechsels der veranderlichen Sterne, zweite 
Ausgabe, Band III (presented by Potsdam Astrophysical Observatory). 


The President announced that Dr Edwin Hubble had found it necessary to 
postpone giving his George Darwin Lecture until 1953 May 8. 
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MEETING OF 1952 NOVEMBER 14 
Professor H. Dingle, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


V. Hugo Benioff, Ph.D., Seismological Laboratory, 220 North San Rafael 
Avenue, Pasadena 2, California, U.S.A. (proposed by B. Gutenberg); 
John Christopher Harrison, B.A., Department of Geodesy and Geophysics, 
Downing Place, Cambridge (proposed by R. I. B. Cooper); 

Neal John Heines, P.O. Box 2353, Paterson, New Jersey, U.S.A. (proposed 
by D. W. Dewhirst); 

Horace James Lewis, 96 Moring Road, Tooting Bec, London, S.W.17 
(proposed by F. M. Holborn); 

Charles Francis Richter, Ph.D., Seismological Laboratory, 220 North 
San Rafael Avenue, Pasadena 2, California, U.S.A. (proposed by 
B. Gutenberg); 

Karl Schiitte, Ph.D., Ortnitstrasse 15, Munich, Germany (proposed by 
F. Schmeidler); and 

William Reginald Slack, B.Sc., Avalon, Oxford Street, Sutton-in-Ashfield, 
Nottinghamshire (proposed by A. J. Ashmore). 


One hundred and twenty presents were announced as having been received 
since the last Meeting, including :— 


H. C. Urey, The planets—their origin and development (presented by 


A. S. L. I. B.); 

H. Mineur, Techniques de calcul numérique (presented by the author); 

lonospheric radio propagation, National Bureau of Standards Circular No. 462 
(presented by the National Bureau of Standards); and 

Charlotte E. Moore, Atomic energy levels, Vol. II, National Bureau of 
Standards Circular No. 467 (presented by the National Bureau of 
Standards). 


MEETING OF 1952 DECEMBER 12 
Professor H. Dingle, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


Hector Macdonald Adams, 2 The Normans, Bathampton, Bath, Somerset 
(proposed by K. Pollock); 

Sivagadadchahkurukal Subramanya Ayer, Madduvil, Chavakacheri, Ceylon 
(proposed by K. S. M. Sarma); 

Noné Chinery Bond, B.Sc., 33 Pereira Road, Harborne, Birmingham 
(proposed by M. C. Johnson); 

Kenneth Brierley, 279 Withington Road, Manchester 21 (proposed by 
W. Porthouse); 

Thomas Noél Burke-Gaffney, S.J., Riverview College Observatory, Riverview, 
Sydney, N.S.W., Australia (proposed by D. O’Connell); 
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William Buscombe, Commonwealth Observatory, Canberra, Australia 
(proposed by A. H. Joy); 

Raymond Coutrez, D.Sc., Royal Observatory of Belgium, 3 Avenue 
Circulaire, Uccle, Brussels, Belgium (proposed by R. H. Garstang); 
John Michael Anthony Danby, Christ Church, Oxford (proposed by 

G. L. Camm); 

Edward Stanley Farr, 64 Langton Court Road, St. Anne’s Park, Bristol, 
Somerset (proposed by J. G. Porter); 

Joseph Junkes, S.J., Ph.D., The Vatican Observatory, Castel Gandolfo, 
Rome, Italy (proposed by F. J. Hargreaves); 

Ronald George Mason, Geophysics Department, Imperial College of 
Science and Technology, London, S.W.7 (proposed by J. McG. 
Bruckshaw); 

Sidney Melmore, 61 Carr Lane, Acomb, York (proposed by H. Dingle); 

Jack Hobart Piddington, Ph.D., M.Sc., B.E., 19 Ocean Avenue, Edgecliff, 
Sydney, Australia (proposed by J. L. Pawsey); 

Colin David Reid, 27 Rake Way, Saughall, Chester (proposed by R. M. Baum); 

Allan Cecil Sanderson, 56 Borrowdale Road, Sefton Park, Liverpool (proposed 
by N. G. Matthew); and 

Paul Taylor, 26 Linksway, Upton, Chester (proposed by R. M. Baum). 


Sixty-eight presents were announced as having been received since the last 
Meeting, including :—- 
E. A. Beet, The sky and its mysteries (presented by the author); and 
A number of ancient books on astronomy (bequeathed by the late Dr Wilfred 
Hall). 


Mr K. C. Blackwell and Mr R. H. Garstang were appointed Honorary 
Auditors of the Treasurer’s accounts for 1952. 





ROTATION AND MAGNETISM IN THE WORLD-MODELS 
OF KINEMATIC RELATIVITY 


P. J. D. Gething 
(Received 1952 November 13)* 


Summary 


Milne’s theoretical relation between rotation and magnetism in the 
idealized world-models of Kinematic Relativity is examined. It is shown 
from considerations of symmetry that the rotation of a body located in a 
“fluid”’ of uniform density should produce a magnetic field of zero 
magnitude; if the mass distribution is asymmetrical with respect to the 
rotating body a non-zero field should be produced, although there appears to 
be no simple proportionality between the field intensity and the angular 
momentum of the body. These results do not appear to support Milne’s 
claim to have established a theoretical basis for a simple formula of the 
Blackett type. 





In a paper written shortly before his death, the late Professor E. A. Milne (1) 
attempted to establish a general relation between rotation and magnetism for the 
idealized world-models of Kinematic Relativity. He claimed that his equations 
could be reduced to Blackett’s empirical relation (2) if certain approximations 
were introduced. When applied to the rotation of the Galaxy, Milne’s relation 
predicted a field of the order of 10~" gauss near the Sun. 

Blackett’s relation now appears to be in poor agreement with observational 
data on the magnetic fields of certain stars. ‘There is, however, some evidence 
(3) for the possible existence of a galactic magnetic field. It is therefore of 
interest to re-examine Milne’s analysis in order to determine whether his equations, 
relating to idealized conditions, can be applied to an irregular distribution of mass 
such as the Galaxy. In particular, any assumption relating to the symmetry of 
the mass distribution must be investigated. 

A restriction on the type of mass distribution to be considered appears to 
have entered Milne’s analysis in equation (37) of his paper, 

dU U 

y der ke (1) 
where U is the angular momentum of the rotating body as measured by an 
observer employing a clock graduated to the ¢t-scale of time. In the derivation 
of this equation, Milne (4) assumed explicitly that any external force acting on 
the body must pass through the origin, i.e. the observer, and implicitly that the 
body was located in the substratum, a featureless fluid of uniform density. 
Thus the angular momentum should be constant according to classical mechanics ; 
equation (1) is not incompatible with the principle of the conservation of angular 
momentum, but it asserts that a quantity which is constant when measured in 
terms of the time scale of classical mechanics should vary in a certain way when 
measured in terms of an alternative scale. 

* Received in original form 1952 April 7. 
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It would appear that, in the presence of local irregularities in density, the 
simple relation expressed in equation (1) is no longer satisfied. Suppose, for 
example, that a massive particle M is introduced near a rotating body. The 
disturbing couple produced by the particle, calculated on the basis of classical 
mechanics, is known to be proportional to the mass of M but to be independent 
of U, the angular momentum of the rotating body. It is difficult to see how 
the relation dU/dt = U/t obtained by Milne could satisfy these conditions, even 
when the approximations which are permissible when distances and velocities 
are comparatively small are introduced. It must therefore be concluded that the 
term U/t does not include the effects of small irregularities in the symmetry of 
the mass distribution surrounding a rotating body on the angular momentum 
of that body. 

It seems probable tha: there should be two separate terms in the complete 
expression for dU/dt. ‘The first arises from the nature of the particular time scale 
employed and is equal to U/t. If an alternative time scale were used, the 
relation dU/dt=U/t would no longer apply, even in a mass distribution of 
perfect symmetry. In fact, when the 7-time scale of Kinematic Relativity is 
employed, we have dU/dr=0. 'The second term depends on the configuration 
of the local mass distribution, and may be expressed in the form 

dU 
dt 
where Fy, is the couple acting on the body due to a particle such as that at M. 

It would be reasonable to suppose that the two terms can be combined in 

the form 


=<ry, (2) 


dU U 
cme tal ee ie 
dt t + =F y. 


~ 


However, Milne’s treatment did not include such a relation and we shall not 
investigate its validity. In effect, he combined equations (1) and (2) into the 
relation 


= =<Pry, 


without, however, attempting to justify this procedure. It is suggested in the 
present paper that it is this step which is not legitimate, and which therefore 
vitiates Milne’s results.* 

If equation (1) is to be used, it appears necessary to postulate that the rotating 
body is located in a mass distribution of such a form that the total force on the 
body acts in the line joining the observer to the centre of mass of the body. The 
mathematical condition is that the sum of the external moments, LT .,4), should 
vanish. Let B denote the body, which may consist of a loose distribution of 
particles, located in a world-modelt W. ‘Then if O is the observer and A is 
the centre of mass of B, we must assume that any gravitational pull exerted by 
W on B acts along the line OA. This condition will clearly be satisfied if the 


* A referee has kindly pointed out to me that Milne clearly believed that equations (1) and (2) 
were alternative forms of the same relation. My interpretation of Milne’s equations must therefore 
be regarded as an expression of a personal opinion; I remain convinced, for the reasons set out 
above, that Milne’s point of view was untenable. 

¢ The distribution of mass surrounding the rotating body is referred to here and subsequently 
as the ‘ world-model ”’ since strictly it is necessary to sum the gravitational forces exerted by every 
particle in the universe. 
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distribution of mass in W possesses symmetry about the line OA, since there could 
then be no “preferred direction’’ in which the force could act except along OA. 
In this connection it is interesting to note that the velocity vector of B with respect 
to O, either alone or in combination with the vector OA, cannot introduce a 
preferred direction since the gravitational forces considered in Kinematic 
Relativity are independent of the velocities of the gravitating bodies. It will 
now be shown how these simple considerations of symmetry make it necessary 
to modify Milne’s analysis and conclusions. 

We begin by quoting four equations from Milne’s paper from which the 
relation between rotation and magnetism is derived. 

The force F, on a particular particle of mass m, belonging to the system B 
due to all the other particles in B and W is given by Milne’s equation (33), 

BS y (P, 0 P,) 

"= is” Bi |P,—P,|°’ 
where y is the gravitational constant and P,, P, are the position vectors of the 
particles m,, m, respectively relative to the observer; the sum is taken over all 
the particles in B and W except m,. The moment I, of F, about O is 


P,=P,aF,. (4) 


The supposed magnetic field H, relative to O arising from the action of a 
distribution of virtual electric charges (typical charge q, at P,) was obtained by 


Milne in the form 
8% ~. (0,—P,) 
eee area LL be 
where qg,=m,(y/2)'”. 
By combining equations (3), (4) and (5) above, Milne obtained 

P, =ctm,H,(y/2)"*. (6) 
It must be noted that equation (6), Milne’s equation (36), applies only to an 
individual particle m, of B. In order to find the total magnetic field of the 
system B it is necessary to sum equation (6) over all the particles of B. We 
obtain 


<m,H, =constant SI,. (7) 
B B 


In the most simple type of mass distribution, in which the system B contains 
only one particle m,, it would appear that the force F, must act along P,. Hence 
the moment of F, about O is zero, 


i.e. IT, =0, 
and from (7) H, =o. 


Thus it would appear that the magnetic field produced by the rotation of a 
system which can be regarded as a single particle is of zero magnitude. This 
result was not stated in Milne’s paper: however, it seems probable that he 
intended his analysis to apply to a more complex system, consisting of many 
individual particles bound together by mutual gravitational forces. For such 
a system, the forces acting on the typical particle m, need no longer lie along 
the vector P,; nevertheless, a result similar to that obtained for the single particle 
might be expected to apply to the more complex system, since the effects of the 
internal forces should cancel out when summed. 
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It is easily demonstrated that this supposition is correct. The gravitational 
force on each particle in B can be divided into two components, Fa) representing 
the internal force due to the other particles in B and F,,,,) representing the external 
forces exerted by the particles in W. Thus 

PL =P int) + Piextys 
where Py cinty = Pa A Ficint) 
and Pscext) = P, A F icext): 


The force exerted by m, on m, is equal and opposite to the force exerted by 
m, ON my, as may be seen from an examination of the form of equation (3); hence 
both EF at) and XP; ing) reduce to zero when summed over the system B. More- 
over, it has already been noted that, if the equation dU/dt = U't is to be used, 
we must confine our attention to world-models in which IF \,.,4)=0. Hence 
k 


LP, = Pay + [PF ceaty 
B b b 
=0, 
From equation (7) it then follows that 


<m,H, =0, 
B 


a result which appears to be at variance with the conclusions of Milne’s paper. 

The argument given above is so crucial in the critical examination of Milne’s 
theory that it is worth repeating in a slightly different form. For the sake of 
definiteness, we consider the rotation of an assembly of “ particles”’ constituting 
a galaxy, as viewed by an observer outside the galaxy. ‘The observer judges 
himself to be at the centre of a universe which exhibits spherical symmetry in 
the large-scale distribution of galaxies. ‘The arguments based on the property 
of symmetry can now be framed in a simple question: why should the observer 
find that the galaxy under consideration has a magnetic field acting in a particular 
direction ? 

It might be argued that the well-known motor and dynamo rules of electro- 
magnetism are examples in which two vectors of an orthogonal triad are sufficient 
to define the third member ; by analogy, the direction of the magnetic field possessed 
by the galaxy might be defined by the vector from the observer to the centre 
of mass of the galaxy together with the vector representing its motion (or possibly 
rotation). However, it has been shown above that this possibility is precluded 
by the form of Milne’s equations, which involve only one “ preferred direction”. 

It therefore appears to be necessary to accept the result that 


<m,H, =0 (8) 
B 


for the symmetrical models considered in the above analysis. ‘The physical 
interpretation which should be given to equation (8) is not clear. Milne attempted 
to compare the magnetic field represented by &m,H, with the strength of an 
equivalent dipole field, and by means of several simplifying assumptions he 
obtained an approximate relation between im,H, and pug, the magnetic moment 
of the equivalent dipole. Interpreted in this way, (8) becomes 


Per =0- 
If the world-model is asymmetric in the sense previously defined, the use of 
equation (1) in the above analysis does not appear to be justified. A non-zero 
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magnetic field could be produced in the rotating body whose magnitude is given 
by equation (7) as before, but from this point no progress can be made in relating 
the magnetic field of the body to the angular momentum. ‘There seems to be 
no reason to suppose that there is any simple relation between angular momentum 
and the strength of the equivalent dipole field produced. 

The results obtained above, for both the symmetrical and asymmetric world- 
models, are inconsistent with those obtained by Milne; it is suggested that 
Milne was not justified in the simultaneous use of the “ classical’? equation, 

dU — 
= 
and the equation of Kinematic Relativity, 
du U 
a te 
Milne’s claim to have established a theoretical basis for a simple formula of the 
Blackett type must therefore be disputed. 


—y 


I am very grateful to Dr G. J. Whitrow for many helpful discussions and 
suggestions during my course of study under his supervision. I am also grateful 
to the Department of Scientific and Industrial Research for providing a 
maintenance grant. 


Royal Observatory, 
Greenwich : 
1952 November 6. 
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ON THE THEORY OF WHITE DWARF STARS 
I. THE ENERGY Sources OF WHITE Dwarrs 


L. Mestel 
(Communicated by F. Hoyle) 
(Received 1952 May 9) 


Summary 

Present theories of the origin of white dwarfs are discussed; it is shown 
that all theories imply that there can be no effective energy sources present in a 
white dwarf at the time of its birth. The temperature distribution of a 
white dwarf is then discussed on the assumption that no energy liberation 
occurs within the star, and that it radiates at the expense of the thermal 
energy of the heavy particles present. In the resulting picture, a white dwarf 
consists of a degenerate core containing the bulk of the mass, surrounded by a 
thin, non-degenerate envelope. The energy flow in the core is due to the 
large conductivity of the degenerate electrons, while the high opacity of the 
outer layer keeps down the luminosity to a low level. Estimates of the ages of 
observed white dwarfs are given and interpreted. Finally, it is shown that 
white dwarfs may accrete energy sources and yet continue to cool off, provided 
the temperature at the time of accretion is not too high; this suggests a 
possible model for Sirius B. 





1. The origin of white dwarfs.—The problem of the structure of white 
dwarts was solved essentially by Fowler (1). He showed that under white dwart 


, 


conditions, the electron gas within these stars is ‘degenerate’, and exerts a 
huge pressure due to its zero-point “‘exclusion” energy. ‘The dominant term 
in the pressure formula depends strongly on the density, but is independent of 
the temperature, so the star is able to withstand its own gravitation without 
possessing a far higher temperature than the observed luminosities would 
warrant. The theory determines the radius of a homogeneous degenerate 
sphere in terms of its mass and the parameter ,, the number of electrons per 
proton mass. 

Stoner and Anderson, and later Chandrasekhar (2), modified the theory by 
using the relativistic Hamiltonian in the Fermi—Dirac partition function, 
and so deduced a more complicated pressure—density relation. The new 
mass-radius relation predicts smaller radii for stars of mass below 5:75 Mo/p,? 
(Chandrasekhar’s limit), while for masses above this limit no static spherically 
symmetric states exist. As all observed white dwarfs are below the limit, the 
existence of the limit is not serious for them, but the result has important 
evolutionary consequences. 

In order to understand how white dwarfs can come about, it is convenient 
to consider the behaviour of a mass of gas, obeying the usual pressure law, 
as it contracts under its own gravitation. The condition of hydrostatic support 
demands a central temperature varying roughly as the inverse of the radius (3), 
and, once the material has become opaque to radiation, the rate of contraction 
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will be determined by the rate at which energy flows down the temperature 
gradient from the centre to the surface. This is the Kelvin-Helmholtz 
mechanism of stellar energy generation. If the mass contains nuclear energy 
sources, however, the contraction is halted when the internal temperatures are 
sufficiently high for the energy liberated to balance the surface loss. (The 
possibility that rotational instability may set in before this stage is reached is 
neglected.) 

Within stars containing hydrogen and the merest trace of carbon and 
nitrogen, the Bethe cycle liberates large quantities of energy per second at 
temperatures above 10’ deg. K. There is also the direct synthesis of helium 
from hydrogen via the proton—proton reaction, which, if quantum-mechanically 
“‘allowed”’, is about as powerful as the Bethe cycle at 2 x 10’ deg. K. ‘The rate 
of energy liberation by both processes increases rapidly with temperature. 
It at once follows that a quantity of gas of the same mass as the observed white 
dwarfs cannot condense to a degenerate state while it contains energy sources; 
long before degeneracy is reached, the central temperature will be high enough 
for rapid energy liberation, and so the star takes up the normal equilibrium state. 
Thus any theory of the origin of white dwarfs must postulate the existence of 
stellar masses containing no energy sources throughout the bulk. As all 
observations of interstellar matter confirm that only the minutest fraction of it 
differs from hydrogen, the direct condensation of hydrogen-free material into 
a degenerate state can be ruled out. 

The first possibility that presents itself is that a white dwarf is the result of 
complete hydrogen-exhaustion in a main-sequence star of the same mass. Once 
all the energy sources have been exhausted, the star again contracts according to 
the Kelvin-Helmholtz theory. As the central density p, increases as 1, R* while 
the central temperature 7, only as 1/R, the criterion for degeneracy (4), 
10°p >T*?, must ultimately be fulfilled, and the star becomes a white dwarf. 
However, this line of evolution has against it one insuperable difficulty. Except 
for the white dwarfs, and a few sub-giants (discussed in the following paper), 
all stars of mass comparable with or less than the Sun’s are in the very early 
stages of their evolution (§), as is shown by the very small scatter about the 
lower end of the main sequence in the Hertzsprung—Russell diagram (6). It is 
quite unsatisfactory to imagine that the initial and final stages in the evolution 
of small stars are present in the heavens today, yet without any of the intermediate 
stages. We are driven to conclude that all the main-sequence dwarf stars are 
still in their youth, and that the white dwarfs, sub-giants and perhaps some 
sub-dwarfs have arisen through indirect lines of evolution. The upper limit 
thus found for the galactic age agrees well with estimates from independent 
evidence, such as the age of the Earth, the rate of disintegration of galactic 
clusters, etc. (5). 

For stars rather more massive than the Sun, this objection disappears 
because of the strong dependence of luminosity on mass: a star of mass 5Mo, 
say, can consume its hydrogen within about 10° years, well within the galactic 
lifetime. Such a mass, however, exceeds Chandrasekhar’s limit for hydrogen-free 
material and so cannot attain equilibrium even when degeneracy has set in. ‘The 
subsequent behaviour has been described in detail by Hoyle (7, 8). As the 
contraction proceeds, the conservation of the star’s angular momentum causes 
its angular velocity to increase, and sooner or later rotational instability sets in. 
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For stars with low angular momentum, this state will not be reached before 
densities and temperatures are so high that heavy elements are synthesized 
within the highly collapsed star. When enough material has been expelled by 
rotational instability for the remaining mass to be below Chandrasekhar’s limit, 
the star can settle down into a white dwarf state. The expulsion of matter can 
be either spasmodic or catastrophic, according to the density when instability 
sets in; these two cases are identified with novae and supernovae respectively. 
Observational support for the theory comes from the work of Baade and 
Minkowski on the Crab Nebula (9, 10). 

This provides a mechanism by which a hydrogen-free star of solar mass can 
be produced within the galactic lifetime. It is to be noted that the theory 
predicts the presence in some of these white dwarfs of considerable quantities of 
elements heavier than helium, owing to the synthesis of the heavier elements 
from the helium in the highly collapsed state reached by stars with low angular 
momentum, 

A rather different theory was put forward by Eddington in 1939 (1x). He was 
impressed by the fact that the observations of Sirius B seemed to demand a 
high hydrogen content for the star in order to account for its large radius. 
In order that the Bethe cycle should not start up during the contraction of a 
mass of gas, rich in hydrogen, it is necessary that carbon and nitrogen be less 
abundant than one part in 10'*; as there is no reason for other heavy elements 
to be present in large quantities if carbon and nitrogen are absent, Eddington 
concluded that the observations were very slightly in error, and that the true 
hydrogen content is 100 per cent. 

A further necessary hypothesis is that Fermi and ‘not Gamow~'Teller 
selection rules apply to the proton—proton reaction, which is therefore 
*‘forbidden”’; otherwise, this reaction will liberate enough energy per second 
when the central temperature is about 107 deg. K to bring the contraction to 
a halt. If this is the case, then “‘ Eddington” white dwarfs may be expected to 
condense, provided that for some period of the galactic lifetime the interstellar 
gas was effectively pure hydrogen. The relevance of Eddington’s theory 
depends therefore on the original composition of the galaxy. Hoyle (12) has 
discussed the early evolution of the galaxy on the alternative assumptions (a) that 
the galaxy originally consisted of pure hydrogen and that the Fermi selection 
rules apply to the proton-proton reaction, (6) the same with Gamow~'Teller 
rules, (c) that the primeval galactic matter was, as today, mainly hydrogen but 
with small quantities of other elements present. He shows that on each 
hypothesis the theory can account for the existence of white dwarfs, nova and 
supernova outbursts, and the subsequent distribution of heavy elements in 
space. The result of interest for the present work is that case (a) predicts the 
existence of ‘‘Eddington” white dwarfs, which should condense in large 
numbers early on in the galactic lifetime, before substantial quantities of heavy 
elements have been distributed in space by supernova processes and so 
adulterated the original hydrogen. 

It cannot be claimed that we have as yet sufficient evidence to decide between 
the alternatives. The problem of the constitution of the primeval galactic 
matter depends partly on the particular cosmology adopted; for example, the 
steady-state theory of the expanding universe, involving continuous creation, 
implies that galaxies condense from partially adulterated hydrogen, and so again 
rules out “ Eddington” white dwarfs. 
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The problem is a real one for cosmogony, in that the number of white dwarfs 
observed is very large, in spite of their faintness, and the indications are that 
their number is only less than that of the red dwarfs. But the proportion of 
massive stars in the solar neighbourhood is very much less, and the number of 
novae and supernovae per millennium observed in our own and neighbouring 
galaxies does not appear sufficient to account for the high proportion of white 
dwarfs. ‘Thus if ‘‘ Eddington” white dwarfs are ruled out for one reason or 
another, we are forced to assume that ail white dwarfs arise through the 
gravitational collapse of burnt-out massive stars. ‘This implies that earlier in 
the galactic lifetime there existed a much higher proportion of massive stars, 
at least in the solar neighbourhood. 

One further possibility might appear open; this is that a small mass of gas 
could contract to a degenerate state, without the internal temperature rising 
high enough for energy liberation to start. The mean temperature T of a 
quasi-static sphere of gas, of mass M, radius R, and mean molecular weight 
must (2) be greater than GuM/5AR, where G is the gravitational constant 
and A# the gas constant. In order that the Bethe cycle should not start up, this 
temperature must not rise beyond 107 deg. K. Degeneracy has set in when 
10%p~107%", where p is the mean density M/{rR*. ‘Together, these values 
give about 2 x 10% g for the upper limit to the mass, which is about one-tenth 
of the solar mass. ‘The time taken for such a mass to contract to a degenerate 
state is determined by the rate of outflow of energy from the centre. The opacity 
of a small dense star is given by Kramers’ law, giving for the luminosity 
4 x 10%( M/Mo)°*/(R/Ro)'*. (The conductivity of the non-degenerate electron 
gas is negligible (4).) The energy (thermal and gravitational) of a star is of 
order —GM?/R (2). The resulting time-scale for contraction to a degenerate 
state is therefore about 10 years. Hence, even if such a small mass could 
subsequently accrete enough matter to become a white dwarf of observable 
size, it seems unlikely that the galactic time-scale will allow this theory. 

2. The thermal properties of white dwarfs.—One fact stands out from the 
above discussion; all the suggested mechanisms for the production of white 
dwarfs demand that such stars are effectively without energy sources at the 
time when they become degenerate. ‘Thus, during their subsequent life one 
would expect them to radiate energy without replenishment. In spite of this, 
previous treatments of the temperature distribution within white dwarfs (13, 14) 
have assumed that the star is in a state of thermal equilibrium, with energy 
liberation balancing the surface loss. One reason for this is perhaps the high 
hydrogen content which the observations demanded for Sirius B. ‘This hydrogen, 
however, must have been accreted by the star after its birth as a white dwarf; 
this is discussed in detail in Section 4 and in the following paper. In any case, 
it is more logical to discuss the star as a cooling system, as it must begin its life 
as such, and then to investigate the effect of subsequent accretion. 

Many of the published treatments of the problem are misleading in that they 
seem to imply that there must be energy liberation within a white dwarf in order 
that it may shine at all. Thus, having built a heterogeneous white dwarf model 
in thermal equilibrium, Schatzman assumes that all white dwarfs must exist in 
this state—he refers in one paper (x§) to ‘‘|’état final hétérogéne existant, comme 
je l'ai montré, dans les naines blanches’’. But in fact neither a white dwarf 
nor a normal star needs nuclear energy sources to make it shine. The Kelvin- 
Helmholtz theory for a normal star is thermodynamically perfectly sound, and 
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predicts a luminosity fixed (at least to order of magnitude) by the mass, 
composition and central temperature of the star; it is only the time-scale difficulty 
which makes us look for nuclear energy sources that will keep the star in a state 
of thermal equilibrium. No similar time-scale difficulty arises for a white dwart, 
as it can supply its very low luminosity from its internal heat for a period 
comparable with 10 years. It will be shown that, as for normal stars, the 
luminosity of a white dwarf is fixed by its mass, composition and internal 
temperature. 

A normal star, in which the ordinary gas laws apply, and which possesses 
nuclear energy sources, always tends to approach a state in which the energy 
generation balances the surface loss. For by the virial theorem (2) 


3(y-—1)U+ Q=0, (1) 


where U is the total thermal energy and the total gravitational energy. The 
total energy E is U+Q, whence 


E= —(3y-4)U=(3y — 4)Q/3(y—1). (2) 


Hence, provided y>§, as is the case for highly ionized matter within stars, the 
thermal energy decreases as the total energy increases, and vice versa. Thus, if, 
say, the energy liberation exceeds the surface loss, the resulting increase in 
energy leads to an expansion and to a decrease in temperature throughout the 
body of the star. Since the energy liberation from nuclear sources is highly 
sensitive to temperature variation, after a small decrease in temperature thermal 
equilibrium will hold. Similar considerations apply if the star is radiating more 
energy than the sources supply—a slight contraction increases the temperature 
to secure thermal balance. A detailed treatment, taking account of the variation 
with temperature of energy outflow as well as of energy liberation has been given 
by Jeans (16). 

Thus normal stars have a “safety valve’? mechanism which prevents secular 
changes in the state of the star. White dwarfs have no such mechanism. The 
pressure of a degenerate gas depends almost entirely on the density, and is almost 
independent of the temperature. This means that a white dwarf whose energy 
sources are insufficient to supply its surface radiation cannot restore the balance 
by contracting ; contraction is prevented by the strong pressure, independent of 
the temperature, of the degenerate electrons. Whereas a normal star heats 
up as it loses energy, a white dwarf cools down, and its energy sources become 
still more inadequate to provide its radiation. Conversely, a white dwarf whose 
energy generation exceeds its surface radiation cannot restore the balance by 
expanding and cooling ; it heats up, continually generating more and more energy, 
until it ceases to exist as a white dwarf. Thus even if a white dwarf has, 
subsequent to its birth, acquired energy sources, there is no reason to equate its 
energy liberation to its radiation. At one internal temperature only does thermal 
equilibrium for a particular star hold, but at this temperature the star is 
‘balancing on a razor’s edge”; as soon as any departure from equilibrium 
occurs, there is no automatic adjustment, and the star departs further and 
further from this unique state. 

A permanent white dwarf must therefore be considered as supplying much 
of its energy-emission from its internal heat. It is convenient first to consider 
a model free of energy sources, and then discuss the effect of accretion. 
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3. The internal temperature of a white dwarf.—All save the outer fringe of a 
white dwarf consists of matter in a highly degenerate state (2). The degree of 
degeneracy is measured by the value of the parameter, A, defined (4) by 

ph* ° dz 
8ar.,mpy(2mkT)*? J, (1+e*/A)’ (3) 
where A is Planck’s constant, k Boltzmann’s constant, and m, mg the electron 
and proton masses respectively. ‘The effects of degeneracy begin to be felt 
when A~1, and p/T?#~10-§. When A~300, p/T?*~10~7; the pressure law 
Kp*® then holds very accurately. ‘Thus degeneracy may be said to be 
“complete” after A=300. Bars will be used to denote values at this point. 

Radiative transfer of energy is negligible compared with thermal conduction 
in the highly degenerate core, in spite of the reduced opacity of a degenerate gas. 
It will be shown that the thermal conduction is so great that, during the cooling 
of a white dwarf, the temperature is always very nearly uniform in the core. 

The thermal conductivity of degenerate matter is very nearly ap7, where 
# = 2°44 108 A/Z* c.g.s. units, and Z and A are the mean atomic number and weight 
respectively of the heavy particles (4). As the energy of a degenerate gas depends 
on the temperature to the second order only, the electronic contribution to the 
heat content is small, as in the theory of metals. Hence the specific heat is 
essentially that of the non-degenerate nuclei; it can be taken as $k/Amg, where 
A is the mean atomic weight of the po" The equation of cooling thus becomes 


I a oT 
Por («ores “> =) > ik yt at’ (4) 
where r is the distance from the centre, ¢ the time. 

The density distribution does not vary appreciably as the star cools; it is 
given in the non-relativistic approximation by the polytropic function for 
index n= #, or in the relativistic case by the Chandrasekhar function 
corresponding to the mass (2). ‘To simplify the problem, we shall use the 
polytropic function for p; it is easy to verify that the conclusions are substantially 
the same if Chandrasekhar functions are used. 

It is impossible to find a general solution for the non-linear equation (4). 
However, we are not interested in a general solution; we are interested in the 
form to which the solution approximates after a time sufficiently long for 
irregularities in the original temperature distribution to be smoothed out. 

Integrating equation (4) with respect to 7, 

3_kp OT : 
or 9 2 Amy rym (5) 





ap Tr? — 


Thus if 7,,, is the maximum value of |87/2t], 


oT - kT, 
te 2 2 
ap Tr? — > < aa in|, pr dr. 


Integrating again 


; 3 kT I . 
T,-T< gf ed [ pr dr bdr. 


Since 1/7 steadily increases with 7, this gives 


3 kT,, os ee dr 
“74% 2aAmyT (for dr) <. 
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The integral can be evaluated in terms of polytropic functions; substitute 


r=ez, p=p,u*?, (9) 
where p, is the central value and 


a i ee) 
= (2-5-0 ) : (10) 


Here K is the constant of proportionality in the non-relativistic pressure law Kp*" ; 
from statistical mechanics (4), its value is known to be 3:2 x (2/,)** x 10!” c.g.s. 
units, where ,, is as previously defined. The variable wu satisfies Emden’s equation 
for index n= 3/2, and can be found from British Association Tables. In terms 


of u and 2, 
T Tr dr 
(Ie) 2_ yz}, (11) 


, 1/2 _ 
T.-T< ¢ Bt I). (12) 


The order of magnitude of 7,, can be inferred from the observed 
luminosity L. The mass of the outer fringe being negligibly small, its 
contribution to the cooling may be disregarded. Thus L arises from the cooling 
of the degenerate part; its value is given by 

‘*3 kp OT 
L=— | 
02 Amg ot 
where 7 corresponds to the boundary of the degenerate matter. Thus if T, is 
a mean value of —d7/dt, 


Hence 


nr® dr, (13) 


3 kT 


c7 
L==-—— | npr? dr 
2 Amy “i , 
3 kT, 
we 2 Amy 
where M is the total mass. Using this equation in conjunction with (12), 
eon ee 
T,— T < —*~ (uw? -1). (15) 
aT MT, 
The maximum value of 7,— T comes at the boundary r =r of degeneracy, where 
we can put T=T7,u=u. Here 
Pe 2LT,, eee . 
T.-T < = (w"* - 1). (16) 
aT MT, 
As is shown in Section 4, the temperature at the onset of degeneracy in white 
dwarfs of observed luminosities is about 107 deg. K. The corresponding 
density p is given (4) by 10%p/u,73?=8-615, and so u can be found. Thus 
taking M=}Mo, L=L¢/400 as ‘ty pical white dwarf values, and assuming that 
T,,/T, is not much greater than unity, it is found that 


' ae a I I - 
T.-F<075 44 oe - cya} +10" (17) 


Hence, taking the mean atomic weight to be about 10, the relative temperature 
drop in the core is not more than about } and this drop decreases rapidly 
inwards. Hence the core is nearly uniform in temperature. (If 7/7) is much 


M, (14) 
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greater than unity, this simply means that near the centre of the star, where 
|@7/dt| is greatest in magnitude, the temperature is rapidly approaching equality 
with that farther out.) 
4. The rate of cooling of a white dwarf.—The thermal energy of the heavy 
particles in the bulk of a white dwarf is given by 
3 kT 
a Amy Gh 
where T is the temperature of the core. The rate of cooling is therefore given by 
3AM dT (19) 
aA a’ ie 
The dependence of L on T is found by integrating the equations to the 
non-degenerate outer regions. Before A= 1, conduction is found to be negligible, 
and the equations of support and energy flow are 
Rad GMp 
mm Ty dr?!) a Si 3 " 
T? dT Kop 1 L 
$ac— , =— T®3" (T) 4nr*’ 
Here f(7) is the guillotine factor, and »(7) the effective molecular weight, both 
functions of the degree of ionization of the material. Dividing (20) by (21), 
and integrating 


M, (15) 


L= 


32nacGM T 
eylLA Jo 
Both u(T) and f(T) are slowly varying functions of T (2), and so they may 

be taken from under the integral sign and replaced by their values ,,, f, at the 
strong maximum of 77°. Then (22) reduces to 
P__ ras 327acG My, f; gf 
i: 25'°5Kyp LA 
This solution holds until degeneracy comes into play, i.e. up toA~1. Then 
L=kT,??, (24) 
where & is constant, and 7; is the temperature at A=1. 
If this solution is continued until A= 300, the temperature increases only 
by a factor of about 2, and this is a gross overestimate, as conduction dominates 
after A~10. Hence the final equation to the cooling of a white dwarf is 


V 
a(3 arr) =47", (25) 


(pT)? = T?>u(T)f(T) aT. 


(23) 


~ dt\2 A 


where T is the internal temperature, given by (23). Numerically, taking the 
guillotine factor to be 10, as estimated by Chandrasekhar (2), T is given by 


y ad L oa 8 4 > 
t=(3) x 108 deg. K. (20) 


Hence the final model for a white dwarf is as follows: a degenerate, nearly 
isothermal core, containing the bulk of the mass, is surrounded by a non-degenerate 
radiative envelope with a high temperature gradient. The luminosity of the 
star, given by (24), is determined by the high opacity of the envelope, which 
acts as a blanket limiting the rate of cooling. 
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These results hold as long as most of the star is degenerate and the mass in 
the radiative envelope is small. It will be shown that this is the case over most 
of the star’s lifetime. The difference between the temperature at the onset of 
degeneracy (A=1) and that of the core will alter estimates of age by at most a 
factor 2, whereas only factors of 10 or more are of interest. 

Equation (25) integrates immediately to give: 

AM I 
ier . kT 
The time for the star to cool down from 7, to 7,, with corresponding 
luminosities L,, Ly, is therefore: 


=t+ const. (27) 


0-6 


A : eed «3 (28) 


Rh kl 

If 7, and 7, differ by a factor only slightly different from 1, the term 7,/Z, 
is negligible. From this, it is easy to estimate the time the star will continue 
to have a luminosity above a certain level, or emit light of a temperature above 
a given minimum. By inserting present values of JT, and Ly, the time that the 
star has been degenerate is found. 

For Sirius B and 40 Eridani B, the only two stars for which we have reliable 
data for both mass and luminosity, the results show that the times taken to 
reduce the present luminosity by a magnitude are about 10!°/2A years and 
101/34 years respectively, while Sirius B will have the same colour as the Sun 
10/4 years. By then, however, the luminosity will have 
dropped by a factor 20, and the star will be very difficult to observe. 

Large numbers of other white dwarfs have been observed, notably by 


RM (2 Ti) 


today in about 2°5 x 10 


Kuiper, but while their magnitudes and effective temperatures are fairly 


accurately known, the masses are seldom easily determinable. However, from 
Land T,, Rcan be estimated, and M then read off from Chandrasekhar’s graph (2), 
The values of M so found should be accurate enough for rough estimates of the 
ages of the stars to be made. 

Schatzman (17) gives a table of visual magnitudes; these must be corrected 
by means of Kuiper’s bolometric table (18). The following table gives the 
visual magnitude, internal temperature and age of a representative selection of 
white dwarfs. ‘The internal temperature depends on the composition of the 
star, in that the opacity of “* Russell mixture”’ is about six times as great as that 
of hydrogen; thus the internal temperature of an “ Eddington’”’ white dwarf 
should be slightly lower than one of the same mass consisting of helium and 
heavy elements. However, this error has only a small effect on the age of a star, 
which is influenced much more by the factor A in the last column. 


TABLE I 

Star Vis. Mag. Int. Temp. Age 
(deg. K) (years) 

Sirius B I1°4 X10’ <10°/A 
40 Eridani B II‘! 3 x10" x 10°/A 
Van Maanen 2 14°3 < 10° 10/4 
Wolf 1346 9'8 3 X10? 10°/A 
Wolf 457 15‘! 107 1014 
Wolf 489 15°1 3-5 x 10% 1011/4 
2 of Hyades 11°3 < 10’ 10°/A 
Ross 627 13°6 x 10" 10°°/A 
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The computed ages are really the times that the stars have been in white 
dwarf states; however, it can be shown that the time interval between the onset 
of degeneracy at the centre of the star and the appearance of degeneracy 
throughout most of the star is very short compared with the time taken by the 
star to cool. For the above analysis will hold as long as most of the star is 
degenerate, so that the mass in the radiative zone makes a small contribution to 
the cooling. If the mass of the degenerate core at any time be M,, then from 
the equation of support for the outer zone we have 


._ Gu (™ MdM_ Gp I (M?— M,?) (2 
eee f ae 9) 


The suffix 1 refers to conditions at the onset of degeneracy, and so 108p,~7;3?; 
also, the radius of the star must be less than the radius of the star at the time of 
onset of degeneracy at the centre of the star. Hence (2y) gives a lower limit to 
the temperature of the star when a certain fraction of its mass is degenerate. 
Assuming a star of solar mass to contract homologously, it is found that 
degeneracy sets in when 7,=2 x 10° deg. K, with the star’s radius then about 
3x10%cm. The available energy of the star in this state is ~24 MT,/A (a factor 
of 2 being introduced to take account of the gravitational energy that will be 
released during the further contraction). With a temperature at the edge of 


the core given by (29), the temperature gradient in the outer zone will be of order 


2/5 
= : fer ° a . (M?— M,?) x 108} (30) 


with a corresponding energy flow 
4ac Tj dT 


3% pi dr’ 


— 


47’. (31) 


When the mass of the core is about three-quarters the total mass, 7 (core) is 
about 108 deg. K, and the consequent flux at least 10 erg/s. Therefore the 
time for the star to become degenerate throughout most of its mass is at most 
2 x 107 years, and this is a huge exaggeration, as the luminosity is higher at higher 
temperatures. ‘Thus it is justifiable to assert that the star spends most of its 
lifetime in a degenerate state, and that the ages in Table I give the time since 
the star reached white dwarf dimensions. 

The tabulated lifetimes of the stars are satisfactory for the cooling theory. 
The results show that white dwarfs exist of all ages up to the age of the galaxy; 
this is to be expected, for the processes which give birth to white dwarfs—novae 
and supernovae—should occur at all times from 10° years after the birth of the 
galaxy (7). But, in addition to this, it is seen that some of the tabulated stars have 
calculated ages longer than the age of the galaxy if the material present does 
not contain substantial quantities of heavy elements. Helium is not sufficient 
to reduce the star’s calculated age satisfactorily. In this way the theory links up 
well with that of the synthesis of elements within collapsed massive stars. 

Finally, it is to be noted that stars in the neighbourhood of Ross 627 in the 
Hertzsprung—Russell diagram have ages of about 10'°/A years. This would 
fit in with the Eddington theory, which requires A = 1, and the age to be not very 
different from the galactic age (cf. Section 1). If the Eddington theory is not 
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accepted, then Ross 627 must be the core of a nova or supernova, which occurred 
about 10° years ago. 


5. The Hertzsprung—Russell diagram.—As a white dwarf cools, its radius will 
not alter by a sensible factor and so the surface temperatures decrease according to 


LeT¢. (32) 
Thus the path of a white dwarf in the Hertzsprung—Russell diagram is given by 
Mag. = —10 In T,,+¢, (33) 


the constant c depending on mass and composition. On this theory, therefore, 
there is no question of arranging the stars in a regular sequence depending on 
mass and composition, as with normal stars in thermal equilibrium. The 
position of a white dwarf of given mass and composition depends on its 
temperature, and hence on its age. The radius of a white dwarf which has not 
accreted any hydrogen is a function of the mass alone, and a star of given mass 
describes, as it cools off, an individual straight line in the Hertzsprung—Russell 
diagram. The greater the mass, the further to the left in the diagram is its cooling 
path. Thus the irregular spread of the white dwarfs over their part of the 
diagram is no reason for surprise. 

6. A white dwarf model with energy sources.—The next problem to be tackled 
is the effect of accretion on a white dwarf. A general discussion is given in the 
following paper; this section deals simply with accretion by a star cool enough 
for energy liberation in the accreted material to be less than the radiation from 
the star. 

The author was led originally to investigate this model by the well-known 
discrepancy in the radius of Sirius B. From its orbital motion relative to 
Sirius A, this star is known to have a mass nearly equal to the Sun’s. The 
radius was first determined by the effective-temperature method; later 
measurements of the Einstein shift were claimed to support not only the order 
of magnitude of the previously found M/R, but also its precise value. Errors 
of a factor 2 were declared very unlikely. A homogeneous star of this radius 
demands a hydrogen content of about 50 per cent if the number of electrons 
per unit mass is to be high enough to provide the necessary pressure. But if the 
opacity is due to heavy elements, the temperature gradient in the radiative envelope 
is high and leads to an internal temperature of about 2107 deg. K, with a 
consequent energy liberation far above the surface loss. 

A possible way out of the difficulty is suggested as follows. As discussed in 
Section 1, the hydrogen in Sirius B must have been accreted, and will therefore 
tend to lie on the surface. Hydrogen is less opaque than “ Russell mixture”’ 
by a factor 6, and so a lower temperature gradient is required to derive the 
observed luminosity. It will be shown that a composite model can be built, 
consisting of a hydrogen-free core (u,=2) surrounded by a hydrogen layer, 
such that the internal temperatures are low enough for energy liberation to be 
less than the surface loss. A necessary condition is that the proton—proton 
reaction be forbidden. 

The first task is to estimate the mass of the hydrogen layer in order that the 
observed and theoretical radii should agree. 

The pressure of a relativistically degenerate gas is given (2) by 


p=Af(x); p= Bx’, (34) 
43° 
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where 
A =6:01 x 10” c.g.s. units; B=9°82 x 105u, c.g.s. units 
and 
f(x) = x(2x? — 3)(x? + 1)"* + 3 sinh“ x. 
Chandrasekhar reduces the equation of hydrostatic support 
1d /(rdP 
a3(Gz) sa 


862)--@-3 


by making the following substitutions : 


y=? +1, 


(%)" I 
=i = “me Ns 
7G) °By,'” (39) 


y=Job> 
| a 2 
Yoo =X" + I. 





The suffix 0 refers to conditions at the centre. In this notation, the mass 
interior to 7 is given by 


A\3/2 PTT 
M(n)=47 (4) ‘ ie .(—7*¢ ). (40) 


Chandrasekhar has tabulated, for different values of yp, a set of solutions 
of (38), which start at the origin. A homogeneous star of given mass and 
composition has a density distribution given by one of these solutions continued 
outwards until its first zero, thus fixing the boundary of the star. In a 
non-homogeneous model, the distribution within the core will again be given 
by a particular one of Chandrasekhar’s solutions of (38) but, after the 
discontinuity in ,, further integrations must be performed. To simplify the 
problem, the gravitational effect of the hydrogen in the hydrostatic equation 
will be neglected ; this is justified by the final result that the mass of the hydrogen 
is well below that of the whole star. 

Then 


dp Se GpM 
y Sumeten (41) 


r? 


holds in the hydrogen zone; M is the mass of the core. Substituting from 
(34-36) this equation integrates to 


.2)1/2 ae s M [i I 
(1 +x") -1=2:66 x10". 57 .(7- 3), (42) 
where the condition of the vanishing of p at the boundary R has been used. 

In equilibrium the pressure, and hence the value of x, must be continuous 
across the interface between the two zones. As y.,=2 in the core and I in the 
layer, this condition leads to a drop by a factor 2 in the density and density 
gradient, thus distending the radius. 

Let 7 be the radial parameter at the interface. Then 7 is given by (39) 
and M by (40). To find M and y, a trial and error method is used. Assuming a 








No. 6, 1952 On the theory of white dwarf stars 595 


particular value of M, (—7?¢’) is found from the mass relation (40). If then 
a particular y, is assumed, Chandrasekhar’s tables supply 7 and ¢, whence x, p, 
and r are found with the help of (34) and (39). The correct yo for the assumed 
value of M is fixed by the relation (42) between r and x. As it is unlikely that 
the correct y, will be among Chandrasekhar’s tabulated set, it is necessary to 
interpolate between the tables to find_the appropriate yp. 

Thus to each assumed value of M, there_can be found a unique set of 
parameters. The final test, to fix the correct M, and so determine the problem, 
is the equality of (M— M) and 

1 
4r-R8 f. p2* dz. (43) 


r/R 
The function p is known from (34) and (42), and this integral can be computed 
in each case. 

An accurate set of results (within the limit of the approximation of neglecting 
the hydrogen mass in (41)) is 


M= 073M, r=0-30R, p(core) = 8-3 x 10° g/cm® (44) 


(using M=1-95 x 10% g, R=1-36 x 10° cm). 

Interstellar matter consists of hydrogen with small quantities of other 
elements intermingled; Dunham (19) gives a proportion of one part in 10‘, and 
certainly I per cent is an overestimate. Such small quantities will not affect 
either the opacity or the conductivity, but any carbon and nitrogen present will 
start up the Bethe cycle at sufficiently high temperatures. It is therefore 
important to estimate an upper limit to the temperature at the interface between 
the core and the hydrogen layer. 

The temperature at the point A =1 is given by (28) of Section 3; the factor Ky 
is one-sixth of that for ‘‘ Russell mixture”’, while the guillotine factor may be 
put equal to I, as pressure ionization sets in early. The value found is 
T=2-91 x 10° deg. K. 

Between A =1 and A= 300, account must be taken of incipient degeneracy and 
thermal conduction. The equations are 
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where v is the thermal conductivity, K the opacity and p, and p, the nuclear 
and electron pressures respectively. The equations must be integrated 
numerically, using previously tabulated values of p,/p, and v (4), and Marshak’s 
expression for the opacity of incipiently degenerate matter (13). Between 
A=1 and A=10, conduction dominates throughout the star. ‘The temperature 
at A= 300 is found to be at most 6-9 x 10° deg. K. 

An upper limit to the temperature at the interface is found by neglecting the 
cooling of the hydrogen and integrating through the degenerate hydrogen the 
equation 


SE ae 
5( apT =) =0, (47) 
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using for p the function given by (34) and (42). The resulting maximum 
temperature is 8-4 x 10° deg. K. 

Energy will be liberated by both the Bethe cycie and the proton—proton 
reaction. On Bethe and Critchfield’s theory (20), the proton—proton reaction, 
if allowed, liberates 2 erg g™1 s™1 at the centre of the Sun, the conditions assumed 
being p=80 g/cm*, 7 =2x 107 deg. K, and hydrogen 35 per cent by mass. 
The energy liberation varies as p7*. From this, a minimum for the energy 
liberation by the proton—proton reaction in Sirius B can be calculated, and the 
result exceeds the surface loss by a factor of about 100. Hence this model 
demands that Fermi and not Gamow-Teller selection rules apply to the 
proton—proton reaction, in which case the energy liberation is reduced by at 
least 10. The evidence from nuclear theory is indecisive as to what interaction 
to assume, and the reaction is far too slow to be measured at terrestrial 
temperatures. Bondi (21%) suggests that the shape of the lower end of the main 
sequence in the L— 7, diagram is best explained by assuming the reaction to be 
allowed, but before this can be decided with certainty, further detailed 
integrations of stellar models of small mass are needed, taking into account 
incipient degeneracy and thermal conduction. 

The Bethe cycle generates 100 erg g™! s™4 at the Sun’s centre, assuming the 
previous data, and a 10 per cent nitrogen concentration (20). The reaction rate 
varies roughly as p7™*. ‘Taking I per cent to be the proportion of nitrogen 
present in the accreted material, the maximum amount of energy produced per 
second in the zone is only zg of the luminosity. Hence on the previously stated 
assumption, the model satisfies the theoretical mass—radius relation and is stable. 

After this work was done a report appeared (20, p. 293) that the generally 
accepted radius for Sirius B is wrong, and that the correct value of the Einstein 
shift is such as to give a radius requiring almost zero hydrogen content. (The 
hydrogen lines visible on the surfaces of Sirius B and 40 Eridani B require only 
a surface layer of negligible mass to account for them.) This would mean, of 
course, that the effective temperature of the star is wrong by a factor 4/2, in spite 
of what the spectroscopists have been reporting during the last thirty years. 
Other workers are of the opinion that the probable errors in such delicate 
spectroscopic measurements are too great for reliance on the values of the radii 
determined from the Einstein shift. The question, therefore, is still open, and 
it seems worth while publishing the details of the above model in order to 
illustrate what happens to a cool white dwarf which accretes hydrogen, even if 
the model should turn out to have no relevance to Sirius B. 

The models considered by Schatzman are similar to this one in structure, 
but he assumes that the temperature is always at the right level for the energy 
liberation in the hydrogen to balance the surface loss. The objections to this 
assumption have already been stated. That Schatzman’s calculated luminosities 
agree with the observed magnitudes is no reason for surprise; for as the energy 
liberation varies as 7™*, while the surface loss only as 77, a slightly higher 
temperature would be sufficient for the energy liberation to balance the 
luminosity, without the luminosity being increased very much. 
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ON THE THEORY OF WHITE DWARF STARS 
II]. THe ACCRETION OF INTERSTELLAR MATTER BY WHITE DWarFs 
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(Communicated by F. Hoyle) 
(Received 1952 May 9) 


Summary 


The behaviour of a white dwarf subsequent to its accreting substantial 
quantities of interstellar matter is shown to depend critically on the star’s 
internal temperature at the time of accretion. A cool star merely develops 
a non-homogeneous structure and continues to cool off, as in the model for 
Sirius B suggested in the preceding paper. A hot star changes its form 
radically; depending on circumstances, it may evolve into a sub-dwarf or 
possibly a sub-giant state, or it may explode like a gigantic hydrogen bomb 
and be observed as a type of supernova. 





1. Introduction.—In the preceding paper it was shown that a white dwarf 
which has acquired energy sources cannot in general be in thermal equilibrium. 
The degenerate part of the star is at an approximately uniform temperature, and 
the generation of energy by the star balances the surface loss only if this uniform 
temperature has a particular value. If the temperature differs from this value, 
the white dwarf has not, like a normal star, the power of adjusting itself to secure 
thermal equilibrium; a star once cooling must continue to cool and a star once 
heating up must continue to heat up. It follows that a white dwarf which has 
cooled to a central temperature below about 107 deg. K is not sensibly affected 
by accretion of hydrogen, as the energy liberation within it is only a fraction of 
the luminosity. On the other hand, if the star’s temperature is high enough, 
accretion may lead to an energy liberation greater than the surface loss. In this 
case the star’s thermal energy increases, and the discrepancy between energy 
liberation and surface loss is increased still further. The star therefore retraces 
its path of cooling. 

In this paper two main cases are discussed in detail: accretion by a white 
dwarf of temperature slightly above 107 deg. K and by a hot collapsed star which 
has hardly begun to cool off into a white dwarf state. These two extreme cases 
illustrate the main alternatives. 

2. The distribution of accreted matter through the star.—Accreted matter will 
initially form a layer on the surface. But a white dwarf which has arisen from 
the rotational instability of a collapsed burnt-out star is likely to be rotating fast, 
and the well-known von Zeipel-Eddington argument applies as much to white 
dwarfs as to normal stars (1). The rotation tends to lead to inequalities of 
temperature over a level surface; meridional currents result, which tend to 
homogenize the star’s material. It is therefore important to estimate the order 
of magnitude of the current velocities. 

Due to the rotation, there is a mean excess or defect of energy of order 
LQ*r/Mg erg g-4s~* over a level surface, where Q is the angular velocity, and 
g gravity (2). This must be balanced by the energy transported by the currents. 
Degenerate matter at the absolute zero is isentropic, and so the degenerate 
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electrons make virtually no contribution to the energy transport. The main 
contribution comes from the non-degenerate heavy particles. ‘The value finally 
found for the current velocities is (2) of order 
L Qr (n+1) (1) 
M* g* *(n—3)’ 
where n is the polytropic index at the point considered. 

If a white dwarf were a polytrope of index n=#% as predicted by the 
non-relativistic theory, the formula (1) would be invalid, as the heavy particles 
would also form an isentropic gas and only the thermal part of the electronic 
energy would contribute to the convection of energy. However, the relativistic 
theory gives a polytropic index between $ and 3, and so (m+1)/(n—3)~1 for 
an order of magnitude estimate. Then on putting in a typical white dwarf 
value for L/M, the upper limit to the current velocities is found to be ~10~" cm/s; 
this is obtained by assuming that centrifugal and gravitational forces are equal, 
implying that the star is on the verge of disruption. Hence the time for matter 
to traverse the star’s radius is at least 10! years, which is far too long for mixing 
to be significant. 

Schatzman’s theory of white dwarfs (3) is based on the sorting of the elements 
in the strong gravitational field of the star. He shows that in statistical 
equilibrium, the hydrogen is completely separated from the other elements at 
the absolute zero, while at finite temperatures there is a zone of mixing, of depth 
about one hundredth of the star’s radius, in which is liberated the energy 
demanded by Schatzman to balance the luminosity. He does not discuss how 
the star evolved into this state nor the stability of the stratification against 
rotation. In the present approach, the hydrogen is initially in this state of 
separation as it has been accreted by a hydrogen-free star. If rotational mixing 
had proved to be important, it would have been necessary to discuss the reverse 
process of the diffusion of the elements to the equilibrium state of stratification 
to decide whether rotational mixing or diffusion dominated. As, however, 
rotational mixing is negligible, and the hydrogen is initially separated, no 
discussion of diffusion is necessary. 

3. A theory of supernovae.—The energy liberated per second in a surface 
layer of hydrogen accreted by a white dwarf is given (4) by 

Er*e, (2) 
where 7 =(150)/(7/10°)**; (3) 
T is the temperature of the degenerate core of the star. The factor E depends 
partly on the proportions of carbon and nitrogen in the hydrogen but mainly 
on the mass of hydrogen which is in the degenerate region and so at temperature 7’. 
However, the energy liberation is much more sensitive to temperature than to 
the other parameters; at 10’ deg. K (2) varies roughly as 7"*. In the model 
discussed in the previous paper, the luminosity exceeds the energy generation 
by about 1000 when the temperature of the core is about 0°84 x 107 deg. K. 
At twice this temperature, the situation is reversed, and energy generation 
outweighs surface loss by a large factor. ‘Thus energy liberation and surface 
loss balance at a little above 10’ deg. K, and this result is not much affected by 
the mass of the accretion because of the strong temperature dependence. 
In what follows, therefore, E is fixed by assuming thermal equilibrium at 


10’ deg. K. 
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If 7 is the temperature of the degenerate core, the density of the material 
at the point of onset of degeneracy is given by 10899~7%*, This density increases 
as the star gets hotter, and the sphere of onset of degeneracy moves inward. 
At 107 deg. K, p=300 g/cm*; at 108 deg. K, p=104 g/cm*. Now the future 
behaviour of the star is strongly dependent on the temperature at which all the 
hydrogen ceases to be degenerate. Chandrasekhar’s table (§) for the density 
distribution within a homogeneous white dwarf of solar mass shows that the 
sphere of density 10* g/cm* contains about 99°5 per cent of the mass, while that 
of density 10° contains 92 per cent of the mass. Thus only a relatively small 
accretion is sufficient for part of the hydrogen to be degenerate even if the core 
temperature exceeds 10° deg. K. The effect of much smaller accretions is 
discussed in the next section. 

Now suppose that the star begins to accrete substantially when its core 
temperature is about 10’ deg. K. The luminosity of the star is given by kT7?, 
while the energy liberation is proportional to 71%. Hence a very small excess of 
temperature above the critical value 10’ deg. K is sufficient for energy liberation 
to outweigh the surface loss, and the star begins to heat up. The equation to the 
heating of the star is then 

AM d 


eed. 4 ae 18 __ 7/2 
y jake iad cilities aioe (4) 


Writing LZ, for the luminosity at 7’=10’ deg. K, and fixing B by equating the 
energy liberation at this temperature to 1), equation (4) becomes 


3 2M aT _ T \38 T \72 
2 4 & ~ (5) - (3) (5) 


Once the temperature is slightly above the critical value 10’ deg. K, the luminosity 
term on the right is negligible. ‘Then (5) integrates to give 


3 2M s _ T's 


1-3 2M) 78 T,\8 6 
2 7A ) Ly (3) i. (3) (0) 


for the time of heating from temperature 7, to 7,. Assuming M= Mo, A~10, 
and Ly, the luminosity of the star in the unique state of thermal equilibrium, is 
L/100, equation (6) gives a time-scale of about 107 years for heating from 
T=10' to 10°deg.K. Of this time, all but about 100 years is spent below 
2 x 10’ deg. K, and all but about an hour below 4 x 10’deg. K. Thus the heating 
proceeds very slowly at first but rapidly accelerates. 

For several reasons equation (5) will not remain accurate as the star heats up. 
First, the stock of hydrogen may give out. The energy liberated per gram of 
hydrogen when transmuted to helium is c?/125erg. Hence the minimum 
accretion necessary for the star to reach 10° deg. K is (@M x 108/A)/(c?/125) g, 
or about }x 10g (again taking A~10), which is only one ten-thousandth 
of the star’s mass. Thus, even a small accretion is sufficient to provide enough 
energy to heat the star to 10° deg. K. 

Secondly, it has been assumed that the mass of hydrogen at the core 
temperature remains constant over the period of heating. In fact, the sphere 
of onset of degeneracy moves inward, so that less and less of the hydrogen 
remains at the maximum temperature. 
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Over most of the period of heating of the star this is not important, as the 
temperature remains near 10’ deg. K, but, in order that (5) should still apply 
at 10° deg. K, it is necessary that enough hydrogen be accreted for a substantial 
quantity of it to remain degenerate at this temperature. An accreted mass of 
one-tenth the stellar mass is more than ample. 

The 77§ law remains valid only near 107 deg. K, but again this restriction on 
equation (5) is unimportant for time-scale calculations, as the star spends most 
of the heating period at low temperatures. It is easy to verify that the times 
taken to heat up from 2 x 107 deg. K and from 4 x 107 deg. K are still very small 
fractions of the total time, even when the lower indices demanded by (2) and (3) 
at higher temperatures are used. 

Finally, it has been assumed that the energy liberated near the edge of the 
star is distributed over the bulk by conduction. The rate of conduction of 
energy from the surface is of order 


.| aT , 
ap 7 Fe 4m R?. (7) 


In order to be able to conduct into the body of the star the energy being liberated 

at the edge, which is L,(7/10’)'*, the temperature gradient inwards must be 

of order 10-*(7/10’)'*._ This shows that the core remains nearly isothermal 

as long as (7/10’)!*<10, which by (6) is for nine-tenths of its heating period; 

after this, the energy liberated at the edge cannot be conducted inwards in time, 

and so the hydrogen heats up even more rapidly than previous estimates suggest 
if the hydrogen is one-tenth of the star’s mass, ten times as fast. 

The time-scale of 107 years is important in that it gives the star enough time 
to accrete a substantial quantity of hydrogen. This time can be increased if 
initial temperature is critical, so that the energy liberation hardly exceeds the 
surface loss and the right-hand side of (5) is a small quantity. In this case the 
heating time can be increased to 10° or 10° years, giving the star even better 
opportunities for accretion. 

Thus a white dwarf of temperature a little above the critical level, which 
accretes a substantial quantity of hydrogen will heat up over a period of at least 
10’ years, first very gently but accelerating later. ‘This is simply a consequence 
of the thermal properties of degenerate matter ; as long as the star is a white 
dwarf, no expansion against gravity can occur to absorb heat energy and cut 
down the rate of energy liberation, as would occur in a normal star. 

For definiteness, the mass of the white dwarf will be assumed to be the same 
as the Sun’s, and it is supposed to accrete hydrogen of about one-tenth its own 
mass. ‘Then degeneracy persists within much of the hydrogen up to temperatures 
in excess of ro*deg. K. But at this temperature the Bethe cycle proceeds with 
phenomenal rapidity. The energy available in the hydrogen layer is 
(Mo/10) . (c?/125)~1-4 x 10% erg. At ro®deg.K this is liberated at a rate 
given by (2) and (3), of order 10“ erg/s ; at 2x 10%deg. K at about 10 erg/s. 
It has been seen that this energy is not distributed over the star’s body but 
remains near the edge, heating up the hydrogen still further: the layer is heated 
from 10° deg. K to 2 x 108 deg. K in about 10 seconds. 

Now as the temperature of the hydrogen increases, it gradually ceases to be 
degenerate. If the temperature were rising slowly, the star could gently evolve, 
with the outer layers expanding and cooling, until thermal equilibrium is reached. 
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But by the time the white dwarf has reached 10° deg. K, it is liberating energy 
much too rapidly. If some of the hydrogen still remains degenerate at this 
temperature—as is the case for an accretion of mass M,/10—then this hydrogen 
is heated up still more. At 2 x 10*deg. K, degeneracy is disappearing from the 
hydrogen, and it can expand and absorb heat energy. The gravitational potential 
energy of a mass M_/10 in the field of the star is of order -GM.?/10R, and as 
R increases this increases at a rate 

Set (8) 

10R? 

In order to absorb just the energy being liberated at 2 x 10° deg. K (10° erg/s), 
the layer would have to expand at about 10! cm/s (even neglecting the increased 
inertia of matter at velocities comparable with c) and this rate of absorption 
merely keeps the temperature, and hence the energy liberation, at a constant 
level, without cutting them down. But the velocity of sound at these temperatures 
is of order 2 x 10°cm/s, and the speed of ejection of the matter can be at most 
of this order as long as pressure gradients are finite. Energy liberation therefore 
continues at a phenomenal rate, and most of the energy available in the hydrogen 
is liberated within a few seconds. Such a quantity is very much greater than the 
negative gravitational potential energy of the outer layer in the star’s field. ‘Thus 
it is concluded that shock conditions must set in, and the outer layers of the star 
are blown off to infinity. 

The velocity of the moving matter must exceed the velocity of escape 
1/(2 GM/R) from the star, which is about 5-2 x 10° cm/s for a star of solar mass. 
Then the energy of the explosion is given by the kinetic energy of the moving 
matter. If a mass of one-tenth the solar mass is blown off, this gives an energy 
of about 10“ erg. 

The nearest terrestrial analogue to this phenomenon is a hydrogen bomb. 
Just as in an atomic weapon, rapid nuclear reactions occur under conditions which 
prevent gradual expansion to meet the rising energy production. ‘The behaviour 
of the star is strictly determined by its composition and mass, and it cannot 
prepare at t= 0 for what might happen to it at t= 107 years and adjust its behaviour 
toavoidanexplosion. By the time degeneracy has disappeared from the hydrogen 
zone and gravitation begins to absorb some of the thermal energy, the 
catastrophe is upon it. 

It is tempting to identify this phenomenon with a type of supernova. ‘There 
is evidence (6) that there are at least two types of supernova, and many pre- 
supernovae appear to belong to the white dwarf region of the Hertzsprung— 
Russell diagram (7). Also, the velocities and kinetic energies of observed 
supernova remnants are of the correct order of magnitude (7). 

Type I supernovae have greater mean absolute magnitudes at maximum 
than ‘Type II, and also differ in light curves and spectra. Type I is found in all 
nebular types, while Type II is limited to systems containing massive O and B 
stars (Baade’s Population I systems (6)). Provisionally, therefore, Type I1 
supernovae are identified with the gravitational collapse of a massive star, as 
discussed by Hoyle, and Type I with the thermal blow-up of a white dwarf. 
This implies that white dwarfs are present in Baade’s Population II systems as 
well as his Population I systems. White dwarfs are difficult to observe in globular 
clusters, which we know to be pure Population II systems, but the additional 
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criterion of stellar velocities shows that some white dwarfs in the Sun’s neigh- 
bourhood have large proper motions and so are included in Baade’s Population 
II systems (6). 

The above theory differs from that put forward by Schatzman (7), both in 
the details of the heating of the star and in the cause of the explosion. He assumes 
that thermal equilibrium holds throughout the lifetime of a white dwarf ; he 
then shows that certain of his heterogeneous models will heat up as the proportion 
of hydrogen decreases, as higher temperatures are required in order that energy 
generation should balance surface loss. At a critical temperature convective 
instability sets in, and Schatzman suggests that a secondary effect of this 
turbulence is to set up a shock wave to the surface of the star, thus giving rise 
to a superficial explosion. The objections against assuming thermal equilibrium 
have been put forward in the previous paper. Convective instability sets in 
when 

1dT /1 dP 

T dr | Pdr (9) 
exceeds 0:4. It is found that this does not occur until temperatures are of order 
108 deg. K, when thermal catastrophe is already upon the star. 

There is one further possibility—a fairly massive white dwarf near 
Chandrasekhar’s limit might accrete sufficient matter to make it approach the 
limit of mechanical stability. This happens when yo, the parameter defining 
the density distribution in the core (§), approaches infinity. In this case, the 
star begins to contract with the temperature rising slowly, as not all the 
gravitational energy released will be absorbed by degenerate electrons. ‘Thus 
the star either becomes rotationally unstable or heats up and explodes. 

4. Non-spectacular evolution.—In the last section it was assumed that enough 
hydrogen is accreted for some hydrogen to remain within the sphere of onset 
of degeneracy until a high temperature has been reached. A very small 
accretion would cease to be degenerate at low temperatures, when catastrophic 
conditions have not been reached. Then the hydrogen will continue to pump 
in energy to the star, but its temperature soon drops as the normal gas law 
applies at the extreme edge of the star. The final equilibrium state reached 
consists then of an isothermal, hydrogen-free core, most of which is degenerate, 
surrounded by a non-degenerate hydrogen envelope generating enough energy 
at the shell surrounding the core to balance the luminosity. 

Such a star can arise in other ways as well. Suppose, now, that a white 
dwarf begins to accrete while it is still hot, perhaps while only a smali part of 
it is degenerate, and with internal temperatures of order 10° deg. K. Hydrogen 
accreted by such a star would be rapidly transmuted to helium and its available 
energy liberated. If the rate of accretion is greater than about 10! grams per 
second, the energy liberated will exceed even the high luminosity of the star, 
and energy will be pumped into it. The important difference from the previous 
case is that the star has no time to acquire a large stock of hydrogen; accretion 
is a continuous process, and at 10° deg. K energy is liberated so rapidly that the 
star changes its state before a stock of explosive hydrogen has been built up. 

In a white dwarf, the thermal energy U and the total energy # increase and 
decrease together, whereas in a normal star 5E and 6U are of opposite sign. 
As energy is pumped into a hot white dwarf, sooner or later U must change sign 
and the star begins to approach a state of thermal equilibrium. When the 
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temperature of the star has dropped to well below 10° deg. K, accreted hydrogen 
will not be transmuted quickly, as the Bethe cycle is now much slower, and so 
the star can acquire a stock of hydrogen. When the hottest hydrogen is at about 
2 x 10’ deg. K, the energy liberation is of the same order as the luminosity, and 
a state of thermal equilibrium is reached. 

Again the question of rotational mixing must be examined. Formula (1) 
again gives the rate of circulation, but here ‘‘ L”’ stands for the rate E at which 
the energy of the star is increasing due to the excess of energy liberation over 
surface loss. ‘Thus the time for the currents to travel a distance of the order 
of the star’s radius is 


2 M 
as : (10) 


while the time which the star takes to change its energy, and so E, by a factor 
of about 2 is 

GM? 

RE © 
But (10) is always considerably greater than (11) as Q?r<g for stability. Thus 
during the time in which the star is evolving towards thermal equilibrium, 
mixing is insignificant, and so the final state the star reaches is again a shell-source 
model. It is well known (8) that no static star can have a non-degenerate, 
isothermal core of mass greater than about one-tenth of the total mass, but if 
degeneracy exists in the core then the star can achieve equilibrium (g¢ The 
final state will consist again of a partially degenerate hydrogen-free core, 
surrounded by a hydrogen envelope. From Miss Harrison’s calculations (9), 
it is seen that stars of this structure all lie below the main sequence, being all 
of small luminosity and nearly all of small radii. The brightest should fall 
in the sub-dwarf region, and the faintest among the white dwarfs. 

The speed of rotational currents again becomes of importance after the star 
has reached its equilibrium state, for the time-scale of hydrogen consumption 
is very long compared with 10’ years. But the time for material to traverse 
the radius of a star of constant angular momentum is inversely proportional 
to the luminosity (2) and, even if the star had in its original white dwarf state 
the maximum angular velocity allowed for stability, it would still take something 
like 10" years to mix its hydrogen thoroughly. Only if the accreted matter 
possesses a higher angular momentum per unit mass than the star can the time 
of mixing be brought down to within the galactic lifetime. If this is so, then 
some of these sub-dwarfs will become homogeneous and would reach a state 
with a convective core and a high » throughout. A subsequent slight accretion 
would then have the effect of distending the radius, according to the Hoyle- 
Lyttleton giant theory (10). 

This rather remote possibility is worthy of mention because there appears 
to exist a small group of stars which have luminosities of solar order but are to 
the right of the main sequence in the Hertzsprung—Russell diagram. Struve (11) 
calls them “sub-giants’’ and gives as examples Delta Eridani and Mu Herculis, 
of spectral type Ko and G7 and absolute magnitude 3-8 and 3°7 respectively. 
The present giant theory (10, 12) can account for their structure by postulating 
a discontinuity in » ; but such stars cannot have arisen by transmutation of 
hydrogen within a dwarf star, again because of the long time-scale demanded. 


(11) 
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If the attempt is made to build a consistent cosmogony, one is forced to look for 
a special line of evolution for these stars. Explanations based on assuming 
that the original material of the stars had a high content of helium and other 
elements face the difficulty of the non-existence (as far as is known) of gas clouds 
containing anything but the smallest traces of elements other than hydrogen. 
So if there is no mixing process (rotational or other) strong enough to allow a 
white dwarf to become a sub-giant, an alternative special explanation of these 
stars must be found. 
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SURFACE PHOTOMETRY OF SOUTHERN ELLIPTICAL NEBULAE 
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Summary 


After the investigation described in Paper I there remained seven nebulae 
south of declination —26° and brighter than integrated magnitude 12:0, 
which have been classified as elliptical. ‘These objects have now been investi- 
gated photometrically. The true elliptical n-bulae are four in number : 
three objects previously so classified are spirals. ‘The four true ellipticals 
conform closely to the type of intensity distribution proposed by Hubble. 
The nucleus of the spiral nebula, NGC 6744, also conforms to this rule. 
Some possible consequences of this fact for our ideas on nebular evolution 
are discussed. 





In a previous paper (1) the surface photometry of a number of elliptical 
nebulae south of declination — 26°, and having an integrated magnitude brighter 
than 12-0 according to the Shapley-Ames catalogue, was described. ‘here 
remained seven nebulae satisfying these criteria, and these have now been 
investigated. In the event, only four of these were found to be true elliptical 
nebulae. Table I sets out the data for these objects, and for the plates (representing 
a selection from a much larger number obtained at the Newtonian focus of the 
74-inch Radcliffe reflector) used in the study. Photography was only attempted 


TABLE I 
Object R.A.(1950) Dec.(1950) Mag. Plate No. Emul xp. Date 
h m ew Y 1950 
NGC 4976 §=13.05'9 —49 14 11°6 A 990 I June 9 
A 992 I June 10 
NGC 5102 13 19°1 — 36 23 10'8 A 975 1.A. ‘ June 2 
A 743 I1A.Z Apr. 17 
IC 42096 13 33°8 —33 43 11°9 A 993 1 A.Z June 10 
A 994 I ‘ June 10 
IC 1459 22 54°5 —36 41 11°3 A 1271 103a-O ‘ Oct. 3 
A 1275 1o3z3a—-O c Oct. 4 
I.A.Z. =Ilford Astronomical Zenith. 
103a—O = Kodak 103a-O. 


on moonless nights with conditions as nearly perfect as possible, and with the 
objects on the meridian. ‘The methods of calibration and reduction were the 
same as those described in Paper I. As before, all the usual photometric points 
have been carefully considered. The calibration spots were imprinted on 
clear glass areas of the plate, using the same tube sensitometer, and, as usual, 
the consistent absence of microphotometer overrun at their edges gives an 
objective demonstration of the complete absence of Eberhard effect. Character- 
istic curves for the same emulsion and exposure time were often completely 
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indistinguishable over the range used in the reductions when different plates 
were compared, and this applied even when the intensity changer on the sensito- 
meter had been used to raise or lower the intensities of the spots as a group. 
The microphotometer work was carried out on the Radcliffe instrument, using 
the same technique as before, the analysing area being a small rectangle with 
effective dimensions of the order of one-tenth of a millimetre, and runs being taken 
parallel to each other with a separation of a quarter of a millimetre. It was not 
thought necessary to make the check previously applied of repeating the whole 
set of runs in perpendicular directions. As before, runs on to clear glass were 
taken every millimetre. The tracings were measured and reduced, the individual 
isophotes being assigned a weight according to the system described in Paper I. 





Imm on plate 


= 2275 


Fic. 1.—IC 1459. 


The final results for the four elliptical nebulae are illustrated in Figs. 1-4 and 
photographs of them are reproduced in Plate 7 A (a), (b), (c),(d). The dotted 
curves in the figures represent the Hubble distributions of best fit. 


44 
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The three remaining nebulae which have, sometimes rather doubtfully, been 
regarded in the past as elliptical are NGC 5253, NGC 7049 and NGC 7213. The 
first has been discussed in a note published separately (2). The other two are 


Imm on plate Imm on plate 
=22°°5 : 
=22°:5 
Fic. 2.—NGC 4976. Fic. 3.—NGC 5102. 


Imm on plate Imm on plate 
Lccevicememaprinnnesl 


222-5 =22°°5 
Fic. 4.—IC 4296. Fic. 5.—NGC 6744. 
shown in Plate 7 B, from which it will be clear why these two subjects were 
misclassified. Both are spirals; NGC 7049 has not very well-developed arms but 
shows a well-defined lane of obscuration round the nucleus; NGC7213 is a 
little more open in its structure. 
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(a) IC 1459 (b) NGC 4976 
1o3a-O, 15 ™ I.A.Z. 10™ 


(c) NGC 5102 (d) IC 4296 
1.A.Z. 30™ I.A.Z. 30™ 


NGC 7049 NGC 7213 
1oza—-O 15 ™ roza—-O 15 ™ 
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The observational material discussed in this paper is completed by a long 
series of plates of the nucleus of the spiral nebula NGC 6744, of which two photo- 
metric plates have been selected for examination by the methods described above. 
The data are as follows :— 


R.A. (1950) Dec. (1950) Mag. Plates 
NGC 6744 19% 05™-o — 63° 56’ 106 A 1623 103a-O —s-115™_~ =1951 Sept. 24 
A 1624 1o3a-O_—s «15™_~ =1951 Sept. 24 


The reasons for the inclusion of this nebula are as follows: the fact that certain 
spiral nebulae possess nuclei which strongly resemble elliptical nebulae has been 
remarked before, but cases suitable for detailed investigation are rare. Such 
nebulae must have a large apparent diameter, and be so oriented that the nucleus 
is open to plain view. Inthe sky south of the Pretoria zenith, there are only about 
ten nebulae which have an apparent diameter greater than five minutes of arc, and, 
of these, only one, NGC 6744, has a nucleus which resembles an elliptical nebula. 
Fortunately this system possesses arms of a very open structure, apparently 
viewed perpendicular to their plane or nearly so. ‘They appear to be defined by 
giant stars, or by bright knots not readily distinguishable from stars, as well as by 
patches of luminosity. If the spiral arms were absent, the nucleus would be 
classified as an elliptical nebula. ‘The nucleus has been investigated by the usual 
methods, and, although the presence of the arms makes the photometry somewhat 
less precise and involves points of interpretation of detail not met with in ordinary 
elliptical nebulae, there seems no reason to doubt the general correctness of the 
results. ‘These are shown in Fig. 5, and it will be seen that they conform to 
Hubble’s rule. ‘The nebula asa whole, photographed with a long exposure, is shown 
in Plate 8 and, on a larger scale, the nucleus, photographed in the blue and red, is 
shown in Plate 9. 


Before proceeding to a discussion of these observational results it is necessary 
to interpolate some remarks on Hubble’s rule, not indeed by way of amplification of 
what was said in Hubble’s original paper (3), but by way of emphasis on certain 
points which do not seem to have been adequately appreciated by later workers. 
In Paper I this rule was quoted in the form: 


m=m,+5log(r/a+1), 
where, clearly, m, is to be interpreted as the magnitude corresponding to the surface 
brightness at the centre of the nebula. The formula has been criticized on these 
grounds: firstly, that it is difficult of physical interpretation; secondly, that it 
contains too many disposable constants; and thirdly, that when integrated to 
infinity it gives an infinite result for the total brightness of the nebula. None of 
these objections seems to have much weight. As Hubble himself pointed out, 
the value of ‘‘a’’ varies along different radii, and must do so in order to make the 
isophotes elliptical, but there is no unnecessary freedom here, only just enough to 
allow nebulae of different ellipticities to be comprehended in the same formula. 
The method of fitting throughout this work has been the following: 
If we write A =10-"e® and X = 10”, we readily transform Hubble’s rule to 
AX=r/a+1 
and can, from the observations, for both the major and the minor axes of a nebula, 
compute the values of X for given values ofr. ‘The problem, then, is to find, ona 
plot of X against r, the two straight lines which fit the data and which pass through 
44* 
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the same point on the X-axis. In fact, within the limits of observational error the 
lines are straight and do pass through the same point on the X-axis. 

We conclude that the observed distributions are of the type required by 
Hubble’s rule, that no violence is done to the observational results by fitting 
the formula to them, and that the formula is not so general that the fact of 
conformity with it is a mere triviality. The definition of the parameter “‘a”’ in 
physical terms can then be made in several obvious ways which have no lack of 
physical significance. 

When all is said and done the formula does not claim more than empirical 
justification. ‘The objection that it gives an infinite total brightness is of a type 
which has beencommoninastronomy. There isnot the slightest reason to suppose 
that elliptical nebulae do extend indefinitely, and certainly no one has pushed his 
observations out sofar. ‘lo maintain this as an objection is to miss the importance 
of Hubble’s rule: in it—or indeed in any other empirical formula which fits 
the results within the observational error (i.e. which is indistinguishable from 
Hubble’s rule over the range of its application)—-we have a means of recognizing 
the presence of a certain type of brightness distribution. We need the rule not for 
mathematical manipulation but as a criterion of classification. 


The statistical position is rather satisfactory, for the fifteen nebulae studied by 
Hubble (3) and the ten true elliptical nebulae described in Paper I and inthe present 
work together include practically all the elliptical nebulae in the sky down to the 
given limits of brightness. ‘That is to say, all the elliptical nebulae which can be 
studied in detail with the largest telescopes now available have been so studied, and 
the two sets of results together permit generalizations to be made which are not 
open to investigations of more limited scope. 


A second digression is desirable before proceeding to discussion. Hubble, in 
a private communication*, remarks: ‘‘ You may be interested in my notes on old 
plates of the five nebulae [described in Paper I] which can be reached with the 
100-inch. ‘The plates are poor in general, because of the low altitudes. NGC 1344, 
1399 and 1404 are noted as normal E with ellipticities roughly estimated as 4 to 5, 
I to 2, and 1totr. All three, together with NGC1380, are considered to be 
members of the South Lobe of the Fornax Group (1344 somewhat doubtfully), 
and Humason has velocities of NGC 1380, 1399 and 1404. ‘The modulus of the 
group is provisionally estimated (from half a dozen velocities) as of the order of 
m— M = 27-4.” 

““ NGC 1380 is noted as not E, but So or very nearly Sa. It shows a central 
lens differentiated from a surrounding envelope which (the envelope) has a vaguely 
thickened rim. The choice between So and Sa depends upon whether or not the 
rim is uniform or broken down into spiral arcs. My plate (of 1938) suggests 
incipient breakdown, but not very convincingly... .” 

“* NGC 1291, also from a 1938 plate, is notedas SBo. Itshows anormal central 
lens about 50 inches in diameter, with a well defined bar running across it from 
E of N to W of S, in agreement with your tilted axes forthe innerisophotes. ‘There 
is no ring, so the lens shows the familiar pattern of 4203, 42607, etc. The inter- 
pretation of the outer regions beyond the lens is rather dubious, but if one is 
‘ sufficiently generous to his imagination’ he may suggest that the extensions W of 
N to E of S form a wide faint bar across a faint tilted ring, perhaps 3'-4 x 2’-0*. 
Such miniature patterns—both spirals and bars—are occasionally found in the 
central lens of SB nebulae. .. .” 


* Quoted by kind permission. 
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“The classes So and SBo are now well established. They fill the range 
between E 7 and nebulae with actual spiral patterns, normal or barred as the case 
may be. I do not recall any normal E nebulae (my own classification) in which 
obscuring patches or additional structural patterns can be detected, with the 
exception of NGC 3115 (in which the 200-inch plates show a couple of small spots 
near one end of the major axis). Otherwise the only evidence of finely divided 
material is the emission at A 3727 found in the nuclei of about 20 per cent of the E 
nebulae whose spectra have been recorded. Patterns of obscuration, indicating 
the presence of dust, are first seen in So, SBo.... ‘Transition cases and, in parti- 
cular, the differentiation of a lens from the envelope as the next step beyond the 
E7, seem to offer a promising field of investigation.” 

It will be seen that in the list of objects under review in Paper I and the present 
work we have a series of nebulae which range from types universally accepted as 
elliptical, through certain transition stages where very faint arms of very rudimen- 
tary type exist, to the case of NGC 6744 with a fully developed set of spiral arms 
and a nucleus of elliptical type. In all these cases the nucleus alone has in the 
past been classed as elliptical, with the exception of NGC 6744, and even here, it 
seems safe to say, the classification in the absence of the arms would have been E. 
In all these cases there is a close conformity with the Hubble type of distribution of 
surface brightness. 

It seems pertinent to enquire what meaning we can attach to the term “elliptical 
nebulae’. ‘I'his has developed as our knowledge has increased, just as the term 
“globular cluster”, for example, has developed. Originally meaning no more 
than the ordinary linguistic significance of the words which compose it, this phrase 
has now acquired far deeper implications. ‘The simple interpretation as ‘‘a 
cluster of stars of globular form”’ has been transcended, and other properties such 
as the possession. of a particular type of density variation, the presence of variable 
stars of a certain type, and the conformity of the star population with certain types 
of distribution are required of such a system. We may even deny membership 
of the class to a system which would readily qualify for the designation as “a 
cluster of stars of globular form” if it does not satisfy these other requirements. 
We might distinguish the class of objects so defined by capital letters, and refer to 
them as “ Globular Clusters”: nor is this distinction merely trivial, for in setting 
up this definition, we almost automatically infer that conformity with it implies not 
only a similarity of structure and properties, but an identity of evolutionary develop- 
ment. It is important to emphasize the change which has taken place because it is 
apt to occur without its being realized. ‘To borrow a biological analogy, we have, 
in setting up the definition of a Globular Cluster, defined a certain phenotype and 
assume that identity of phenotype implies identity of genotype. 

What then do we mean by an “elliptical nebula’? Until recently the sole 
ground of classification was appearance, supplemented by radial velocity data 
which proved that the object was extra-galactic. ‘The term “ 
meant an extra-galactic system of elliptical shape with a smooth gradation of 
intensity from the centre to the periphery and little or no surface detail. ‘The 
ellipticity might have any value, and deviations from ellipticity towards a lens 
shape were permitted. Until recently no elliptical nebula had been resolved into 
stars, but the work of Baade (4) suggests that nothing but the limits of present 
telescopes stands in the way of the resolution of all ellipticat nebulae, and that, were 
they so resolved, they would be found to comprise a special type of stellar population 
(Population II). There is nothing to show that this is not correct, but it is an 

* 


elliptical nebula” 
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illustration of the way in which conceptions develop almost unnoticed that, in spite 
of the limitations of the observational evidence, it is generally assumed that the 
populations of all elliptical nebulae are of type II. Lastly, the work of Hubble and 
the present papers show that obedience to Hubble’s rule of distribution of surface 
brightness is a characteristic of elliptical nebulae. We have now gone as far as our 
observational evidence will carry us and have tentatively formulated a definition 
of an “‘ Elliptical Nebula”’, implying, as in the case of the term “‘ Globular Cluster”’, 
an identity of evolutionary origin. 

There are several interesting points about this definition. It denies classifi- 
cation as “Elliptical Nebulae”’ to objects such as NGC 1553 (classified as So) 
which are certainly elliptical nebulae. What, in view of Hubble’s remarks, we 
should call the nuclei of NGC 1380 and 1291 would seem to qualify as “ Elliptical 
Nebulae”’: there seems no difficulty about NGC 1380, while in the case of 
NGC 1291 the only question is whether we are to overlook the development of 
twist in the isophote axes. It is especially interesting that even in cases where 
the eye construes the density distribution into a bar structure, the isophotes 
remain sensibly elliptical. Finally, in the nucleus of NGC6744 we have a 
structure which retains the Elliptical Nebula characteristics right into the spiral 
sequence. 

If speculation is to be eschewed this paper should end here, but it may be 
profitable to discuss briefly the evolutionary implications of a situation in which it 
is possible for the Elliptical Nebula characteristics to be retained by nebulae and 
nuclei over a wide range of types of object. Most discussions of the evolution of 
nebulae have been based on dynamical considerations: the evolution of a bar 
structure under the influence of rotation (§); evolution as described by non- 
Newtonian mechanics (6); the condensation of primordial matter into stars (7) ; 
and so forth. 

The objects we have been discussing are so close to each other in structural 
type, and yet show such a sequence of progressive differentiation, that it seems fair 
to infer that they represent stages of nebular evolution which are rather close 
together. One problem is the temporal order in which they are arranged. ‘This 
part of astronomy is full of such conceptions as the condensation of stars from 
primordial gaseous or finely divided matter, and yet at the same time the elliptical 
nebulae are regarded as preceding the spirals in the evolutionary sequence. ‘These 
ideas are difficult to reconcile, for the elliptical nebulae are characterized by an 
almost complete absence of dust and gas and presumably by a complete resolution 
into stars. 

On the other hand, it is difficult to see how a configuration such as an elliptical 
nebula can be the last stage of a process in which spiral arms have condensed into 
the nucleus. Considerations of angular momentum alone would seem to provide 
an insuperable barrier. In addition a widely dispersed system such as an open 
spiral involves relative velocities of the component parts of quite a high order, and a 
gravitational field at the periphery which is relatively feeble. It seems hard to see 
how the dispersive tendencies can possibly be overcome by any forces tending to 
aggregate the material together. 

A way out of this impasse would perhaps be provided by the assumption that in 
elliptical nebulae there are forces which produce two effects: first, a tendency to 
differential evolution of stars causing the production at a relatively rapid rate of a 
class of blue giants, and second, the ejection of these giants from the nebula.* 


* Or, of course, the ejection of objects destined to evolve into blue giants. 








No. 6, 1952 Surface photometry of southern elliptical nebulae 613 


The explanation for these processes would have to be sought, partly in knowledge 
of the evolutionary processes of individual stars and partly in knowledge of the 
dynamical evolution of systems like elliptical nebulae. If we assume these two 
effects we can see how evolution would proceed. Rotational effects would arise: 
arms would appear, marked out by giant stars, and possibly an interstellar gas 
cloud would make its appearance in consequence of the ejection of lighter atoms 
from the surfaces of the giant stars by radiation pressure. In some ellipticals, as 
has been remarked by Hubble, some interstellar gas is already present. The 
nucleus of the nebula would remain as a quasi-elliptical nebula, and it would 
consist of a star population representing the residue (but possibly itself somewhat 
modified by stellar evolution over a long period of time) after the removal of the blue 
giant stars. ‘That is to say, there should arise a differentiation in type of stellar 
population between the arms and the nucleus. Moreover, if we have assumed a 
differentially rapid rate for the evolution of the stars forming the spiral arms it 
seems reasonable to suggest that the arms are a relatively temporary phenomenon, 
and might disappear completely leaving behind an elliptical nebula reduced in size. 
This is not a necessary consequence of the ideas suggested, but if it should be 
found that various classes of elliptical nebulae exist having sharply differentiated 
sizes it would be an observational confirmation of them. 

In the observational part of this work I have been much assisted by Mr J. 
Churms of the Royal Observatory staff who made the tracings of several of the 
plates using the Radcliffe Observatory microphotometer, and by Mr R. D. Cherry 
of the University of Cape Town, a volunteer worker at the Royal Observatory 
during 1952 January and February, who carried out most of the reductions of the 
tracings. ‘The facilities made available at the Radcliffe Observatory during the 


latter part of this investigation have been afforded under the recent cooperation 
agreement between the Royal Observatory at the Cape and the Radcliffe Observa- 
tory at Pretoria. 


Royal Observatory, 
Cape of Good Hope: 
1952 June 13. 
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Summary 


A set of tables is given from which it is possible to obtain quickly the 
relative transition probabilities in band systems of diatomic molecules if the 
spectroscopic constants r,, w, and x, are known. T'wo methods of using the 
tables are described. One is applicable if the electronic transition causes a 
small or moderate change in the equilibrium internuclear separation; and 
the other is applicable if the change is moderate or large. The first method is 
not significantly less accurate than the Morse approximation. ‘Though the 
second method does not yield results of such high precision, it should be 
useful for many purposes. 





1. Introduction 
One of the difficulties in the full exploitation of quantum mechanics to 
problems of astrophysical interest is that many of the calculations required are 
elementary. For example, it is often desirable to know the relative transition 
probabilities in band systems of diatomic molecules (cf. Swings (1)). These 
probabilities can be computed quite accurately. All that it is necessary to do is 
to derive the vibrational wave functions and from them evaluate the relevant 
nuclear overlap integrals. However, because of its straightforward character 
the investigation makes little appeal to the theorist, and consequently only a 
small proportion of the band systems of importance in astrophysics have been 
treated. ‘The purpose of this paper is to rectify the position at least partially 
by presenting a general set of tables from which results for a wide range of cases 
can be read. To ensure that the procedure to be followed is clear, some examples 
are worked through in full in Sections 2.1.1 and 2.2.1, which have been made 
as self-contained as possible for the benefit of readers concerned only with the 
practical application of the tables. 
2. Theory 
The electronic part of a molecular transition matrix element depends on r, 
the internuclear separation, but provided the dependence is not too marked 
it is a fair approximation to ignore it (2), and to take the spontaneous transition 
probability between the vibrational levels v’ and wv”, of two electronic states 
A and B, as 
Cv(v', o")(v', o"), (x) 
where @ is a constant for a given band system, v(v’, v”) is the frequency of the 
radiation emitted, and p(v’, v”) is the nuclear overlap integral 
Jf P(r|v’, A)P(r|v", B) dr, (2) 
P(r|v’, A) and P(r|v’,B) being the normalized vibrational wave functions. 
The important dimensionless quantity p(v’,v”)? does not seem to have been 
assigned any convenient name. Because of its connection with the Franck—Condon 
principle we will refer to it as the Franck—Condon factor of the particular band 
concerned. 
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In detailed calculations on a limited number of band systems the wave 
functions, P, are best represented by the Morse approximation (3), but 
unfortunately such representation is not convenient in a more general 
investigation. ‘The researches of Gaydon and Pearse (4), Pillow (§), and 
Nicholls (6) show that the simple harmonic approximation of Hutchisson (7) 
is satisfactory if allowance is made for distortion. However, the procedures 
they introduce are also unsuited to systematic tabulation. 

2.1. Method A—change of equilibrium internuclear separation small or 
moderate.—'Vhe basic problem is to take proper account of the main molecular 
parameters on which the Franck-Condon factors depend without making the 
presentation of the results unacceptably lengthy. Now to obtain sufficiently 
accurate wave functions it is necessary to include terms involving the reduced 
mass, the equilibrium internuclear separations, the vibration constants and the 
anharmonic constants, or, in the standard spectroscopic notation (8), 4, 7., We 
and x, respectively. ‘Thus in addition to the vibrational quantum numbers 
there are seven essential parameters. ‘The tabulation of a quantity dependent 
on so many is usually impracticable but as will be seen it can be carried out for 
Franck—Condon factors by expressing their square roots as the sum of three 
parts, one being a function of two variables and the others being similar functions 
each multiplied by a further variable. 

It will be convenient to distinguish both the parameters and vibrational 
quantum numbers of the electronic states concerned in the transition by fixing 
the subscript S to all quantities associated with the state having the smaller w, 
value, and the subscript L to all quantities associated with the state having the 
larger w, value. Thus r,s, w,g, X,g and vg are the parameters and vibrational 
quantum number of one state; r,,, w,z, X,z, and vz, are those of the other state ; 
and by definition 

Wg SWez- (3) 
The nuclear overlap integrals will always be written as p(vg, vz), that is, the 
vibrational quantum number of what will be called the S state will be placed 
before that of the ZL state. Particular care should be paid to this necessary 
convention as it is not consistent with the common practice of placing the 
vibrational quantum number of the upper electronic state before that of the lower. 

We will first take x, and x,, to be zero, that is, regard the molecule in both 
electronic states as a simple harmonic oscillator. Hutchisson (7) has given 
formulae for the overlap integrals, g(vs, vz), corresponding to this model. ‘These 
are rather cumbersome and need not be displayed here, it being sufficient to 
note that from them it is possible to express q(vg, vz) for any chosen vg and vz 
values as a function of w,;/w,g and a quantity A defined by 


A = 2r(uw,z/h)"* {rg —Ter}; (4) 


h being Planck’s constant. We will always measure the reduced mass, p, on 
the chemical scale (O=16), the vibrational constants, w,, in cm}, and the 
equilibrium internuclear separations, r 
becomes 


, in angstroms (10°8cm), so that (4) 


€ 


A =O°172 (jew,z)"* {7,5 —T-z}- (5) 


It may be remarked that (vg, vz) is either a symmetrical or an anti-symmetrical 
function of A according to whether (vg+v,) is even or odd. ‘This makes it 
sufficient to consider positive A only. Cases of negative A are of course rare 
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because of the Birge-Mecke empirical rule (8) that the product r,2w, has about 
the same value for all electronic states of any one molecule. 

The double entry tables * at the end of the paper give q(vg, vz) as functions 
of w,,/w,gs and A for all combinations of vg and vz, such that 


Ugt Up S4 
and, in addition, for Ug=0, v0, =5,6, 
U,=0, Ug=5,9. 
They cover the ranges Wp /Wes = 1-0 (0°1) 1°6, 
A=0°0 (0:2) 5:0, 
which are ample for most transitions. The signs of the figures in the tables 
were determined by the convention of making all wave functions positive near 
the origin. 
To obtain a better approximation to p(vg, vz) than is provided by q(vs, vz) 
it is necessary to allow for the fact that the potential of the electronic state of a 
molecule differs from the potential of a simple harmonic oscillator. ‘The matter 
has been investigated by Hutchisson (g) but his treatment is too complicated to 
be used without modification. 
In the region of the minimum the difference between the potentials may, 
with sufficient accuracy, be taken as 


Rrr2 \ 3/2 
“a “ ( +) w,> o air i, Px (6) 


c being the velocity of light (8). Since this is small we may use first-order 
perturbation theory to expand the true wave functions, P(r|v), in terms of the 
corresponding simple harmonic wave functions, Q(r|v), finding after some 
trivial algebra 

»| {x(v—1)(v—2)}? 
P(r|v)=Q(r|v)+x2 f ( xe } 





ne govt 
O(r|v—3)+ =o Brie —1) 


3(v + 1)9? {((v+1)(v+2)(v+3)}? 
eae coment ok inthe’ 2am QO(r|v+3) |. (7) 
4 I2 
There is, of course, no first-order change in the energy levels. It may be noted 
that the dimensionless anharmonic constant, x,, is usually of order 107. 
Substituting (7) in (2) and neglecting a term involving the product x,," x,,¥? 
we immediately obtain the approximate formula 


P(g, Vz) = (Ug, Uz) + Xes'"S(Ug, Up) + Kez" "(Vg Vz) (8) 





with 
{vg(vg—1)(vg—2)}? 30692 
5(Ug, vz) = |! s\ , A . qvg < 3 Uz) + : Yvg- I, Uz) 
3(vs + 1)°” {(ug + 1)(vg + 2)(vg + 3)}"? 
— Es ges +1,01)- 4 S s (vs } 








vst: 
and 


fo. (9, ~ 310, —2))"2 » 9/2 
vs, v;) = | I)(v, —2)} = 





12 4(Ug, Uz — 3) + q(vg, Uz —1) 


3(vz, +1)*? {(v, + 1)(v_ + 2)(vz, + 3)} 
se aa — 909 es+ 1)- + ms £2) z+3)} (Us, Vz, + 3) |- 
( 





10 





* Thanks are due to the departmental computer, Miss Irene E. Montgomery, for her able 
assistance in the preparation of these tables. 
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In order to facilitate the use of method A a set of double entry tables of 

$(vg, vz) and (vg, vz) covering all vg and v, satisfying 

Ugt Vy <3 
is presented. The range of w,;/w,s is from 1-0 to 1-6 as in the main tables; 
that of A is reduced to from 0-0 to 2-2, but is still sufficient for many transitions 
of astrophysical interest. 

Inspection of the tables shows that the anharmonic correction is relatively 
small when vg and wv, are equal, and the difference between x, and x,,, is not 
very great. A useful consequence of this is that the simple harmonic approxima- 
tion may generally be used for the (2,2) band (for which the anharmonic 
correction has not been computed). 

Before concluding this section it should be pointed out that the expression 
for the wave function P(r|v) given in (7) should strictly be multiplied by the 
normalization factor 


[ ae I)(v-—2) ge O(v+1) (v+1)/v+ Met} | “1/2 
1+%.§—— + St + - ea 
144 16 16 144 
Because of the smallness of x, this is actually not very different from unity. 
The experience so far gained in the method indicates that it is best omitted. 
2.1.1. As an example we shall consider the calculation of the Franck—-Condon 
factors for the bands of the First Negative system of nitrogen, B *X,*+— X #X,*. 
The reduced mass » is 7-004 on the chemical scale. For the spectroscopic 
constants we will adopt the values used in some earlier work (2) on the Morse 
approximation.* ‘These are as follows: 


B 2° 


~u> 





r,=1°074 A, 


w,=2418 cm}, _ - P 
giving x,//2=9:3 x 1077; 
w X, = 20°95 cm"|!, 


ry 25° + == T° 
&% °Z,"*, r,=11I7A, 


. yo 1 
wW,= 2207 cm”, hg 2 
giving x, 


wx, = 16°14 cm", 
Since w(X 2X 5*)<w,(B #Z,*) 


the subscript S will be affixed to the parameters and vibrational quantum 


numbers associated with the X *X,* state, and the subscript L to those associated 
with the B?X,,* state (cf. (3)). We thus have 

® -1/Weg = 2 418/2 207 = 1°00 
and from (5) A=0-172,(7-004 x 2 418)"*{1-117 — 1-074} =0°96,. 
Consider a particular band, say the band for which the vibrational quantum 
numbers, vg and wz, of the X*X,* and B#X,,* states are 2 and I respectively. 
Using linear interpolation we find from the tables that when w,,/w,g is 1°0Qg, 


and A is 0-964, then 
S WW q(2, 1) = +0°OI9, 


(2, 1) = — 0°48, 
(2, 1) = — 0°36. 





* The values recommended in Herzberg’s recent book (8) are slightly different. 
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Hence by formula (8) we obtain * 
pl2, I)= +0°61,+(0-085)( — 0:48) + (0°093)( — 0°36) 
oa 


and so the Franck—Condon factor 


p(2, 1)? =0-29. 
Results obtained for other bands of the system are displayed in Table I together 
with the corresponding results obtained from the simple harmonic and the 
Morse approximations (2). ‘lable II gives a similar comparison for the First 
Positive system of nitrogen, B 311 ,—A*%2X,*, for which w,z/w,s, A, x,g'" and 
x,y)" are I°18,, 1°574, 9°8 x 10°? and 9-1 x 10° respectively. ‘Taking the Morse 


TaB_e I 
Franck—Condon factors for bands of First Negative System of nitrogen 
Figures in round brackets, ( ), are based on the simple harmonic approximation, those in curly 


brackets, { }, on the present approximation (A), and those in square brackets, [ ], on the Morse 
approximation. 





(X *Z yt) 


(-30) {:26} [-26] (06) {-08} [-o7] | (00) {-o2} [-o2] 
(:20) {-20} [-22] | (-37) {29} [-29] 
(-30) {-40} [41] | (:03) {:03} [-o4]T 























Tasie II 
Franck-Condon factors for bands of First Positive System of nitrogen 
Figures in round brackets, ( ), are based on the simple harmonic approximation, those in curly 


brackets, { }, on the present approximation (A), and those in square brackets, [ ], on the Morse 
approximation. 





(A *2,,*) 


(B *IIy) 
Cp 
(-32) {-31} [-32] | (-39) {-32} [32] | ¢ . & (06) {+12} [ 
(-33) {40} [41] | (or) {or} [-o1] 
(:20) {+22} [-21] | (-10) {+17} [-19] 
(09) {07} [ J 


























approximation as standard, we see that the present approximation is superior 
to the simple harmonic approximation, and indeed does not cause any significant 
inaccuracy. ‘lo verify this conclusion some other systems which have been 
treated by the Morse approximation were investigated. ‘The detailed results 
need not be presented, for in each case the essentially complete success found for 
the First Negative and Positive systems was repeated. 

* Attention is again called to our convention of placing vg before vz in p(vg, vz) irrespective 
of which electronic state lies highest. 

t See comment in text Section 2.1 concerning the (2, 2) band. 
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2.2. Method B—change of equilibrium internuclear separation moderate or 
large.—Cancellation effects in the overlap integrals are unimportant if the 
difference between the equilibrium internuclear separations of the electronic 
states is sufficiently great. Useful estimates of the Franck-Condon factors can 
therefore be derived with comparatively crude wave functions. ‘The simple 
harmonic approximation is indeed often adequate. It can readily be applied 
with the aid of the q(vg, vz) tables which, as has already been noted, cover a 
wider A range than do the s(vg,v,) and U(vg,vz,) tables. However, increased 
accuracy can be obtained with but little additional labour. 

The major contribution to each overlap integral arises from the region 
between the equilibrium internuclear separations of the two electronic states. 
Consequently it is natural to attempt to make partial allowance for the asymmetry 
of the molecular potentials by introducing in the simple harmonic approximation 
to each wave function an effective vibrational constant w,’, which is greater or 
less than the actual vibrational constant w, according to whether the electronic 
state concerned is the one with larger or smaller equilibrium internuclear 
separation. Several rules may be devised for determining the optimum value 
of this new constant (8), but to avoid any undue increase in the numerical work 
incurred in the application of the tables we will merely adjust the simple 
harmonic potentials so that at the point midway between the equilibrium 
internuclear separations they become equal to the corresponding Morse 


potentials. Writing 
27? cuw,x,\12 
a= ——|} d, (1 


where d is the modulus of the difference between the equilibrium internuclear 
separations, we readily see that this is equivalent to taking 


pe I—exp( +a) ; 
4 =H] Fa } (13) 


w {1 2 ha}, (14) 


the upper signs referring to the state with the larger r, value, and the lower to 
the state with the smaller. A more convenient form of (14) is 


w, =w,{I + 0°06 d(w,x,u)""}, (15) 


d being in angstroms and the other parameters being in the units already specified. 
It may be remarked that the final results are not critically dependent on the 
precise value of the factor multiplying d(w,x,.)". 

A further improvement can be gained by making the central turning points 
or zeros of the simple harmonic functions coincide with those of the corresponding 
Morse functions. ‘This requires (cf. Pillow (§)) the introduction of an effective 
equilibrium internuclear separation r,’ such that 


(16) 


with 


6e (ee 
Sm“ UC We, 


x, \ 
=4°1I | — '+ 1) angstroms. 
pin 
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The method of finding the Franck—Condon factors based on the modified 
simple harmonic wave function is clear. Briefly, we first label the electronic 
states S or L so that 

Weg Spy, (18) 

and then use the q(vg, v,) tables as before, taking as the variables w,;'/w,,' and 

A’(Ug, Up) = 0°1729(Hmz')"*{(re5 +55) — (Ten + 51)}- (19) 

In general w,;'/w,s' is less than w,,;/w,g which has the fortunate effect of bringing 
additional band systems within the tabulated range. 

2.2.1. To illustrate the practical procedure we will calculate the. Franck— 
Condon factors of some bands of the Vegard—Kaplan system of nitrogen, 
A *X,+—X'Z,*, which is a typical system involving a large change in the 
equilibrium internuclear separation.* We will consider in particular the 
progression arising from the zero vibrational level of the ground state as 
information on this is of some interest in absorption spectroscopy. 

The values of the molecular parameters are as follows: 

#=7°004, 
A *h,*, r,=1:293A, 
w,=1 460 cm", be 
giving x,=9°5 x 107%; 
wx, = 13°93 cm™, 

r,=1:095 A, 

w, = 2 360 cm~, nee 
giving x,=6-1 x 1073; 
wx, =14°44 cm", 

d =(1-293 — 1095) A=o-1908 A. 
Noting that the equilibrium internuclear separation of A °X,*+ is greater than 
that for X 'X,+, we see from (15) that 

A 8,,*, w,’ =1 460{1 + 0-06 x 0-198 (13°93 x 7°004)"?} cm™ 


=1631 cm", 
and 


X12 ,+, w,’ =2 360{1 — 0-06 x 0-198(14°44 x 7:004)"?} cm=! 
=2078 cm}, 
Clearly A *X,,* is the S state and X 'X,* is the L state (cf. (18)). Thus 


wp /wW.g = 2078/1 631 =1°274, 
and from (17) we see 


9°5 x 10°% 


1/2 
8e=4°11 < ———— v I)=4°0 x 1073 (vg+ 1 
8 Sera (Ug+1)=4 (vs+1), 


5 se 6-1 x1078 )1/2 ee nner 
= A° ee ee Ye = = 2°5 X z va) I 
L=4 7-004 X 2 360 (vy +1)=2°5 x 10-* (vz +1), 





* The main reason for the choice of the Vegard—Kaplan system is that the Franck—Condon factors 
based on the Morse approximation are known from some unpublished calculations by the writer. 
It must be emphasized, however, that the system is one in which the electronic part of the transition 
integral may vary so rapidly with internuclear separation that the Franck—Condon factors may lose 
much of their usual significance. 
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so that substitution in (19) yields 


A’(vg, vy) = 0°1729(7°004 x 2 078)"*{0-198 + 0-004 0(vg+ I) — 0-002 5(vz,+ 1)} 


= 20°7,{0°198 + 0-004 0(vg+ 1) — 0-002 5(vz,+1)}. 


With the aid of the g(vg, vz) tables we can now determine the Franck—Condon 
factors. ‘The values obtained, together with the values based on the simple 
harmonic and Morse approximations, are displayed in Table III, a logarithmic 
scale being used. As can be seen, the present and the Morse approximations are 
in fair accord ; and even the simple harmonic approximation gives results correct 
to well within an order of magnitude (which in some instances is all that is 
required). 
Taste III 
Log, of Franck—Condon factors for bands of the Vegard—Kaplan system of nitrogen 

(X 129+), up=0; (A*2y*),vug= 0 I 2 

Approximation 

Simple harmonic 

Present (B) 

Morse 


2*4o 
2°05 


Department of Applied Mathematics, 
Queen’s University, 
Belfast, Northern Ireland : 
1952 September 30. 
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Description of Tables 

Of the following tables, ten groups relate to transitions in which the combination 
(vg, vz) of the vibrational quantum numbers is (0, 0), (0, 1), (1,0), (©, 2), (2, 0), (0, 3), 
(3, 0), (1, t), (1, 2) or (2, 1). The first table in each of these groups gives the values of 
the simple harmonic overlap integrals q(vg, v,), and the second and third give the values 
of the anharmonic corrections s(vg,v,) and (vg, vy) respectively. ‘The nine separate 
tables give the values of q(vg, v,) for the transitions specified by the combinations (0, 4), 
(4, 0), (0, 5), (5, 0), (0, 6), (6, 0), (1, 3), (3, 1) and (2,2). In reading from the tables it is 
sufficient to employ linear interpolation. 
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q(0, 0) 
1°*3 


0'996 0’990 
‘987 . 982 
‘962 ° “959 
‘921 ‘ "921 
‘866 ° ‘871 
‘B01 : ‘SIO 
. 28 ° "742 
‘650 . 669 
‘571 : 593 
“417 " “445 
"348 ; ‘376 
"285 . 313 
229 ’ °2506 
‘181 . ‘206 
"141 "152 "164 
"107 . ‘125 
‘O81 ° ‘095 
‘060 ' "074 
043 ’ 055 
O31 ‘ ‘O40 
"022 025 029 
‘O15 ‘O15 ‘O21 
"O109 . "O14, 
‘006, ° "000, 

0'0043 . 0-006; 


s(o, ©) 
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qo, 1) 
- : ; 1-3 


“000 "000 0000 
‘140 ‘133 — ‘121 
272 "259 . ‘237 
"388 ‘371 : "340 
“482 462 : "426 
“551 *530 , *493 
"592 573 ; "537 
‘606 ‘SOI 

‘597 "585 

566 “560 

"520 “519 

464 468 

"402 "410 

339 "350 

‘279 ‘291 

‘224 ‘237 

‘175 “188 

"134 146 

"100 "LI 

073 "082 

"052 ‘060 

036 "042 

025 029 

"016 "020 

"O10, O13, 

006, —0-'008, 
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q(1, ©) 
13 


0000 
138 
"270 
"387 
“486 
"562 
613 
638 
“640 
‘621 
"585 
"536 
‘479 
‘417 
‘356 
296 
‘241 
‘192 
‘150 
"115 
086 
063 
046 
032 
022 

O°Ols 


s(1, 0) 
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q0, 2) 


I'2 fe 


0064 0092 
O75 ‘102 
*107 "130 
156 174 
‘217 228 
283 288 

‘347 
‘407 ‘401 
‘453 443 
‘483 ‘472 
‘496 ‘485 

482 
‘471 465 
‘438 ‘435 
"396 "396 
"348 ‘352 

"304 
‘248 ‘257 


sim nO re) 
Anya eh O 6 


S6uuGs 
 ~ 
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q(2, ©) 
1°3 


—0°092 
‘077 
"035 
"030 
"II3 
"205 
‘2907 
"383 
"455 
"509 
"542 
553 
544 
"517 
‘477 
427 
"373 
‘317 
263 
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Op™3 qo, 3) 
W,7/Wog I ‘oO . . I 3 


°'o 0'000 . 0°000 
o'2 — ‘ool . — ‘020 
0°4 -009 ; : 044 
°°6 ‘028 . . 073 
° 063 . . ‘110 
‘II2 , , "155 
‘174 ‘ ‘ 205 
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+312 A : 312 
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425 “415 ; 399 
‘458 < ‘ “426 
473 : ‘ ‘441 
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’ “427 

"402 

368 

*328 

"285 

‘241 

"199 

‘160 

"126 

*097 

073 

—0'054 


coAR N OG Co 


~ NO 


s(o, 3) 
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q(3, 0) 
1°3 


0’000 
-— S83 
"034 
"034 
— ‘O15 
+ +023 
‘079 
"149 
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"305 
‘376 
‘435 
‘477 
“500 
"502 
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‘170 
"132 
‘100 
+0'074 
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v,p=!I q(1, 1) 
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168 : 208 

128 . 164 

127 

"095 

070 

"050 

035 


°° 


on = me Hm Ht Hm mH O 
= HH eH om HH OO 


OF He eH me He mH mw me O 


| 
° 
| 
° 
| 


W.1/%g 1 


°o’o 
o'2 





D. R. Bates, Relative transition probabilities 
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"403 
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Ug =2 

vp=! q(2, 1) 
W,,/%5 10 , . 1°3 
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oo 0*000 . 0*000 
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qo, 4) 
1°3 14 


0'O10 O'OI7 
013 "020 
"020 028 
033 "042 
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“080 ‘OQ! 
‘IIS 125 
‘157 "165 
‘203 ‘210 

‘255 
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‘378 ‘371 
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"339 ‘338 
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‘225 


~~ 00000), 
vO WwOO+ N O 


AEN GO OAS 


NNNNN SH HH 
i) 





No. 6, 1952 in band systems of diatomic molecules 


val ined 


Up=5 G0, 5) 
1/5 1°0 : 13 
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3°0 


5 
° 


Wepl%g 
A 


°° 
-t- 


‘OO 
“OI! 
"033 
069 
‘117 
‘174 
236 
‘298 
"354 
‘398 
426 
438 
‘432 
‘411 
‘378 
‘336 
"290 
‘242 
‘197 

+0156 





D. R. Bates, Relative transition probabilities 


qo, 6) 
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Ug=I 
vU,>=3 qi, 3) 


@,,/@.5 10 ’ : 1°3 
0°000 . +o'158 
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"092 ? : ‘211 
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‘397 "386 
468 ° : "424 
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g(2, 2) 
1°3 1*4 1°5 16 


+0°970 +0'952 +0°930 +0908 
"929 ‘912 "894 ‘873 
"809 ‘800 “788 ‘174 
627 628 626 621 
"407 “419 ‘427 "432 
176 "197 ‘215 "230 
039 O12 + ‘o12 + 034 
"214 "187 ‘160 "134 
"334 ‘311 286 ‘261 
‘392 ‘376 ‘357 ‘338 
"387 ‘382 ‘373 ‘361 
"330 ‘337 ‘339 "336 
235 252 264 ‘272 
"116 "142 *162 ‘179 
008 - 022 - "049 — ‘O71 
123 + 092 + 064 t+ "039 
‘219 ‘192 "165 "140 
"290 268 246 ‘224 
‘333 "319 "303 ‘286 
"350 ‘344 ‘336 "325 
"344 ‘347 "346 ‘342 
"321 "331 "336 ‘339 
301 "313 321 

‘279 ‘292 

"241 256 

+0'201 +O'217 





PLANETARY MODELS OF TERRESTRIAL TYPE 
K. E. Bullen and A. H. Low 
(Received 1952 October 7) 


Summary 

Curves have been constructed showing the variation of mean density 
with mass of planetary models of terrestrial type. ‘The models are taken to 
have the same composition at zero pressure as the Earth, and are based on 
the pressure—density relations of an Earth Model B and on the hypothesis 
that the outer part of the central core of this model is a pressure modification 
of the mantle. ‘Two curves are shown, corresponding to the assumptions 
that the Earth’s inner core is caused by chemical and pressure changes, 
respectively. The curves are compared with Ramsey’s original curve 
which was based on an Earth Model A and assumed the density distribution 
to be determined essentially by pressure alone. 





In the course of his work on the constitution of terrestrial planets, W. H. 
Ramsey*, on certain assumptions, obtained a curve giving the mean densities 
p of possible planets of terrestrial type in terms of the masses M over the range 
— 2-0 < log ,9(M/ Mz) <0, where My is the mass of the Earth. ‘The assumptions 
were essentially that : 

(a) Planets of terrestrial type are of the same primitive composition. 

(b) The density distribution throughout the Earth and any of the planets is 
broadly determined by pressure rather than by changes of chemical composition. 

(c) The density distribution in the Earth is that of the Earth Model A set 
up by one of the present writers. + 

Ramsey has compared his curve with observational data on the masses and 
diameters of the Moon, Mercury, Mars and Venus. ‘The agreement is excellent 
with the Moon, but this only checks the compatibility] of the Moon’s mean 
density with the density in the upper part of the region B§ of Model A. It is 
far from certain that the Moon is a typical terrestrial planet, and theories other 
than the assumption (a) would lead to the agreement found. Thus the Moon 
does not provide an important test of Ramsey’s assumptions. 

The uncertainties in the data quoted by Ramsey for Mercury were so high 
that Mercury could not be used as a test. Urey|| has since pointed out that 
more recent data of Rabe would yield a mean density for Mercury of order 
1-8 g/cm* higher than the most probable value quoted by Ramsey. Urey has 
interpreted the new data as strong evidence against the assumption (4) above, 
but one of us{] has put forward further reasons to show that Mercury should not, 
on present evidence, be used as a test; the exceptional features of Mercury will 
be discussed in a further paper. 

* W. H. Ramsey, M.N., 108, 409, 1948. 
K. E. Bullen, Theory of Seismology, Cambridge, 1947. 


t 
t H. Jeffreys, M.N., Geophys. Suppl., 4, 50, 1937; 4, 62, 1937. 
§ K. E. Bullen, loc. cit., p. 218. 


I 
H. C. Urey, The Planets, Their Origin and Development, Yale, 1952. 


” 


“| K. E. Bullen, “‘ Cores of Terrestrial Planets ’’, Nature, 70, 363, 1952. 
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In the case of Venus, Ramsey showed that the assumptions (a), (), (c) led to 
fairly good agreement with the observational data, whereas Jeffreys* had shown 
that the replacement of (5) by the assumption that the Earth’s whole central core 
is chemically distinct from the mantle is incompatible with the assumption (a). 
Ureyt considers that Venus is so similar to the Earth that it can be made to fit 
any theory that accounts for data on the Earth. ‘The present writers, while not 
in full accord with Urey on this point, nevertheless hold that the measure of 
agreement obtained by Ramsey in the case of Venus is far from sufficient to 
substantiate the set of assumptions (a), (5), (c). 

This leaves only Mars for crucial testing purposes, and the agreement obtained 
by Ramsey was much less satisfactory for Mars than for Venus; his diagram 
shows a discrepancy of 10 per cent between the calculated and observed values 
of the mean density of Mars. He would rectify this discrepancy by modifying 
the Earth Model A by increasing the density gradient steadily throughout. ‘This, 
however, is an ad hoc procedure, and considerably reduces the value of Mars 
as a crucial test of (a) and (0). 

On grounds quite independent of astronomical data, one of the present 
writers set up a second Earth Model BY, and was able to show§ that, in conjunction 
with the assumptions (a) and (b), this model led to agreement with data on Venus 
comparable with that found by Ramsey and, at the same time, rectified the 
discrepancy found by Ramsey with Mars. 

In these circumstances the present writers have thought it worth while to 
repeat Ramsey’s work, taking the assumptions (a) and (4), but using the Earth 
Model B instead of Model A. ‘The result is shown in the broken curve (I1) of 
Fig. 1. 

A second curve (1), shown by an unbroken line in Fig. 1, was prepared on 
the same assumptions as in the preceding paragraph, except that (6) was modified 
by assuming that the inner core of the Earth arises from changes of chemical 
composition and not as a consequence of the pressure. The effect on the 
calculated mean densities of thus modifying (4) is shown separately in Fig. 2. 

For the sake of comparison, the original curve (III) of Ramsey is shown, by 
a dotted line, also in Fig. 1. The points plotted for Venus, Mars and the Moon 
correspond to observational data of Russell, Dugan and Stewart.|| Mercury 
has not been included for reasons already stated. 

Scrutiny of Fig. 1 shows that Ramsey’s curve fits the Moon best. ‘The 
discrepancy between the curve II and the data for the Moon arises from the 
density changes at a depth of 80km in the Earth Model B, but is not significant 
since an arbitrary feature in assigning these density changes makes the curve II 
least reliable when p is in the vicinity of 3-3 g/cm*®. ‘The curve I, however, makes 
it unlikely that the Moon has a core that would correspond even to the Earth’s 
inner core. 

In respect of Venus, Fig. I confirms the conclusion in Paper 3 favouring the 
presence of chemically distinct inner cores in the planets when Model B is used. 
The present calculations have incidentally brought to light a slip in Paper 3, in 
which the value 6312km in the last line of page 256 should be increased to 
6399km. ‘The effect of making this correction is, owever, to strengthen the 

* Loc. cit. tLoc. cit. 
¢ K. E. Bullen, M.N., Geophys. Suppl., 6, 50-59, 1950 (to bx _cved to as Paper 1). 
§ K. E. Bullen, M.N., 109, 688-692, 1949; 110, 256-259, 1950 (to be referred to as Papers 2, 3, 


respectively). 
|| Astronomy, Ginn and Co., 1945. 





No. 6, 1952 Planetary models of terrestrial type 639 


conclusion of Paper 3 by widening the gap between the observational data and 
the results calculated on the assumption that Venus has no inner core. 

Both I and II fit Mars far better than Ramsey’s curve III, and again a chemically 
distinct inner core is favoured, a result which has been further supported (Paper 2) 
by observations of the ellipticity of Mars. While the agreement of the quoted 
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Fic. 1.—Curves (based on three sets of assumptions) showing the variation of mean density with 
mass of p'anets of terrestrial type. 
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Fic. 2.—Proportionate increase in calculated mean density of a planet when it is assumed that 
the Earth’s inner core is due to a change of chemicai composition instead of to pressure, 


observational data on Mars with the curve I is remarkably good, it needs to be 
mentioned that more recent observations* leave some considerable doubt as to 
the best value to take for the diameter of Mars. Urey inclines to the view that 
the mean density of Mars may be appreciably greater than that shown in Fig. 1. 
Should this view be substantiated, the astronomical evidence in support of the 
assumption (6) or a modified form of (6) would be much reduced. But, in the 
meantime, the curve I provides a useful model. 

In preparing the curves I and II, ten planetary models were worked out in 
each case. ‘The crustal layers were taken as in Paper 1 and, as in Paper 2, the 
inner core of the Earth was taken to consist of the region G together with half 

*See Urey, loc. cit. 
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that of F. In preparing II, selected central pressures were taken and numerical 
integration outwards applied in conjunction with the pressure—density relation 
of Model B. The procedure in the case of I was by successive approximation 
from the results in II, following the method in Paper 2 for determining the 
pressure—density distributions in the inner cores. 
TaBLe I 
Planetary models (with chemically distinct inner cores) 
7 Mass of Radius of Central 
Radius P : 
(km) inner core inner core pressure 
(10 g) (km) (10%* dyn/cm*) 

6 371 15°4 I 320 3°94 

6 367 15°3 I 315 3°82 

6 360 15°I I 310 3°67 

6 287 12'9 I 300 2°52 

6 285 12°8 I 290 2°47 

6 304 12°9 I 300 2°49 

5 797 9°2 I 200 1-38 

4 452 3°92 950 0°68 

2 473 0°62 540 0°173 

1 980 o'31 430 0'106 

1 267 0'07 270 0°042 


TaBLe II 
Planetary models (with inner cores formed by pressure) 
Mass Radius _Mass of Radius of Central 
25 inner core inner core pressure 
(10% g) _ (10% g) (km) (10! dyn/cm’) 
597°6 6 371 15°4 I 320 3°94 
592°3 6 382 6°75 I 025 3°63 
583°6 6 390 oh) ° 3°28 
498°8 6 356 2°00 
495°1 6 426 1°50 
498°5 6 453 1°33 
357°0 5 832 1:00 
152°1 4477 0°50 
23°9 2 485 0'130 
119 I 990 0080 
2°83 1 268 0°025 
Tables I and II show data for the various planetary models that were con- 
structed and also for the Earth Model B. The set of data in the first three lines 
of entries in each table is relevant to planets which would have inner cores on either 
hypothesis (the third line of Table II corresponds to a planet whose central 
pressure would just reach the pressure at the outside of the Earth’s inner core). 
The set of data in the next three lines shows how the assumptions underlying 
Table II give nearly the same mass (near that of Venus) under widely different 
central pressures. ‘This corresponds to Ramsey’s suggestion that Venus may have 
had an early catastrophe. 
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INTENSITY AND RADIAL-VELOCITY MEASUREMENTS ON THE 
SPECTRUM OF ZETA AURIGAE AT RECENT ECLIPSES* 


A. McKellar and R. M. Petrie 
(Received 1952 September 5) 


Summary 


Fifty-nine spectrograms of { Aurigae have been measured, 25 near the 
1950 eclipse, 29 near the 1947-8 eclipse, and 5 near earlier eclipses. 

The equivalent width of the chromospheric K-line of Call was measured 
on the 1950 series of plates. When these intensities were plotted against 
phase from mid-eclipse, the curves showing growth and decay of the line 
did not follow either of the two general courses found for previous eclipses, 
but fell between them. Furthermore, the curves for ingress and egress in 
1950 differed appreciably from one another. 

The radial velocities given by chromospheric lines of Call, of other 
metals, and of hydrogen, were measured, both on spectrograms taken near 
the 1950 eclipse and on suitable, available, earlier plates. The deviations 
of these radial velocities from the normal radial-velocity curve of the K-type 
star of the binary were plotted. It was found that the deviations were 
definitely different from eclipse to eclipse, and also that significant variations 
sometimes occurred from one day to another. Further, the behaviour of 
the deviations before and after eclipse was of a type not explicable by ordinary 
stellar rotation. In a few cases, the chromospheric K-line was found to be 
complex, the main line having a weak companion. 

It is concluded that the various effects, including variation in extent of 
chromosphere, irregular behaviour of the various chromospheric lines with 
respect to radial velocity and the existence of complex line structure, can be 
best explained by postulating the existence of varying prominence-like 
activity on a large scale in the extensive outer atmosphere of the giant K-type 
star. 





Introduction.—The observations and measurements discussed in this paper 
fall into two parts. First, measurements were made of the intensity of the 
chromospheric absorption lines of ionized calcium (A 3933 K-line) as produced in 
the spectrum of ¢ Aurigae before and after the total eclipse of 1950. Secondly, 
these results led to wave-length (radial-velocity) measurements on the K-line and 
other chromospheric lines on all suitable Victoria spectrograms and toa rediscussion 
of the available relevant data for recent eclipses of this interesting binary star. 

The eclipsing system ¢ Aurigae has been described many times during the past 
two decadesf, and it is necessary only to point out that spectrographic studies just 
preceding and following the geometrical eclipse of the relatively small hot B-type 
star by the large cool K-type star provide a powerful means of investigating the 
atmosphere of the latter. As will be apparent from the subsequent discussion, the 
general conclusion reached from our study is that the extensive outer chromosphere 
of the giant K-type star must be in a state that is far from static. The behaviour 
of the ionized calcium lines in the spectrum gives evidence of large-scale motions 

* Contributions from the Dominion Astrophysical Observatory, No. 29. 

t See, for example, reference (x) in which references to numerous previous papers are given. 
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of bodies of gas, possibly of the general nature of solar prominences, but on a much 
greater scale. ‘This conclusion is in accord with those reached in a short note on 
{ Aurigae by W. H. Christie (2) and in recent work (3) on the similar star, 31 Cygni. 

Information concerning the Victoria spectrograms used in the present work is 
given in Table I. The plates were obtained at the Cassegrainian focus of the 
73-inch telescope with the single-prism form of the universal spectrograph. All 
except two of the spectrograms were photographed with the camera of focal length 
71 cm (designated IM) and nearly all the rest with a camera of focal length 42 cm 
(designated 1S). ‘These two cameras gave linear dispersions at the K-line (A 3933) 
of 20 and 33 A/mm respectively. ‘The focal length of the f/18 collimator was 
114m; the spectrographic slit-width, usually 0-05 mm, was sometimes narrowed 
to 0-025mm. ‘The other two spectrograms, taken during total eclipse and used 
in the determination of the normal radial-velocity curve of the main component 
star, were obtained with the two-prism spectrograph and long-focus camera, which 
gives a dispersion of 8-5 A/mm at A3933. Before the 1947-8 eclipse the observa- 
tions obtained at Victoria near eclipse times were comparatively few and scattered. 
Hence, for comparison with our more recent data, earlier observations made 
elsewhere, including those of the 1934 eclipse by Christie and Wilson (4) by Beer 
(5), and by Guthnick, Schneller and Hachenberg (6) have been used. 

The phases from mid-totality given for the observations in Table I were 
calculated from the expression recently recommended by F. B. Wood (7), 
namely, primary minimum =JD 2 432553-666+972-162E. The length of 
totality is about 37-1 days, and of partial geometrical eclipse of the B star is 
about 1-3 days, hence the phases corresponding to second and third contact 
would be +18-6 days and to first and fourth contact would be +19°9 days (8). 
The nature of the “‘edge’’ of the disk of the K-type giant is undoubtedly not 
sharp; indeed the length of eclipse is different in different wave-lengths of 
light (7), so the figures in the preceding sentence are not specified with a precision 
better than one-tenth of a day. 

The intensity of the chromospheric K-line of tonized calcium.—The K-line of 
Call at A 3933 is a particularly favourable chromospheric line for study in the 
spectrum of { Aurigae and similar stars. Among the reasons are: (1) it is an 
ultimate line in the spectrum of singly ionized calcium, which is a comparatively 
abundant metal; (2) the very strong and wide K-line in the spectrum of the late- 
type star of the binary affords a window through which the presence of the chromo- 
spheric line on the effectively continuous spectrum of the B-type star may be 
detected and studied; and (3) ionized calcium is apparently supported at a very 
high level in the chromospheres of the stars (as it is in the case of the Sun). The 
other component of the ultimate doublet, the A 3968 H-line, is not as favourable for 
studies of line intensity, since it is immediately adjacent to the He line of hydrogen 
(A 3970) which is present in the spectra of both the K- and B-type stars of the binary. 

The method of measuring the total intensity of the chromospheric K-line on 
the spectrograms taken near the 1950 eclipse was that normally followed at this 
observatory. All the plates had been calibrated using a rapidly rotating step-sector 
and auxiliary spectrograph. Microphotometer tracings of the plates were made, 
and these were transformed, using the semi-automatic intensitometer, to log 
intensity versus wave-length tracings. On the logarithmic tracings the background 
intensity was drawn and the tracings then run through the intensitometer again to 
produce the final intensity versus wave-length profiles, having unit intensity (the 





TABLE I 
Radial Velocity of Various Lines and Intensity of Chromospheric K-Line in the Spectrum of € Aurigae near Eclipse 
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* These values of the intensity of the chromospheric K-line were measured on Victoria plates by 
H. L. Welsh. The intensities listed for the 1947-8 spectrograms are also taken from Welsh’s work (1). 

+ These are spectrograms upon which the K-line was noted and measured as complex; for details of 
these measurements see Table II. 
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continuum) a horizontal straight line. Unlike the situation usually met in the 
violet region of the spectra of late-type stars, in this case, for the K-line, there was no 
difficulty in locating the effective continuum. It was essentially, as previously 
noted, the continuous radiation of the B-type star in the A 3933 region, as seen 
through the wide deep K-line of the spectrum of the late-type giant. The latter 
line could be seen on each tracing as an extensive, apparently shallow dip in the 
composite spectrum. ‘This dip had the same depth and width on all the tracings. 
Only on plates very near totality, between first and second contacts, and between 
third and fourth contacts, when part of the disk of the B-type star would be cut off, 
would there be any difficulty in deriving directly the profile and equivalent 
width of the chromospheric line. In the present paper we are not concerned with 
the line at these phases, hence no such problem was encountered. The equivalent 
widths of the chromospheric lines were measured by obtaining their areas on the 
intensity tracings, using a planimeter. ‘These equivalent widths, expressed in 
angstrom units, are listed in Table I, column 5. 

Several investigators have made series of measurements of the equivalent width 
of the chromospheric K-line for the atmospheric phases of eclipses of { Aurigae 
between 1934 and 1948. H. L. Welsh (x) has given measurements for the 1947-8 
eclipse and, in addition, has summarized and discussed the previous results. For 
the ingress of the 1947-8 eclipse he found that the K-line was visible for at least 
fifty days before totality, whereas during egress the chromospheric K-line was 
present for only ten days. ‘These observations confirmed the conclusion of 
earlier investigators that the chromosphere of the giant star varies in extent. 
Calling the lesser extent “normal’’ and the greater extent “abnormal”, Welsh 
considered that all the observations up to 1948 could be described fairly satis- 
factorily by these two cases. When the equivalent width of the chromospheric 
K-line was plotted against phase from mid-eclipse, there resulted the two shapes of 
curve labelled “‘ normal” and “ abnormal”’ as shown in Fig. 1. Only two or three 
observations failed to fall near these curves. The available points fitted the 
“normal” curve for the ingress of 1939 and for the egresses of 1934, 1939-40 and 
1947-8. Points fell on the “abnormal”’ extended curve for the ingresses of 1937 
and 1947-8 and for the egress of 1932. 

When the equivalent widths of the chromospheric K-line for the ingress and 
egress of the 1950 eclipse are plotted in the above way (shown as filled circles in 
Fig. 1), it is found that the total intensities do not fall along either the “‘ normal” 
or “abnormal” lines. In both cases they define fairly definite curves between the 
two types previously encountered. Furthermore, the curves for ingress and egress 
are notthe same. For example, during ingress the equivalent width of the K-line 
was 0°35 A at 18 days befere totality began, while it had the same equivalent width 
only ten days after totality ended. The equality in the intensities on these 
particular days may be seen clearly in Plate 10, showing reproductions of a series 
of representative single-prism spectrograms. We conclude that for the 1950 
eclipse the plotted curve of the intensity of the chromospheric K-line against 
phase, both at ingress and egress, falls between the two previously recognized cases 
and that, for ingress and egress, the curves are appreciably different. 

Possible explanations advanced to account for the different extent of the 
atmosphere of the K-type component of { Aurigae have included: (1) regular or 
irregular pulsation of the star, (2) a permanent asymmetry in the atmosphere of a 
rotating star, and (3) prominence-like activity on a large scale. The previous data 
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indicated that, as regards ingress and egress, there was no apparent regularity in the 
difference in atmospheric extent. The present results for the 1950 eclipse go 
further and suggest that the irregularity of the effect includes many different extents 
and not only the two previously recognized. While present data do not entirely 
rule out explanations of types (1) and (2) above, yet taken in conjunction with the 
earlier material discussed briefly by Christie (2), with results to follow on complex 
structure in the chromospheric K-line for { Aurigae and particularly with recent 
observations on the rapidly changing structure of the chromospheric K-line in the 
spectrum of 31 Cygni (3), they support more strongly the conclusion that the effects 
observed arise from prominence-like activity on a vast scale. This subject will be 
discussed more fully later in this paper. 
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Fic. 1.—Equivalent widih of the chromospheric K-line of Ca II plotted against phase from mid-eclipse. 





The radial-velocity measurements.—Since the behaviour of the chromospheric 
K-line as regards variation in intensity with phase for the 1950 eclipse differed from 
that at previous eclipses, it was decided to measure the chromospheric lines for 
radial velocity to see if any clue to the source of the difference could be found. 
When the measurements of radial velocity were made, it was apparent that the 
variation of radial velocity with phase was not the same as that shown by the only 
previous extensive series known to us, namely, those of Christie and Wilson (4), 
of Beer (5), and of Guthnick, Schneller and Hachenberg (6), for the 1934 eclipse. 
This result, in turn, led us to measure for radial velocity the good set of spectro- 
grams obtained near the 1947-8 eclipse at Victoria, as well as the several other 
plates of earlier eclipses, with the results described below. 

The measured radial velocities and data concerning the spectrograms are given 
in Table I. The weights attached to the mean velocities were assigned at the time 
of measurement with due consideration being given to the dispersion. In the 
column headed Ionized Calcium, K was always measured and, whenever possible, 
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H was also measured; for hydrogen Hd was always measured and He included 
when possible. 


The lines measured for radial velocities are separated into four groups as follows : 
(a) The lines K and H of ionized calcium; 
(6) the lines He and H8 of hydrogen ; 


(c) a number of extremely sharp and prominent lines, mostly of neutral iron, 
which appear at phases just outside totality and which are referred to as chromo- 
spheric metals ; 

(d) anumber of lines, chiefly due to iron, in the composite spectrum, measured 
on plates obtained during totality and outside eclipse, to obtain the relevant part 
of the velocity curve of the K-type star. 


The first task was to determine as precisely as possible the part of the normal 
radial-velocity curve of the K-type component of the binary near eclipse. On 
plates taken during totality the K-type spectrum was measured in the usual way. 
In the composite spectrum the radial velocities of the K star were found from 
spectral lines with wave-lengths greater than A 4200. ‘To the violet of this the loss 
in contrast caused by the increasing dominance of the B-type spectrum makes the 
measurements less certain. ‘The velocity curve of the K star is essentially linear 
over the phase interval covered by our observations. It was decided therefore to 
adopt the slope of the velocity curve at eclipse, as given by Harper’s revised 
calculation of the orbital elements*, and to fix the position of the line on the 
velocity axis from our results. In this way we hoped to take advantage of the many 
earlier observations used to define the velocity curve and at the same time to avoid 
systematic error by using our measurements based upon recently revised and tested 
wave-lengths (9). It can be shown that the slope of the velocity curve at mid- 
eclipse may be computed conveniently from the formula 

< = SK {1 + (ze =) cos M + $e* cos 2M +e cos 3M + eas \, (1) 


where P=period, K =semi-amplitude, e=eccentricity and M=mean anomaly. 
This converges sufficiently rapidly for moderate eccentricities. A linear relation 
with the computed slope was fitted to our velocities from the nineteen plates 
measured for the purpose (shown in Table I, column 12) and gave for the velocity 
over the region near eclipse, 

V = +19°0+ 0°283 (t—ty), (2) 
where ¢—?, is the interval in days from mid-eclipse. ‘This line lies between 
Harper’s curve and that of Christie and Wilson (4) for the 1934 eclipse. On the 
average it is I-7 km/sec more positive than Harper’s curve and about 0-8 km/sec 
more negative than the Mount Wilson result. We have made this new deter- 
mination in order to have a firm basis for comparing our velocities from the calcium 
and hydrogen chromospheric lines with those of the K star. In what follows all 
comparisons are with respect to the line represented by (2) above. 

The radial velocities derived from K and H are shown in Fig. 2, for the different 
series of observations, plotted against phase from mid-eclipse. Individual plate 
velocities are plotted in all cases. Fig. 2 also contains a similar plot in which the 
radial velocities from the hydrogen chromospheric lines, He and Hé, are shown. 
In both cases the full straight line is the portion of the normal radial velocity of the 


* Unpublished. 
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K-type star, determined as described above. The broken lines in the lowest 
section of the two sides indicate the run of expected velocities if the deviations of 
the calcium lines after the 1950 eclipse were to be interpreted as shifts arising from 
the outer atmosphere rotating with constant angular velocity, i.e. a normal 
rotation effect. 
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Fic. 2.—Radial velocities from the chromospheric lines of ionized calcium and of hydrogen 
shown with relation to the radial-velocity curve of the K-type star. 


The results from the chromospheric metallic lines are not given in a diagram, 
They are few in number because of the relatively short duration of their appearance 
but the velocities should be very reliable because of the number and sharpness of 
the lines. The radial velocities from this group of lines are slightly positive with 
respect to the radial-velocity curve of the K star before totality, and slightly 
negative after totality. The average value of the difference is 2-5 km/sec, agreeing 
almost exactly in direction and amount with the behaviour found by Christie 
and Wilson, who attribute the small amount of rotational velocity exhibited by these 
lines to a blending with the K-type lines of the composite spectrum. ‘The chromo- 
spheric lines seen on our spectra are so narrow and intense that we find it difficult 
to believe that they are seriously affected by blending. If our opinion is correct 
the true rotational period of the K-type star is much longer than has been previously 
accepted, being about 4 000 days. 

The chromospheric lines K and H, and also He and Hé, are well defined, 
and one gains the impression that they can be measured very accurately. It is 
of interest to obtain a numerical estimate of the precision of the velocities from one 
or two lines since the interpretation of our plots must depend upon this accuracy. 

We have estimated the precision in two independent ways. First, acomparison 
of the velocities of the K star with the velocity curve gives the probable error, based 
upon 16 plates, of a single-prism plate of average weight as + 1°35km/sec. We 
then find the probable error of one line of unit weight, reduced to the dispersion at 
A3970, to be + 4°6km/sec. A second estimate is obtained from the differences 
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between repeated measures of a number of spectra. From 13 such pairs we find 
the probable error of a single line of unit weight to be + 3-3 km/sec. 

The first estimate contains the effect of instrumental and guiding errors as well 
as errors of measurement and is to be preferred over the latter, which measures only 
the errors of setting. On the other hand, it is considered that the calcium and 
hydrogen lines are more suited to accurate measurement than the general metallic- 
line spectrum of the K star, and that the proper value may lie between the two 
derived. We shall, however, adopt + 4-5 km/sec as a safe estimate of the probable 
error of a single line of unit weight in our discussions. It may be noted that a 
weight unity, on our system, is applied to a poor line and that the chromospheric 
lines of hydrogen and ionized calcium were assigned weights on the average of 
51 and 5°8 respectively. ‘Thus the velocities from the chromospheric lines of 
hydrogen and calcium, for an average plate, have probable errors of about 
+ 2km/sec. 

Results of the measurements of radial velocity —The most important aspects 
of the radial velocities appear to be, first, their departure from the symmetrical 
distribution to be expected on the supposition of a single rotation of the atmosphere 
of the K star and, secondly, the large amount of scatter in the individual points and 
the differences between eclipses and between the numerical results from hydrogen 
and ionized calcium. 

On the hypothesis that the K star rotates with a period constant at all distances 
from its centre, the radial velocities before first contact, and after third contact, 
should lie very closely upon a straight line intersecting the orbital velocity curve at 
the point corresponding to mid-eclipse and inclined to it at an angle depending upon 
the speed of rotation. At egress the radial velocities from K and H conform to this 
expectation in 1950 but show a departure in 1947-8. Furthermore, the velocities 
at ingress in 1950 depart widely from the line established by those at egress. Our 
observations, at ingress, at the earlier eclipses are too few to give a reliable result but 
they do not contradict the rotation explanation in the 1939 and the 1947-8 eclipses. 
In 1937, and again in 1950, the radial velocities from K and H lie below the orbital 
curve whereas they should lie above it. This behaviour is followed by the hydro- 
gen velocities, at ingress in 1950, but to a less marked degree. 

It is apparent from Fig. 2 that the hydrogen and calcium velocities agree at 
egress in 1950 but differ from one another at egress in 1947-8. On the hypothesis 
of a uniform angular velocity, the component of velocity in our line of sight 
through the K-type atmosphere would be constant for a given phase, regardless of 
the height to which any particular element extended. ‘The observations are not in 
agreement with this supposition ; the calcium velocities show the greater departure 
from the orbital velocity. 

It appears then that some influence, other than the rotation of the K star, is 
affecting the line-of-sight velocities of the chromospheric gases. ‘This makes the 
deduction of a velocity of rotation from the lines of calciumand hydrogen of doubtful 
meaning. However, if we were to accept the calcium velocity-displacements at 
egress in 1950 as indicating rotation (which we do not believe to be the case), the 
rotational velocity at phase 22 days would be 11 km/sec, and the graph of velocity- 
effect would be that shown by the broken line in the bottom sections in Fig. 2. 

Further information can be obtained about the nature of the motions by 
considering the scatter of the individual points in Fig. 2. Christie and Wilson 
remarked that the scatter they measured was “ much more than one would expect 
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from the character of the lines”, and we were at first inclined to the same opinion 
after measuring the plates. 

We considered first the velocities from spectrograms taken during a single 
night and have adopted, as explained above, a probable error of + 4-5 km/sec for a 
line of unit weight. A velocity-difference is considered to be significant only if it 
is at least three times its probable error. On this basis there is little evidence that 
single-night differences may be considered as real velocity changes. ‘Thus, for 
example, a range of 13 km/sec is shown by the calcium lines on JD 2 432 574 but 
the deduced probable error is + 4-8km/sec. Again, the hydrogen velocities at 
JD 2433552 show a variation of 8-8km/sec but the associated probable error is 
+2:8km/sec. When all the velocities are scrutinized, we find that the errors of 
observation can probably account for the variations between the radial velocities on 
one night, and thus a mean value for each date may be taken. 

A different situation exists when one considers the nightly means. ‘Thus, 
on the successive nights JD2432580-1, the calcium velocities give 
AV =7-5 + 2:0km/sec, a difference which is probably significant. ‘There appears 
to be an abrupt change at JD 2 432580; the mean velocity from six plates on 578 
and 580 is + 15:4 +0°8km/sec, while three nights later (mean of 581 and 583) we 
get +6-3+1-7km/sec from five plates. The difference, 9 + 1-9km/sec, appears 
to represent a real change. It is as if a calcium “ wind” were blowing steadily at 
11 km/sec up to 50 x 10°km above the photosphere and increasing then to reach a 
gale of 24km/sec at a height of 80x 10°km. These changes are not shared by 
He and Hé but between the nights J D 2 432 574 and 581 a difference of 12 + 2 km/sec 
is observed. In 1950 the changes are smaller but some significant variations are 
found. It should be recalled that the motions discussed are only those in the line 
of sight, i.e. tangential to the K-type star. 

Perhaps the best evidence for the existence of irregular motions comes from a 
comparison of the measures with the velocities expected assuming simple rotation. 
The calcium velocities at egress in 1947-8 give a dispersion about the rotational 
velocity of +8-7km/sec per unit weight, instead of the expected value of + 6-7; 
while the 1950 velocities give a normal value, +6-1. Evidently irregular dis- 
turbances of considerable magnitude are present in the velocities for 1947-8, which 
result in a dispersion of +5-6km/sec. Results of the same kind but with even 
greater dispersion are found from the hydrogen velocities both in 1947-8 and 1950. 

Structure in the calcium chromospheric lines.—An important feature of our 
observations is the detection of faint “‘satellite’”’ lines beside the main K absorption. 
These were seen on several spectrograms and were measured for radial velocity 
by both of us on all three plates of JD 2 432580. Asummary of the measurements 
for the cases where the structure was seen with certainty is given in Table II. One 
might think that we are observing the remnant of a line in the K-type spectrum but 
the appearance of this satellite line and its position one angstrom unit from the 
centre of K, where there is practically no light in the spectrum of the K-type star, 
argue against this explanation. Nor can this be an absorption line of ionized 
calcium in the normal spectrum of the B star; its radial velocity and short life 
dispose of this interpretation. We regard the “ satellites’”’ as strong evidence of 
the chaotic nature of the motions of the chromospheric gases. It is believed that 
this phenomenon may be rather general and that we have seen only the most 
spectacular manifestation of it, the less easily seen satellite lines being blended 
with the main absorption line to produce the irregularities in radial velocity 
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discussed above. ‘This conception receives support from the analogous and much 
more firmly established effects found recently in the spectrum of 31 Cygni (3) 
during atmospheric eclipse. It is of interest here to draw attention to a note 
by Christie (2) on the observations of faint absorption lines at K in spectra 
of ¢ Aurigae at phases far removed from the time of eclipse. Also, evidence of 
structure in the chromospheric K-line during the eclipse of 32 Cygni in 1949 was 
reported by D. B. McLaughlin (10) and has been supported by K. O. Wright (11). 


TABLE II 
Observed Cases of Structure in Chromospheric Calcium Lines* 





Radial Velocity 





Date Phase 
Main Line | Satellite Line 





km/sec km/sec 
1948 Jan. 30 +-26°925 +11°'0 — 98 
1948 Jan. 30 > OeF7 +15°5 — 61 
1948 Jan. 30 + +989 +16°1 — 84 
1950 Sept. 22 +21°053 +12°6 +105 




















*In addition to the four cases listed, a faint but definite line was measured about 
110 km/sec to the violet of the K-line on one of the plates obtained 1947 December 14, 
and on the three plates of 1950 September 21. These plates were all at phases from 
20°0 to 20°2 days, and showed well-developed chromospheric spectra. It is believed that 
the faint component was a chromospheric line, probably A3932-007 of Till. 


Conclusions.—Our general conclusion from both the intensity and radial- 
velocity studies is that the upper atmosphere of ¢ Aurigae contains moving clouds 
of atoms, perhaps similar to prominences, which markedly influence the intensities 
of the calcium chromospheric lines and dominate the radial velocities measured 
from the calcium and hydrogen chromospheric lines. ‘The evidences in support 
of this hypothesis are: 

(a) the varying rates of growth and decay of the chromospheric calcium lines 
at the various ingresses and egresses of different eclipses ; 

(6) the failure of the calcium and hydrogen velocities to give a symmetrical 
rotational effect before and after totality ; 

(c) the differences in radial velocity of chromospheric lines at different eclipses 
and between different elements at equal phases ; 

(d) the existence of sudden changes in velocity and the excessive scatter about 
a mean velocity when many good observations are available; and 

(e) the existence of satellite lines corresponding to velocities sufficient to show 
them separately with single-prism dispersion. 

Because of the variable behaviour of the radial velocities and intensities of the 
chromospheric lines, it is considered premature to attempt to draw more specific 
conclusions. Also, it appears highly desirable that observations of radial velocity 
and intensity similar to those described should be continued at future eclipses as a 
powerful tool for investigating the mechanical state of the chromosphere of the 
late-type giant star. Further, the comparative scarcity of spectrographic 
observations during ingress as compared to egress, at all eclipses studied up to the 
present, emphasizes the need of early and intensive observation preceding eclipse, 
in order to procure well-balanced data. 
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A COMPARISON OF THE INTENSITIES OF INFRA-RED AND 
VIOLET RADIATION FROM THE SOLAR CORONA AT 
THE ECLIPSE OF 1952 FEBRUARY 25 


D. E. Blackwell 
(Received 1952 September 24) 


Summary 


Measurements made in Khartoum at the total solar eclipse of 1952 
February 25 show a considerable excess of infra-red radiation in the solar 
corona at 2°5Ro from the solar centre. The ratio 

I(1 9m, aS Nea 15Ro) 

1(0°43p, 2°5Ro)/ 1(0°43p, 1'5Ro) 
equals 2-17, where I(p, q) is the intensity of the corona at the wave-length p and 
distance g from the Sun. The measurements were made with a 20-in. mirror 
of 45 in. focal length, using a lead sulphide cell and a photomultiplier. 
It is shown that the infra-red excess supports the current theory of an 
F-corona caused by diffraction of sunlight by a cloud of interplanetary dust 
particles. ‘Three models with widely differing distributions of dust are 
calculated but it is not possible to distinguish between them by using 
existing infra-red and visible data alone. 








Introduction.—It has been supposed by several authors (1, 2, 3, 4) that 
two scattering mechanisms contribute to the light of the solar corona. The two 
coronas that they form have been designated F and K, the former being due to 
diffraction by interplanetary dust particles and the latter to scattering by free 
electrons at a kinetic temperature of about 10° degrees. ‘The observations made 
in Khartoum at the eclipse of 1952 were designed to test this theory and provide 
further data for determining the size and distribution of the dust particles. 
They consisted of the measurement of the ratio of coronal brightness at 1-gu 
to that at 0-43 at as many positions as possible. 

Equipment.—Radiation from the corona was received by an equatorially 
mounted 20-in. reflecting telescope of 45 in. focal length, at the prime focus of 
which was the detector unit shown in Fig. 1. In this unit, radiation from a 
selected portion of the coronal image was focused directly on to a lead sulphide 
cell by a lens of infra-red transmitting glass. A part of the radiation was also 
reflected to one side by a partially aluminized glass mirror and focused on to the 
photomultiplier tube. Although about 15 per cent of the infra-red radiation and 
85 per cent of the ultra-violet radiation is lost, this arrangement has the great 
advantage that almost simultaneous measurements can be made at the two 
wave-lengths by switching from one detector to the other. The loss of ultra-violet 
radiation is of no consequence because the limits of the experiment are set by the 
amount of infra-red radiation available. The radiation entering the detector 
unit was interrupted at a frequency of about 750 cycles per second, and the 
signals from the two detectors amplified with a wide band-tuned amplifier, that 
from the lead sulphide cell passing first through a cathode follower and a pre- 
amplifier. The amplified signals were fed into a phase sensitive rectifier, the 
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auxiliary signal for which was obtained from the same chopping disk with a 
6-volt lamp and photoelectric cell. The output appeared as a meter reading with a 
response time of about rsec. Both the lead sulphide cell and the photomultiplier 
tube were used uncooled, but this did not matter because the signal/noise ratio 
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Fic. 1.—Optical arrangement of detector unit. 
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Fic. 2.—Infra-red spectral response of detector unit. 
Ordinates: Spectral response. 
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of this particular lead sulphide cell increases very little on cooling and there was 
ample ultra-violet light for the photomultiplier. A mica infra-red filter was 
placed before the lead sulphide cell and a Wratten No. 47 filter, with a trans- 
mission peak near 0-43, was before the photomultiplier tube. Fig. 2 shows 
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the spectral sensitivity of the apparatus to sunlight in the infra-red region, this 
curve being a combination of the spectral sensitivity of the cell, the transmission 
of the mica filter and the optical components of the detector unit, and the 
transmission of the atmosphere ; the effective wave-length was taken to be 1-9 p. 

In order to reduce errors which would arise from a non-linear response, the 
detectors were used in a balanced arrangement. The inner corona was measured 
first with a diaphragm over the main mirror. ‘This screen, which was pierced 
with a number of holes 9-5 mm in diameter and transmitted 11 per cent of the 
incident light, had been chosen so that a half-scale deflection would be obtained 
with the lead sulphide cell. Then the output from the photomultiplier was 
observed and the amplification of the signal from it altered until approximately 
the same deflection was obtained : this amplification was not changed again during 
the eclipse. In subsequent observations farther from the Sun the screen was 
removed. Extensive observations on the Moon before the eclipse had, in fact, 
shown the response of both detectors to be linear, but it was thought wise to use 
this extra precaution. 

The selecting aperture was 30 minutes of arc long by 5 minutes of arc broad 
and was maintained constant in size throughout the experiment. Such a large 
aperture was needed only for examination of the outer corona, but it was kept 
the same for the inner and outer corona because a change in size during the 
experiment might have caused a spurious result through variations of sensitivity 
over the receiving areas of the detectors. For convenience in setting, the aperture 
was placed perpendicular to the celestial equator on the E. side of the Sun. 

Observations.—At the time of the eclipse the sky was fairly transparent and 
almost free from cloud. Owing to the climate and the nature of the surrounding 
country, dust deposits formed very easily on the optical parts and were a great 
nuisance, but they have not interfered with the observations. 

It was not possible to make observations at many points during the three 
minutes of the eclipse. Sixteen readings are needed to make a reliable comparison 
between only two points of the corona—two for each wave-length in the two 
positions and eight zeros—and in addition to this a balance point must be found 
and the telescope directed in turn to two measured places in the corona. Also, 
if the points are not in an intense part of the corona a further eight readings are 
needed to measure the sky background. In this experiment thirty readings 
were taken at distances from the solar centre of up to 4:5Ro (all quoted distances 
are measured from the solar centre), but the last two sets of readings must be 
rejected. The last set, which was intended to give a measure of the sky brightness, 
is rejected because the sky was beginning to brighten when they were taken; 
the preceding set because they show a brightness only a little greater than that 
of the sky and are unreliable without an accurate knowledge of the sky background. 
The readings that are left compare the colour of the corona at 2:5Ro from the 
centre with the colour at distance 1°5Ro from the centre. The fact that no 
photoelectric measures of sky brightness are available is not of great importance 
because the corona at these places is many times brighter than the sky and an 
estimate only of the sky brightness is needed. An estimate of the brightness 
at 0°43 was obtained from a coronal photograph, taken by Dr von Kliiber, giving 
a value of 4 of the coronal brightness at 2-5Ro. The sky brightness at 1-9 
was assumed to be zero, an assumption supported by measurements made 
previously on the daylight sky and the region of the night sky near the full Moon. 
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If I(0-43 4, I-5.Ro) is the observed intensity of the corona, in arbitrary units, 
at a distance of 1-5Ro from the solar centre and at a wave-length of 0-434; then 
the mean of the two sets of observations gives 

U(t-94,2°5Ro) [(t-9u,15Ro) _ 

1(0°43 4, 2°5Ro)/ 1(0°43 #, 15Ro) 
the individual measurements being 2:20 and 2:13. This means that the corona 
is considerably reddened at increasing distance from the Sun. 

In general, one of the difficulties of photometric observation of the outer 
parts of the corona is to correct satisfactorily for light scattered by the instrument 
from the intense inner corona. Experiments on the full Moon in a clear sky 
before the eclipse showed that in this particular case such scattering was negligible 
at both wave-lengths, nor would much scattering be expected because there was 
only one optical surface and a diaphragm to exclude all but the radiation actually 
measured. If scattering of this kind had occurred it would have caused an 
excess of blue in the outer corona, not the observed great excess of infra-red. 

Another difficulty is the elimination of the atmospheric background, which 
is particularly difficult for the fainter parts of the corona. But such errors should 
not affect these observations appreciably, for if the intensity of the background 
brightness were in error by 10 per cent, the corresponding error in colour could 
be no more than I per cent. The magnitude of the infra-red excess is a 
sufficient argument against the supposition that it originates either in the Earth’s 
atmosphere or in instrumental scattering. 

Other observations.—These are the first coronal observations to be made with 
a lead sulphide cell, but other observations of a similar nature have been made 





= 2°17, 


with photoelectric cells, thermocouples and photographic plates. Abbot (5), 
Stetson and Coblentz (6) and Pettit and Nicholson (7) have all made measurements 
of the infra-red radiation from the parts of the corona near to the limb, with 
conflicting results. Allen (2) has made photographic observations of the 
variation of colour with distance from the limb. Allen finds 


1(0°65 u, 2°4Ro) /1(0-65 u, 1-2Ro) 

1(0°40 2, 2°4Ro)/ I(0°40 p, r-2Ro) 
Both Ludendorff (8) and Grotrian (g) find no change of colour with change of 
radial distance, but their measurements were over a smaller range of wave-lengths 
and radial distance. It will be shown later that these last three results are not 
incompatible with the present observations. 

Discussion of results.—Allen (2) and van de Hulst (3) have summarized the 
evidence for the existence of the F-corona and shown how it may be accounted 
for by diffraction of sunlight by interplanetary dust particles. The strongest 
evidence so far has been the existence of Fraunhofer lines in the coronal spectrum, 
and the falling off of the percentage polarization below the value predicted by 
a model consisting solely of free electrons. In both these discussions it is 
assumed that the F-corona can be identified with the zodiacal light, and this 
assumption is used in the reduction of the observations because otherwise there 
would be insufficient data to determine both the size of the diffracting particles 
and their spatial distribution. A more satisfactory proof of their identity would 
be the calculation of the zodiacal light intensity from data provided by observation 
of the corona alone, and it was originally hoped to achieve this, but there are still 
insufficient data. 
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In the following analysis we show how it is possible to account for the infra-red 
excess on the basis of diffraction of the continuous solar spectrum by interplanetary 
dust. We neglect the contribution from infra-red emission lines, but this seems 
reasonable because their contribution in the visible region is less than I per cent. 
We follow the method of calculating the diffracted light intensity which was 
originally given by Allen, with the refinement of a satisfactory correction for the 
angular size of the Sun andthe limb darkening. The latter correction is important 
because of the large range of wave-length, and because we shall be considering 
parts of the corona very near to the solar limb. 

According to Allen, the total amount of sunlight diffracted into unit solid 
angle by one particle of radius p is 


4 
p I 
1012 5 7B 
I+—s 


where /, is the surface brightness of the Sun expressed as the rate of radiation 
of energy by unit area into unit solid angle, Q is the solid angle subtended by the 
Sun at the particle, A is the wave-length of light and @ the angle between the 
directions of the observer and the Sun at the particle. 





Fic. 3. 


Hence the intensity of the light diffracted by a cloud of particles, of concen- 
tration N(L) per unit volume per unit solid angle viewed from the Earth, is 


rop*!, : oe .. dL 


_ » a ML) — ae ae 
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where L is the distance of a particle from the Earth, D is the distance between 
Sun and Earth, and s the radius of the Sun (see Fig. 3, which has been taken 
from Allen’s paper). The upper limit is chosen to be 4s because within this 
limit all particles are vaporized by solar radiation. 

Putting dL = —dr and 6=R/r this expression becomes 


romp‘ ,s? t N(r)rdr 
2 s 47p3R3’ 
r+ x 


where N(r) is the number of particles per unit volume at a distance r from the 
Sun, and R is the projected distance of the line of sight from a surface element 
of the Sun. To a good approximation, J, is the rate of radiation from the solar 
surface in the direction of the observer, so that it is necessary to make a correction 
for limb darkening as well as for the finite size of the Sun. Putting in these two 
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corrections we find that the radiation from the dust corona at a distance R from 
the centre of the Sun is 





, 


ef ° rNir)I,(s)drdA 
D 2, 47PR 
ae ial 
the first integration being taken over the projected area of the Sun. 
It is convenient to define an effective scattering coefficient og , which gives 
the brightness of the dust corona, at distance R from the solar centre, in terms of 
the mean brightness of the Sun. 


_ 1omp's? ( i rN(r)I,(s) dr dA 
°R, A= A? | D 47 47*p* RS 
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O's 1o rs 20 
Fic. 4.—Colour curves for model No. 1, particle size p= 10~* cm. 





Ordinates: Intensity ratio 0 3/09.45 y+ 
Abscissae: Wave-length in microns. 


opr, , has been evaluated graphically for a range of particle sizes, wave-lengths 
and values of R, taking N(r) to be a constant between r = 4s and r=D, with the 
results shown in Figs. 4, 5,6. In these curves the ratio 04/0943 has been plotted 
as a function of A for various values of R. This ratio is a measure of the colour 
of the corona in terms of the colour of the Sun, a value of unity meaning a colour 
equal to that of the Sun averaged over its apparent disk. In interpreting these 
diagrams we bear in mind that they refer to the dust corona only ; the real corona 
has an admixture of electron scattering for which o,/943, is equal to unity 
independent of wave-length and distance from the limb (neglecting a small effect 


47* 
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due to limb darkening). A striking feature of the diagrams is the strong 
dependence of colour, in the wave-length interval 0-43 .-I'9, upon the size of 
particle. With a particle radius of 10-*cm there is an infra-red excess at 3‘5Ro 
and an infra-red deficiency at 2-5Ro and closer to the Sun, although at wave- 
lengths shorter than 1-2, there is a slight infra-red excess at 2°5Ro. Passing 
to particles of size 2 x 10-*cm there is an overall increase of reddening, and at 
10-*cm radius there is a reddening at all distances from the Sun. The last 
diagram shows the way in which the curves for all distances tend to converge to 
a limiting curve as the particle size increases. Fig. 7 demonstrates the superiority 
of I-9p radiation over 0-65 u radiation for the determination of particle size. 

+o 





T T ' 35Ro 


MODEL NOI 
PARTICLE SIZE p= 2x10 CM 


j t i i 
os ro 's zo 
Fic. 5.—Colour curves for model No. 1, particle size p=2X 10~* cm, 











Ordinates: Intensity ratio 03/0943 u- 
Abscissae: Wave-length in microns. 


As we have reliable infra-red data for only two places in the corona it is not 
possible to calculate the particle size without recourse to other data. The other 
data that have to be used are determinations of Fraunhofer line intensities in the 
coronal spectrum, for these give directly the relative proportions of dust diffraction 
and electron scattering at short wave-lengths, and a measurement of the surface 
brightness of the corona in terms of the brightness of the Sun. We also make 
the initial assumption that the radiation from the inner corona, distant 1-5Ro 
from the solar centre, is identical with solar radiation. This is a reasonable 
assumption because the Fraunhofer lines are weak in comparison with the 
continuous spectrum here, and the reddening of the dust corona near to the solar 
limb is slight, 
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Let o, and a, be total effective scattering coefficients at red and blue wave- 
lengths respectively. Each scattering coefficient is the sum of two scattering 
coefficients, o, due to dust and o, due to free electrons. 


o Cga-ete 
Thus — = htt 
Cy Sgn ten 

O,g"eto 
= £5" (because o, ,=0, ») 
Santee 


=2:17 (by observation). 


(=), -27/(Gt -217) 


Fan t+,» 
og la, gt 
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Og ,=1'17 
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Fic. 6.—Colour curves for model No. 1, particle size p=107* om. 





Ordinates: Intensity ratio oj) 09.43 4 
Abscissae: Wave-length in microns. 


oq,» at a distance of 2-5Ro from the solar centre can be calculated in terms of 
the particle concentration from the expression already given, (o4 »+¢,, ,) is the 
brightness of the corona at this point in terms of the brightness of the solar disk— 
according to Allen this is 0-078 x 10-*—and a, ,/a4 , can be read off for various 
particle sizes from Fig. 7. In this way there can be constructed the relation 
between particle size and concentration which is shown in Fig. 8, a curve which 
ends at the particle size p= 2-3 x 10-*cm because particles smaller than this size 
give a reddening at 2°5Ro of less than 2:17. This does not give the distribution 
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of particle sizes in the interplanetary dust, for so far these have been assumed to 
be of one size only, but if the particle size can be fixed it will give the concentration. 
The particle size can be determined from the measurements by Allen and Grotrian 


of the Fraunhofer line intensities, these giving directly (. a) from which 
d e/ b 

@q -/%4,, and then p can be calculated. The results are given in Table I. The 

p, N curve slopes steeply in this region so that the large uncertainty in the intensity 

of the Fraunhofer lines does not cause a large difference between the two estimates 

of the particle radius. Indeed, as will be clear later, all the uncertainty arises 

from the choice of model and not in the choice of supplementary data. 
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2 + 6 8 10 x io CM 
Fic. 7.—Comparative sensitivities of 1°9 4 and 0°65 radiation to particle size. 
Ordinates: Intensity ratio 
A. (0y-94/o-asule-sRo 
B. [09-65 4/%0-asule-8RO- 
Abscissae: Radius of particles. 


Tasie I 
Observer ( <¢ ) p N 
Fate] b 
Grotrian o-71 2°78 X 107% cm 15°6 X 1071*/cm! 
Allen 0-61 3°07 X 107* cm 12°0 X 107!*/cm! 


As this corona model depends upon only one infra-red comparison it is 
desirable to check it with other observations. Suitable observations are those of 
Allen, Grotrian and Ludendorff on the variation of corona colour with distance 
from the Sun, those of Allen on the variation of the Fraunhofer line intensity 
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intensities of infra-red and violet coronal radiation 


with wave-length, and those of Ohman and Allen and others on the variation of 
polarization with wave-length. Results of the calculations are given in Table II, 


together with the observations. 


TaB.e II 





Radius 
of particles 


[-65u/%o-su]2'4RO 
[ Fo-45u/Fo-a0u] 1'2R 0) 


Measured by 
Allen: 1°48 


Ratio of 
Fraunhofer 
line intensities 
at o-6041/0-40u 
at 1°5RO 
Measured by 
Allen: 1-o 


[%0-65/o-sau] 
at 1°35RO 
Measured by 

Grotrian : 1-0 


| 





Polarization at 4600A 
Polarization at 6250A 
at 2-0oRO 
Measured by 
Allen: 1°o 





3°07 X 10-° cm 
2°78 X 10-8 cm 











1°30 
1°28 


\ 1°03 
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; 
5 6x10 CM 


Fic. 8.—Permitted relation between concentration of particles and radius of particles. 


The agreement with observation is, at first sight, fair. 


Ordinates: Number of particles/cm® x 10", 
Abscissae: Radius of particles. 


However, Allen’s 


measurement at 2°4Ro is almost certainly too high because he obtains a 16 per 
cent variation of colour between 1-2Ro and 1-6Ro, even though at 1-6Ro the 


F-corona contributes only 10 per cent to the combined F- and K-coronas. 


It is 


not easy to make a comparison with Allen’s measurement in the third column as 


he has averaged the intensity ratios over all distances from the Sun. 


In making 
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the comparison we have supposed that this is equivalent to a measurement at 
1°5Ro. Although the calculated result is 30 per cent higher than the measure- 
ment, it is fair to say that a result of 1-3 would fit Allen’s observations as well as, 
or better than, a result of 1-0 (see Fig. 7 of Allen’s paper). The fourth column 
shows that the colour effect to be expected is very small near the limb and that such 
measurements can give little information about the F-corona. 





25 T T T 


MODEL NO. 5 
PARTICLE . SIZE p= 0.25xI0 CM 
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Fic. 9.—Colour curves for model No. 5, particle size p==0°52 X 107° cm. 
Ordinates: Intensity ratio 04.94/7-43 
Abscissae: Wave-length in microns. 


It is also possible to obtain an estimate of particle size from the decrease in 
polarization of the corona below that calculated from the pure electron scattering 
model, assuming that the F-corona is unpolarized. It is well known that the 
results do not agree well, and van de Hulst has already given a discussion of 
possible causes for the discrepancy. But the fact that there is no detectable 
change of polarization with wave-length (4, 10) has to be explained on the basis 
of our model. Some effect would be expected because a dust corona with a 
particle size of 3 x 10-*cmisstrongly reddened. Theeffect can be easily calculated 
from the data of Figs. 4-6 and knowledge of measured polarization, with the 
result that for a distance of 1Ro from the limb and a wave-length range of 4600 A 
to 6250A (corresponding to the experiments of Allen in 1940) there is a variation 
of no more than 4 per cent in the polarization. This is smaller than the scatter 
of all polarization measurements so far made, and it would be difficult to measure 
it with certainty. The change of polarization with wave-length increases with 
increasing distance from the limb, becoming 6 per cent at 2-5Ro from the solar 
centre, but already at this distance there are great difficulties in allowing for the 
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polarization of the background radiation—which must really be the background 
radiation and not simply the dust corona. 

Other corona models.—This is a naive model and we should neither complain 
if it does not fit the facts nor be content if it does. A completely general model 
would have a distribution of particle sizes, and the space distribution of each 
size would vary in a specified way with distance from the Sun. As more data 
accumulate it may be possible to elaborate the model, but at the moment we can 
only calculate other models to see if they agree with the scanty data better or 
whether they can be rejected by reason of complete disagreement. Two other 
models have been examined ; model 3 in which the particles are nearer the Earth, 
extending with uniform concentration from distance %D to distance D from the 
Sun, and model 5 in which the particles are near the Sun, extending from distance 
4s to ;.D from the Sun. The representative curves of Figs. 9, 10 show that 
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Fic. 10.—Colour curves for model No. 5, particle size p= 10~* cm. 











Ordinates: Intensity ratio 4.9 p/% 9-43 
Abscissae: Wave-length in microns. 


to obtain the required reddening a larger size is needed for particles near the 
Earth than for particles near the Sun. The calculated data for these two models 
are given in Table III. 

It is clearly impossible to distinguish between these widely differing models 
by appealing to data relating to the inner part of the corona and a short wave-length 
region, although this comparison with observation makes it likely that there is 
an increase of concentration towards the Sun rather than a decrease. We lack 
sufficient information at this stage to discuss models having a range of particle 
sizes. 
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Further experiments.—The problem is largely solved if by an independent 
method we can measure the effective particle size without appealing to zodiacal 
light observations. ‘The curves given show that observations near the Sun using 
visible light can contribute little more. But from simultaneous measurements 
of the infra-red excess at three or more points and the polarization in the infra-red, 
which gives directly the ratio o, ,/(o4,,+,, ,), both the particle size and the spatial 
distribution may be deduced. ‘The latter measurement can be made with an 
infra-red quarter-wave plate, infra-red polaroid and a lead sulphide cell. The 
problem is simplified in the intermediate portion of the corona, extending for 
several degrees beyond 6Ro from the solar centre, because this is a dust diffraction 
pattern without the complication of the electron corona. However, the inter- 
pretation is not easy because as it is necessary to integrate over a large range of 
diffraction angles, the rate of decrease of intensity is not sensitive to the choice 
of particle size. 

Tase III 





Ratio of Fraunhofer 


Model} Observer 


Radius 
of particles 


[0-05u/70-40u]2"4RO 
[0-65 ,/ F040] I -2R © 


[20-65u/%o-say] 
at 1°35RO 


line intensities at 
0-6op/0-40u 


(cm) Measured by 


Grotrian : 1°0 


at 15RO 
Measured by 
Allen: 1-o 


Measured by 
Allen: 1°48 





Allen 
Grotrian 
Allen 
Grotrian 


4°04 X 107% 
3°74 X 107% 
0°31 X 107% 
0°28 x 107% 


1°45 
1°39 
1°38 
1°34 


1°35 
1°32 
1°25 
1°22 
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FURTHER OBSEPVATIONS OF MAGNITUDE WITH A 
COARSE GRATING 


R. v. d. R. Woolley, K. Gottlieb and A. Przybylski 
(Received 1952 February 6) 


Sumunary 


Further determinations of magnitudes centred on 5980A _ and 
4480 A with the coarse grating attached to the 30-inch reflector on Mount 
Stromlo are given. The relation between the Mount Stromlo scale and 
the International scale is discussed and the position of the Sun considered. 
The colour-luminosity diagram for some of the stars is shown. 





1. Further observations have been secured with the coarse grating attached 
to the 30-inch reflector on Mount Stromlo. ‘These have used the technique 
already described *, except that the plate is now driven relative to the telescope, 
which goes at normal speed, whereas formerly the telescope was trailed slightly. 
The present arrangement is therefore a one-dimensional schraffierkassette with 
dispersion perpendicular to the trailing motion. ‘The arrangement gave good 
results on slow-moving stars near the (South) pole when the star was so bright 
that a fairly short exposure was possible. During longer exposures the telescope 
wandered enough to cause trouble, and defects in the images occurred frequently 
in near-equatorial stars. ‘They were attributed to want of uniformity in the 
rate of movement of the telescope. Since rate and not mean position is concerned, 
it is very difficult to cure these defects by guiding. Some images were rejected 
on account of this fault. 

The central images and spectra of the stars (dispersion 9800 A/mm in the 
first order) were weakened by a rotating sector just in front of the photographic 
plate. Usually two sectors were employed such as would weaken the star to 
about 6™-o and 6™-5, if only one trail were made: but where fainter stars were 
measured the trails were repeated, up to thirty times for the faintest star measured. 
The plate was driven through 4 mm in thirty seconds and then reversed. Each 
plate was calibrated by exposing it in the calibrating spectrograph as already 
described (except that the plate is now driven instead of the pinholes); but the 
method depends for its accuracy on the angles of the sectors used to weaken the 
stars rather than on the plate calibration. Each plate contains a number of 
rows of spectra of the same star (usually weakened by at least two different 
sectors) obtained at one setting of the telescope. The plates usually begin and 
end with a row of spectra of a standard star (i.e. one observed photoelectrically 
on Mt. Stromlo). ‘The standard stars are, however, too bright for comparison 
with the faintest stars, as the heaviest sector only cuts down the light by five 
magnitudes. In measuring faint stars it was therefore necessary to use an 
intermediary comparison star. ‘The stars used for this purpose are indicated 
in the tables of results. 

* R. v. d. R. Woolley and K. Gottlieb, M.N., 110, 206, 1950. 
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No attempt was made to determine absorption. Instead, stars were observed, 
as far as possible, with 1-2<sec z<1-4, and corrections of o™-1 A sec z in the 
red and o™-3 A sec z in the blue were applied. 

The spectra were measured, as previously described, at 5980 A and at 4480 A. 
They were reduced by taking as standard the magnitudes of the stars whose 
names are shown in heavy type in Table I. The blue standard magnitudes 
were adopted from the work of Hogg and Mrs Hall *, applying a small correction 
according to the relation 


m(4550) — m(4480) = 0-061(m(5980) ~ m(4480)), 
to transfer the photoelectric results from 4550A (used by Hogg and Mrs Hall) 
to 4480 A. 

The red standard magnitudes were found from a least squares solution of 
the measured differences between the pairs used as blue standards. 

To fix the zero of the red magnitude scale, gradients were calculated (with 
a preliminary red zero) from G= — 1-645 (m(5980) — m(4480)) for eighteen stars 
whose gradients had been measured by Gascoigne.t The zero of the red 
magnitude scale was adjusted so as to make the mean difference between the 
gradients and Gascoigne’s gradients vanish. 

Table I gives results for 67 stars, including the 16 standard stars. Successive 
columns give (1) the HD number, (2) the BS number, (3) the star’s name, 
and (4) and (5) the coordinates for 1950-0. ‘These are followed by (6) the 
magnitude at 5980 A and (7) the magnitude at 4480 A. Column (8) shows 
the spectral type. If a determination of spectral type and absolute magnitude 
by Miss Hoffleit f is available, it is shown first, and followed by a spectral type 
and class by Miss Woods § where this exists. If neither of these is available the 
HD spectral type is shown. 

Column (9) shows the gradient G from G= —1-645{m(5980) —m(4480)} 
and column (10) shows the photovisual magnitude m,,, which is interpolated 
between m(5980) and m(4480) to 5430A according to the interval in 1/A, 
so that 

M yy = M5430) = m(5980) + 0-184 G. 

Table II shows magnitudes and gradients of 43 stars in the Harvard 
regions F1, F2 and F3. The comparison stars used are shown. ‘Table III 
shows results for 23 stars in the Hyades, compared with « Hydri, 8 Hydri and 
B Octantis (and, on two plates, with 8 Carinae). 

2. Order of accuracy of the results.—Results were only kept if a star was 
obtained on at least three plates and if the plates were in reasonable agreement. 
The weight attached to a determination on a single plate is 2/{1/m+1/n}, where 
m and n are the number of spectra of the star and comparison star respectively. 

The internal consistency of the observations can be judged from the least 
squares solution for the red magnitudes of the sixteen comparison stars. The 
average (unweighted mean) of the standard errors from the least squares solution 
is +0™-05, 

Turning to the blue magnitudes, differences between the adopted magnitudes 
{based on the photoelectric work of Hogg and Mrs Hall) were compared with 

* A. R. Hogg and Beryl Hall, M.N., 111, 325, 1951. 
+ S. C. B. Gascoigne, M.N., 110, 15, 1950. 


} D. Hoffleit, Annals of Harvard College Observatory, 105, No. 3; 119, No. 1. 
§ Miss Woods, C.S.O. Memoirs, to be published. 
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Further observations of magnitude with a coarse grating 


TABLE I 
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18 866 
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43 834 
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78 045 
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104 902 
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109 026 
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119 164 
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141 891 
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B Cari 
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Crux 
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18 200 
18 259 
18 384 
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18 522 
a Circ 
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6°54 
3°99 
5°37 
2°52 
4°50 
4°66 
0°44 
2°59 
4°77 





7°72 
4°58 
5°37 
3°22 
6°11 
5°23 
0°26 
3°02 
4°96 
5°69 
4°73 
4°47 


6:00 


4°56 
4°21 
4°29 
5°13 
401 
3°38 





Gs 
Go, 4°4; GoV 
Bog 
G1, 3°5; GrIVv 
K2, —o'1 
F7, 3:0; F6V 
Bs IV 
Fo, 2:1; Fo V 
As, 0-7; Ag III 
F2 
K2, 0-2 
Mo, o-o; Ma II-III 
Ko 
Gg, o’o 
Ko 
Go, 4°3 
F7, 0°6 
Bs V 
F3, 2°8 
K1, 08 
2°0; As Ill 
—o8; Ks I 
Ao, III 
1°7; FolIb 
B8 
K6, o-o 
Bg III 
K3, ov1 
K2, 0°7 
Bs IV 
B3 V 
Bs IV 
B3 IV 
B3 V 
Bo III 
Ao 
Go 
F8 
F8, 3'9 
Bs 
Ko 
Ko 
F2, 2°4 
Ks, 0o*0 
G8, 1°1 
AoV 
Ki, o°5 
Fi, 2:1; F2V 
G2, —1°0; G2Il 
Ms, —1°1 
Go, 4°0 
Gg, o'o 
Ko, 0*4 
Ks, —o-4; Ka III 
B8 V 
Ki, 0°3 
Ao 
Bg IV 

















Woolley, K. Gottlieb and A. Przybylski 
TaBLe I—Continued 


Vol. 112 





HD 


2 
BS 


3 
Star 


4 


“1950-0 


5 


81950. 0 


6 


Mso80 


7 


mM 4480 


8 
Spectr. 





171 759 


6 982 


¢ Pavo 


18h 37m 


—71° 28’ 


3°65 


4°89 


Ko, 0°5 


A Pavo 
€ Pavo 
p Pavo 
p*® Pavo 
§ Pavo 
¢' Pavo 
B Pavo 
a Pavo 
a Octn 
o Pavo 
y Pavo 
v Octn 
a Grus 
% Octn 
«a Tucn 
v Indi 
8 Tucn 
Grus 
Octn 
GC 33 061 
y' Octn 
y*® Octn 
7 Tucn 
e Tucn 
6 Octn 


18 
18 
19 
19 
20 


173 948 
188 228 
188 584 
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201 371 
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211 416 
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212 581 
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214 846 
223 444 
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224 362 
224 392 
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7 603 
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7 665 
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twenty measured differences, of average weight 6-1. ‘The root mean square of 
the measured difference minus adopted difference was +0™-04. 

The gradients of eighteen stars can be compared with Gascoigne’s gradients. 
In these the mean difference has been adjusted to be nearly zero (actually 
— 0°00, + 0°02) by a choice of zero for the red magnitudes. The root mean 
square of the individual gradient minus Gascoigne’s gradient is +0-08. Since 
G=1-645(Am,—Am,), an error of +0°08 in the gradient corresponds to 
+0™-05 in the magnitude difference or +0™-03, in either magnitude if 
independent. It is likely however that the red and blue magnitudes are not 
quite independent and that the standard error of a comparison between two 
stars with the coarse grating method, of weight 6, is about +0™-05 in either 
colour. 

The systematic errors of the method, judged by comparison with Eggen’s 
work on the Hyades, are much smaller. Since the magnitudes m(4550) by 
Hogg and Mrs Hall are related to Eggen’s Pg, and C, by 


m(4550) = Pg, —0-16 C, 
and since m(4480) = m(4550) + 0-037 G, the relation to be tested is 


m(4480) — 0-037 G= Pg, —0-16C,,. 


We have in fact, from 19 stars in Table III (excluding those of type Ko), 


mean (m(4480) — 0-037 G) — mean (Pg, —0-16C,) = +0™-00, + 0-01; s.e., 
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TABLE II 


3 4 5 6 7 
1950-0 019500 ™5980 M4420 Spectr. 


I 
HD 


F Region 


Comparison stars : 


2151 
12311 
21024 
24 512 
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29 769 
21 166 
24 636 
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25 191 
25 224 
25 269 
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TaBLe [I1—Continued 
4 5 6 
8950-0 ™s5980 M4480 


F3 Region 
Comparison stars : 


158 094 17526™ —60°39’ 3°55 = 345 
188 228 18 54 —73 3 3°81 3°82 
190 248 20 3 —66 19 3°22 4°02 
197 O51 20 40 —66 23 3°28 3°46 
214 846 22 41 —81 39 §= 3393 4°19 


196 051 20 35 —76 22 5°73 6-22 
186 117 19 43 “32 T 7°24 
188 229 19 55 —74 9 6-10 7°26 
186 502 19 45 —72 55 6:98 7°49 
187 915 19 53 —74 16 7°20 7°71 
184 997 19 38 —74 35 7°28 7°70 
189 899 20 3 —74 22 7°21 7°83 
189 723 20 2 —74 54 681 7°99 
186 219 19 43 —72 38 5°25 5°48 


HD 196051=BS 7 863 =p! Octn. 
HD 186 219=BS 7 498. 
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meso 4480 Yerkes - Eggen 
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a highly satisfactory result. The comparison is of course made through the 
connection by Hogg and Mrs Hall of their system with Eggen’s. The four 
stars of type Ko give m(4480) — Pg, = —0™-24 + 0™-02. 

The relation between Eggen’s colour C, and the Greenwich gradient G is 
not linear, but is nearly so between G= +0:1 and G= +1-0, when C,=0°41 G 

a oe 
ee a m,,+.0°387 G = Pv, + 0°84C,, 
and with Cs — G, Mov + 0°043 G= Pov,. 
Applying this result to the nineteen Hyades stars (excluding the K stars) 
mean (m,, +0°043 G)— mean Pv, = —0™-01, + O™-01,, 

again a very satisfactory result, in view of the rather long chain of comparisons 
involved. 

HM—mpv 





+ 








| | | | l | m 
0-00 0-50 1-00 1-50 2-00 2-50 3-00 G 
Fic. 1.—Comparison with Harvard Mimeograms. 
+ Stars from Table I. @ Stars from Table II. @ Stars from Table III, 





3. Comparison with other results. (a) Harvard magnitudes.—One hundred 
of the stars observed (60 in Table I, 17 in Table II, and 23 in Table III) have 
magnitudes assigned to them in Harvard Mimeograms. Differences HM—®m,, 
are shown plotted against G in Fig. 1. Excluding stars for which G>2-25, the 
results for 77 stars in Tables I and II are represented by 

HM-—m,,= +0™-22—0:16 G, 
the standard deviation for individual stars being +0™-12. For the 23 stars in 
Table III (the Hyades), 

HM—m,, = +0™-17—0-12G, 
with a standard deviation of +0™-09. 

* O. J. Eggen, Ap. F., 113, 664, 1951. 
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Thirty-one stars in Table II have photovisual and photographic magnitudes 
assigned to them. Mean residuals and standard deviations from them are shown 
in Table IV. 

In forming Table IV, the stars L, M and O in F1, and D in F2, which show 
large residuals in the red, were omitted from the photovisual results (but 
retained in the photographic). According to the Mt. Stromlo measures the 
Harvard magnitudes make the stars too faint in the red, especially in region F2. 
This is reflected in the negative Harvard colour indices shown in this region, 
even for stars of type F5. 

(b) Holmberg’s magnitudes of Hyades stars.—Fourteen of the stars in ‘Table IV 
have values of JPv and JPg listed by Holmberg.* The values shown are the 
results of photographic comparisons with the North Polar Sequences. ‘The 
relation between Holmberg’s JPv and m,, is 


IPv—m,, = +024 + 0™-03 
and Holmberg’s colour index C, is related to the Stromlo gradient by 
C,=0°6 G—0-07, 
the standard deviation of individual comparisons from the latter relation 
being +0°II. 
TaBLe IV 
Mean Residuals in F Regions 
Harvard minus Mt. Stromlo 
Region Photographic Photovisual 
H—m(4480) 8.e. H—my,» s.e. 
m m m m 
Fi +0°'07+0°'07 G +0°17 +0°37—0'16G +0-09 
F2 +o'13+0:°07G +0°'12 +o'58—-0'16G +016 
F3 +0:'20+0°'07 G +o'1r +0°28—0'16 G +012 
(c) Southern bright stars measured by Cousins.—T welve of the stars in Table I 
appear in Cousins’ f list of photovisual magnitudes. Comparison shows 
Cousins — m,, = 0™-04(m — 4:7) — 0°01 G+0™-19, 
the standard deviation being +0™-10, The magnitude term in this relation 
agrees with the magnitude term shown by Cousins in comparing his magnitudes 
with those of Williams.f 
4. Position of the Sun on the Mt. Stromlo and IPv scales.—According to 
Woolley and Gascoigne§ the magnitude difference between the Sun and Sirius, 
interpolated to 5430A, is 25™-54+0™-04. Hogg and Mrs Hall give —1™-49 
for the magnitude of Sirius at 4550A. Assuming that the gradient of Sirius is 
zero, the m,, of Sirius is also —1™-49, and then for the Sun 
My, = —27™-03, M,,= +450. 
The connection between HM and JPo is 1Pu=HM+0™-10+0™-06C. We can 
make use of the relations between HM and m,, given in the preceding section, 
to find a relation between JPv and m,,. It will be sufficiently accurate to take 
C =0°'5G, as the term in C is small: and then 
m,, =1 Pv —0™-32 + 0™-13 G from stars in Tables I and II, 
= 1 Pv—0™-27+0™-09 G from stars in the Hyades. 
E. Holmberg, Medd. Lund Obs., Ser. I1, No. 113, 1944. 
A. W. J. Cousins, M.N., 103, 154, 1943. 
E. G. Williams, M.N., 102, 226, 1942. 
R. v. d. R. Woolley and S. C. B. Gascoigne, M.N., 108, 441, 1948. 
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Taking G= +1-1 for the Sun, we have for stars of this colour m,, =1Pv—o™-18 
from the first relation and m,,=IPv—o™-17 from the second. Preferring the 
result from the Hyades, we have for the Sun 


I Pv (apparent) = — 26™-86, I Pv (absolute) = + 4™-73. 


If we work through Holmberg’s magnitudes (Section 3) the results are —26™-79 
and +4™-80 respectively. Kuiper* gives —26™-84 and +4™-75 from a 
discussion of a number of observations. The agreement between all of these 
results is as good as could be expected. 


G 
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Fic. 2.—Relation between gradient and spectral type. 


Hoffleit O Dwarf Woods @ Fi © 
@ Giant Harvard © F290 
© Early-type stars Hyades x F3 @ 


5. The connections between gradient, spectral type and luminosity.—The 
gradients of all the stars observed are plotted against spectral type in Fig. 2. 
Unfortunately the spectral types are not on a uniform system. For stars in 
Table I, Miss Hoffleit’s spectral types are used, failing them, Miss Woods’: and 
if neither is available the spectral type is taken from the Henry Draper 
Catalogue. ‘This was also used for the stars in Table II. For stars in the Hyades, 
the Yerkes spectral type (as shown in Eggen’s paper) t was used. 


* G. Kuiper, Ap. 7., 88,°429, 1938. 
t,.O. J. Eggen, Ap. F., 111, 65, 1950. 
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The pecked line shows the relation between gradients of unreddened stars 
and Yerkes spectral type, according to Gascoigne.* It may be noted that 
Miss Hoffleit’s classifications of F-type stars are usually later than Miss Woods’, 
and that if the latter had been used for F-type stars the points representing 
F stars would lie closer to Gascoigne’s line. 

Some of the early-type stars are conspicuously redder than their spectral type 
would suggest, notably HD 25 423, 25 224, 28 457 in the Fr region and HD 106 759 
and 104 685 in the F2 region. 

A notable feature of the gradient-spectral type diagram is the wide dispersion 
in gradient of stars whose types are close to Ko. This is partly due to an 
absolute magnitude effect. 

Mov 0 1-00 ? . G 


~~ 





Sub- giant 


Owarf 
Sub-dwarf 








1 i at I f 


—0-08 40-44 40-90 +10 +125 Cy (Eggen) 
Fic. 3.—Colour luminosity array—Hoffleit’s luminosities. 





The gradient-luminosity relation is probably more significant than the 
relation between gradient and spectral type. It is shown in Fig. 3, in which 
Miss Hoffleit’s luminosities are used for the stars in Table I, and Eggen’st 
moduli for stars in the Hyades are adopted. The pecked lines show Eggen’s 
relations for various sequences, and Fig. 3 confirms these generally without 
exhibiting the fine structure of his classifications of early-type stars. The late- 
type giants shown in the diagram establish reasonably well the relation 


M,,.= —1™-12(G — 2-35) 





*S. C. B. Gascoigne, M.N., 110, 15, 1950. 
t O. J. Eggen, Ap. F., 112, 141, 1950. 
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between absolute photovisual magnitude and gradient. The root mean square 
deviation in magnitude is +0™-45. Schlesinger* states that the average 
probable error of spectroscopic absolute magnitudes is about + 0™-4. Accordingly, 
the relation between absolute magnitude and gradient in the subgiants may 
well be of the close kind found by Eggen for other classes of stars. 

Absolute magnitudes and parallaxes may be found for some late-type stars 
in Tables I and II by assuming them to be giants and to obey 

M,, = — 1™-12(G — 2°35). 

The likelihood of these stars not being giants must of course be considered. 
There is not much risk that they are dwarfs, because dwarfs have conspicuously 
smaller gradients than giants in the spectral types considered. For example, 
HD 22676, Ko, G=1-45 and HD1o01805, Go, G=0-94 are both probably 
dwarfs and have not been included in Table V. There remains the risk that 
some of the stars are subgiants. But late-type subgiants are rare. Miss Hoffleit 
does not give any among the stars in Table I and Eggen comments on the rarity, 
which amounts to complete absence in types later than K2. The possibility 
remains, however, that when spectral classes are available, one or two of the 
stars in Table V will be found to be subgiants and therefore about two magnitudes 
fainter than is shown in the table. 


TABLE V 


Absolute magnitudes and colour parallaxes of late-type stars 
Star Star 
HD Sp. G My» p HD . G Myo p 


” m 


m ” 
323 1°94 “tT 0°46 104752 2°14 + 0°23 0°005 7 


23 128 1°88 +o'§2 108 558 2°52 —O'IQ 00023 
24940 2°66 —0°35 117 000 ; 2°04 +0°35  0'0059 
25 269 2°84 —0'55 119 076 2°12 +0'26 0:0048 
28 778 y 189 +0°51 119 164 c 2°25 +O'II 0°0039 
33 285 1°71 +o-71 188 229 < 1'9g! +049 0°:0064 
37 763 ‘ 1'79 «= +: 0°63 188 584 < 1°84 +0°57 0°0095 
98 627 3°09 — 0°83 188 887 2°27 +0°09 0'0093 
103 925 < 21906 +018 189 723 1°94 +0'46 0°0045 
104585 K2 2:07. +0°31 223 444 2°50 —O'I7 0'0042 


Acknowldgements.—Mr R. H. Brockman assisted in securing the observations 
and Miss Janet Lamb measured all of the plates (many of which were measured 
independently by the authors). The authors are indebted to Dr Olin J. Eggen 
for many very valuable discussions about magnitudes and colours of stars. 


Commonwealth Observatory, 
Mount Stromlo, Canberra : 


1952 January 23. 


* F. Schlesinger, Yale Cat. Parallaxes. 





MICROMETRICAL MEASURES OF DOUBLE STARS 
F. Schmeidler 


(Communicated by the Director, Cambridge Observatories) 


(Received 1952 September 1) 


During the author’s stay at the Cambridge Observatories some micrometric 
measures of double stars were made in the spring of 1951 with the 25-inch Newall 
refractor. The following list contains only 40 measures of 13 pairs; extremely 
bad weather conditions made a greater number of measures impossible. 

On some nights a mechanical defect which could not readily be rectified 
caused sticking of the wires of the bifilar micrometer. The defect is quite obvious 
and observations were always stopped when it was noticed, but it is possible that 
it was occasionally overlooked, so that smaller weight should be accorded to the 
measures of distance than to the position angles. 

For each pair in the following list of observations the first line gives the ADS- 
number, alternative designation, coordinates for 1950-0 and the brightness 
of the components; the subsequent lines give for each observation the date, 
sidereal time, position angle, distance, position of the eyes and quality of the image. 
For the position of the eyes, .. means eyes parallel, and : perpendicular to the line 
joining the stars. The quality of the images is given on an arbitrary scale, I 
meaning extremely good conditions and 4 very poor. Distances of only a few 
tenths of a second of arc are estimated. The last line gives the arithmetic means 
of all measures of the pair. 


The author owes acknowledgment to Professor W. Rabe for encouraging him 
to publish these measurements, and also for putting at his disposal computed 
places for those pairs for which orbits exist. He is also indebted to Professor 
Redman for permission to use the Newall telescope. 


6650 X1196 8hooM-3  +17° 48’ A 5™-o0, B5™-7, C 5M-5 
h ° 


” 


AB 1951°160 9°2 32°1 18 
284 8:8 32°9 19 
305 98 35°4 "12 
308 9°5 31°3 ‘10 
51°26 32°9 Irs 
O-C according to Makemson, A.7., 42, 153: +2°°6, 0”:00 
h ° ” 
AC 1951°160 94 g2°6 5°98 
284 9'0 92°7 6-01 
305 9°8 93°4 6:23 
308 9°6 94°3 6:06 
51°26 93°2 6:07 
O-C according to Makemson, A.7., 42, 153: +3°°8, +0”:40 
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7284 X3121 


1951°160 
395 
51°23 


O-C according to Van den Bos, Union Observatory Circular 99 : 


7307 21338 
19515157 
284 

393 

308 
51°26 


Ho 369 


1951°160 

393 

395 
51°26 


OX 210 
1951°303 


308 
51°31 


A 2145 


1951°*160 

395 

308 
51°26 


x 224 
1951°160 


395 
51°23 


O-C according to Ekenberg, Lund Medd., Series II, Nr. 116: 


8094 X1517 


1951 ‘160 
308 
51°23 


O-C according to Arend, A.7., 45, 63 : 


8140 X1534 


1951°281 

303 

395 

308 
51°30 


8148 x1536 


1951°305 


O-C according to Rabe (unpublished) : 
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g'14™-9 +28° 47’ 
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10°3 213°6 
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gh 49m 2 
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95°6 

, 98 59™-4 
11°3 260°8 


10°4 261°5 
261°2 


10% 06™-6 
h ° 


10°7 25°6 
10°9 40°2 
10°7 32°5 

32°8 


10h 37m +1 


209 
228°6 
218°8 


rr yyM.-y 
h ° 
r1°6 203 
11'5 194°7 
198°8 


11 ygM.2 
h ° 
11°9 321°2 
12°3 321°3 
I2°1 3211 
11°6 323°4 
321°8 


11h 27M.3 


h 9 
11°6 246°4 
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+4°°O, 


0°58 


0°58 


gh 17M-q +38" 24’ 
‘ * 


1°39 
1°39 
1°42 
1°40 


+ 36° 43° 
elongated 
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+ 46° 36’ 
1'22 
I'12 
I°l7 
+20° 35’ 
O'1s 
o'2 
O'1s 
O'17 
4 9° 06’ 


elongated 
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+-20° 25’ 
O'ls 
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5°66 
5°61 
5°25 
5°85 
5°59 
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8355 OL241 r1Ms3m-7 +35 44" Os my 
h ° ° 


1951°281 12°3 136-2 Ig! 
303 12°7 1351 1°57 

3°95 12°3 139°5 1°95 

308 12°2 139°7 1°72 
51°30 137°6 1°79 


8446 1606 12% 08™-3 = +40° 10’ 
h ° ” 
1951°303 12°9 307°! 0°73 
305 12°5 307°6 080 
308 12°5 3090°5 0°76 
51°31 308°1 0°76 


8539 21639 r2221M-g = +25°52’ 6m -7 7-9 
h ° ” 
1951°305 12°8 330°5 2a -: 4 
308 12°8 332°8 3 er 3-4 
51°31 334°6 1°33 2n 
O-C according to Schrutka v. Rechtenstamm, A.N., 279, 81 (orbit A): +4°:3, +07:16 


Grossfriedrichsburgerstrasse 21, 
Munich, Germany: 
1952 August 28. 
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A PHOTOMETRIC SURVEY OF SOLAR FLARES, PLAGES AND 
PROMINENCES IN Ha LIGHT* 


M. A. Ellison 


(Communicated by the Astronomer Royal for Scotland) 


(Received 1952 November 18) 


Summary 


Accurate Ha line profiles have been obtained for flares, plages and 
prominences of different types. ‘The instrumental equipment used for 
photographing the spectra consists of a combined spectrohelioscope and 
spectrograph which was designed for this purpose and incorporates a 
54-inch plane grating. 

Special attention has been given to the problem of eliminating sources 
of systematic error (a) by measurement of the scattered heterochromatic 
light within the spectrograph; (b) by measurement of the total intensity 
of the grating ghosts; and (c) by determination of the instrumental profile. 
As a check upon the effectiveness of these methods and corrections, the 
profile of Fraunhofer Ha has been re-determined for the centre of the disk, 
giving a corrected central intensity of 0-158, in good agreement with the 
earlier measures of Thackeray who used. a monochromator. 

The profiles of intense flares are discussed and estimates are given for 
the total radiation rate of Ha light in erg/s for five flares near their times of 
peak intensity. A combination of the photographic profiles with visual 
measures of central intensity and line-width yields the average values and 
the approximate ranges of these quantities at peak intensity for flares of 
importance I, 2, 3 and 3+. 

Eleven profiles of representative plages give central intensities ranging 
from 0°27 to 0°45 with a mean of 0°35, expressed as a fraction of the local 
continuum outside the plage. 

Ha profiles have been obtained for quiescent prominences by photo- 
graphing the same prominence at the limb and on the disk. ‘The following 
results are of interest: (a) the corrected limb profiles yield kinetic 
temperatures in the range 10000-20000 deg. K, in agreement with the earlier 
results obtained at Edinburgh by Dr M. Conway (now Mrs H. A. Briick); 
(6) some prominences yield double-peaked or flat-topped profiles for the 
light which has been emitted from their central and brightest parts. It is 
considered that this effect arises from self-absorption and that it also 
indicates a much higher temperature in the central regions of these prominences 
than in the outer layers; (c) a prominence of normal brightness on the limb 
(central intensity ~o-10) has approximately the same brightness when 
seen against the disk. 

The central intensities of the brilliant surge prominences are found to 
lie within the range 0:15—0-60, again expressed as a fraction of the adjacent 
continuum. 


* The full text of this paper appears in Publications of the Royal Observatory, Edinburgh, 1, 
No..5, 1952- 
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ERRATUM 
M.N., 112, No. 5: A. E. Roy, Models for unmixed stars. 


P. 491, Fig. 1.—Please substitute the figure reproduced below for the earlier 
reproduction, from which the numerals on the right-hand ordinate were omitted. 
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Fic. 1.—Values of the radius, luminosity and central temperature of a star with Kramers’ opacity, 
in terms of the Cowling continuous model, at various values of the discontinuity in molecular weight. 
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